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Abstract 



Models with extra dimensions have attracted much interest recently. The rea- 
son is that string theory, a serious candidate for the theory of unification of all 
interactions, predicts the existence of up to seven extra dimensions. These mod- 
els have many attractive features and may turn out to provide the solution for 
many long standing problems in physics. One interesting and very attractive 
idea is that our visible universe is confined to a four- dimensional hypersurface 
in a higher-dimensional spacetime. This membrane like universe was dubbed 
brane- world. 

The main goal of this thesis is the study of the four-dimensional (4D) effective 
theories and their observational consequences in the brane-world universe. 

After introducing the brane-world idea with more detail we shall use the gra- 
dient expansion method to obtain the 4D effective theories of gravity for several 
higher-dimensional theories with different numbers of extra-dimensions. 

In chapter [21 we shall consider a five-dimensional two brane model with a 
scalar field living in the extra-dimension. This model reduces to the Randall 
and Sundrum model and to the Hofava-Witten theory for a particular choice of 
our parameters. We will show that the 4D effective theory can reproduce exact 
dynamically unstable solutions present in the literature and we will identify the 
origin of this instability. We will argue against a claim that the 4D effective 
theory allows a wider class of solutions than the 5D theory. 

In the following chapter, we shall derive the low energy effective theory on 
a brane in 6D supergravity. In order to consider matter other than tension we 
will regularize the singular brane by a codimension one brane. We will show that 
the 4D effective theory is a Brans-Dicke (BD) theory with the BD parameter 
uiBD = 1/2. We will prove that time- dependent solutions in the 4D theory can 
be lifted up to exact 6D solutions found in the literature. We will also discuss 
the possibility of stabilizing the BD field in order to recover standard cosmology. 
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At the end of chapter [3] we comment on the imphcations of these models for the 
cosmological constant problem. 

In chapter HI we will present a systematic way to derive the 4D effective 
theory for warped compactifications with fluxes and branes in the lOD type IIB 
supergravity. We will obtain the 4D effective potential for the universal Kahler 
modulus. 

In the second half of the thesis, after introducing the concept of brane- inflation 
we will focus on some observational consequences of these low energy effective 
theories. In particular, chapter [6] is devoted to the study of non-Gaussianities 
in general single field models of inflation. We will compute the fourth order 
action for scalar and second order tensor perturbations for a fairly general model 
that includes the brane-inflation scenarios as particular cases. This will enable 
us to calculate the trispectrum at leading order in the slow-roll expansion. We 
will point out that in order to obtain the correct leading order result using the 
comoving gauge action, one has to include the second order tensor perturbations. 
This fact has been wrongly ignored in the literature. We will also provide the 
necessary formalism to calculate the next-to-leading order contribution that might 
become observationally important in the future. 

The final chapter before the conclusion studies the bispectrum in general mul- 
tiple field models of inflation. These general models include multi-field /T-inflation 
and multi-field Dirac-Born-Infeld inflation as particular cases. We shall derive the 
exact second and third order actions including the metric perturbations. In the 
small sound speed limit and at leading order in the slow-roll expansion, we shall 
derive the three point function for the curvature perturbation which depends on 
both adiabatic and entropy perturbations. We will show that the contribution 
from the entropy perturbation has a different shape dependence from the adi- 
abatic perturbation contribution if the sound speeds of these perturbations are 
different. This provides a way to distinguish models. 

The final chapter is for conclusions and discussions. 
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Chapter 1 

Introduction to the brane-world 



Superstring theory, or better, the union of the five different string theories known 
as M-Theory is regarded as a promising candidate for the "Theory of Everything" . 
Scientists hope that, once completely formulated, this theory will not only be able 
to reproduce, in certain limits, the two theories that are the pillars of the 20^^ 
century physics, quantum mechanics and general relativity, but it will also be able 
to explain new phenomena. One of the most interesting and surprising aspects 
of this theory is the fact that it can only be correctly formulated in a higher- 
dimensional spacetime, to be more specific a 10 + 1 dimensional spacetime. But 
it is also true that almost all of our daily experiences can be perfectly explained 
using just a four- dimensional spacetime, 3 + 1 dimensions. Therefore we need 
a mechanism of compactification of these seven extra dimensions, so that they 
become "invisible" at least at low energy scales. 

In 1996, Horava and Witten [2] came up with a model where this M-Theory 
compactified in a certain way would reduce to a known string theory [Eg x 
Heterotic String). Basically, this model can be realized, from an effective point of 
view [3], as a five- dimensional (5D) spacetime, the bulk, with two boundary hyper 
surfaces, the branes. All the Standard Model particles and fields are confined to 
these hyper surfaces while gravity is the only interaction that can occur in the 
bulk. 

In 1999, motivated by these ideas, Randall and Sundrum [1] proved that in 
this scenario one of the biggest problems in physics, the hierarchy problem could 
be explained in a simple way. I shall describe their model in section 11.11 In 
the same year, Randall and Sundrum published another paper [5] where they 
proved that Newton's law of gravity could be retrieved in these models with 
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extra dimensions. Although this was a very important paper, it was somehow 
incomplete because the recovery of Newton's theory is not enough when we know 
that the correct theory is general relativity. The gap was filled with the work of 
Garriga and Tanaka [6] (see also [7j) and the subject of linearized gravity will be 
discussed is section [L2l 

The effective non-linear Einstein equations of the 3-brane world were derived 
in [S] and they will be under consideration in section II. 3[ I shall say a few 
words about the cosmological implications [9l [10] of this new set of equations 
in subsection 11.3.11 The effective Einstein equations are very difficult to solve 
analytically. In fact, it was shown that in general one needs to solve the full 5D 
bulk equations in order to find the 4D geometry, so it is crucial to develop tools to 
tackle this problem. Kanno and Soda [TTl [T2] developed such a tool and section 
11.41 will be dedicated to the study of their method, the low energy expansion 
scheme. This perturbative method solves the full 5D equations of motion using 
an approximation and after imposing the junction conditions, one obtains the 4D 
effective equations of motion. 

In the next three chapters of this thesis we will use the low energy expansion 
scheme extensively in order to obtain the four- dimensional (4D) effective theories 
of gravity for three different models with several extra-dimensions. 

The first model under consideration will be a 5D dilatonic two brane model 
that interpolates between the Randall and Sundrum model and the Hofava- 
Witten theory, which we shall introduce in section II. 5[ 

The second model that we will study is the 6D supergravity scenario. As 
an introduction, in section 11.61 we will present the idea of self-tuning in a six- 
dimensional brane- world model (co-dimension two brane- world) . This self-tuning 
property is remarkable and it might prove to be the solution to another one of the 
biggest problems in physics, the cosmological constant problem. We will obtain 
the 4D effective theory and we will show that it is a Brans-Dicke (BD) theory 
with the BD parameter equals to one half. 

It would be desirable to realize the Randall and Sundrum mechanism of gen- 
erating large hierarchies in a more fundamental theory like string theory for ex- 
ample. This was achieved in [13] for the type IIB string theory. In section [TTTl we 
will present the lOD type IIB supergravity model that is a low energy approx- 
imation to the type IIB string theory. Section 11.71 is of particular relevance for 
the understanding of Chapter HJ where we will obtain the 4D effective theory in 
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warped compactifications of the lOD supergravity model. This constitutes the 
last of the three higher- dimensional theories that we will discuss. 

In this thesis, we shall follow the sign conventions and definitions of the Rie- 
mann and Ricci tensors of Wald [H] . In particular, we will use the metric signa- 
ture — !- + + ■■■ We also choose to work in an unit system where the speed of 
light and the reduced Planck constant are equal to one. 

1.1 Randall and Sundrum two brane model 

This model [1] was proposed as an explanation to the hierarchy problem. So, 
first I shall clarify what exactly is this problem and why it is one of the most 
important problems in physics. 

The hierarchy problem highlights the enormous difference in strength of grav- 
ity and the electroweak interaction. The electroweak energy scale is rriEw ~ 
lO^GeV while the Planck energy scale (the scale at which general relativity breaks 
down and a theory of quantum gravity must be used) is MpL = G^^"^ ~ lO^^Ge^. 
On denotes the 4D Newton's constant. It is problematic to have two (so differ- 
ent) fundamental energy scales when searching for the theory of unification of all 
interactions. 

However, there is a crucial difference between the electroweak scale and the 
Planck scale. While we have experiments that probe the electroweak interaction 
down to m^Jy ~ distances, we do not have experiments to test gravity 

at distances of order of Mp^ ~ 10~^^cm. The smallest scale that we have ac- 
cess through gravity is of order 10~^cm [15], so our assumption that Mp^ is a 
fundamental energy scale relies on a huge extrapolation. 

This reason led Randall and Sundrum (RS) and others before them to infer 
that MpL might not be a fundamental scale and to proposed mechanisms to 
generate it. Before RS, the higher-dimensional mechanisms [H] which existed 
solved the hierarchy problem relying on the formula 

Mj,^ = M"+2v;. (1.1) 

This formula is valid for a 4-|-n dimensional spacetime, where the n extra dimen- 
sions are compactified, i.e., the full spacetime is equal to a AD spacetime times a 
nD compact space, being Vn its volume. M is the {A + n)D Planck mass and we 
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Figure 1.1: The brane-world universe 



suppose M ~ Mew- If Jt- = 2 then ~ Imm. 

The ingenious step in RS work [1] was to consider a non-factorizable metric 
(11. 2p instead of a f actor izable one: 

ds^ = QABdx^dx^ = e-^^y^/^7]f,^dx^'dx'' + dy^, (1.2) 

where is of order of M and y is the coordinate for the extra dimension, 
^ y ^ L. rj^jj denotes the Minkowski metric. 

Their set-up is represented in figure 11.11 It consists of two 3-branes at the 
fixed points, y = (the hidden brane. A) and y = L (the visible brane, B), of the 
orbifold S^/Z2. Effectively the branes are the boundaries of the 5D spacetime 
(bulk). 

The action is 




where / is the curvature radius of AdS^, ayis, Chid are the brane tensions and 
is = Sh^^^Gn = M~^. denotes the 5D Ricci scalar, ^Jenotes the 

induced metric on the brane. 

They have proved that (11. 2p is a solution of the Einstein's equations that 
result from the action (11. 3p . without the matter terms, and obtained the following 
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relations for the brane tensions: 




(1.4) 



G 



It can also be shown that for a negative tension brane observer 



Mj,L = MH [e^^'^ - l] . 



(1.5) 



So if we take M ~ mEw ^.s a fundamental scale of the "hidden" brane (positive 
tension) we can easily (without introducing other large hierarchies) obtain a scale 
for gravity in the "visible" brane (negative tension) of the order lO^^GeV. This is 
achieved thanks to the exponential "warp" factor a{y) = e~^/'. 

1.2 Linearized gravity 

Once we have seen this new solution for the hierarchy problem it is relevant to 
discuss the characteristics of gravity in the RS model. RS also realized that 
their model could provide us with an alternative to compactification, i.e., we can 
reproduce Newtonian gravity even with an infinite extra-dimension. 

Let us begin with the same set-up of the previous section but with the negative 
tension brane {y = L) removed. We can do this simply by allowing L +oo. 
Now, an observer living on the brane will calculate the Planck mass as 



As we can see the Planck mass is finite in spite of an infinite extra-dimension. 
This is in contrast with the factorizable metric's models where the Planck mass 
would be infinite. As mentioned before, it is the fact that the extra-dimension 
is curved that acts as an effective compactification. It could be said that the 5D 
graviton only "sees" an extra-dimension of size /. 

In order to determine if we recover the Newtonian potential we need to study 
the behaviour of tensor perturbation /i^j, around the metric (11.21) . 



ds^ = [qab + hAB) dx'^dx'' = (e-^IJ'l/'r/^, + h^,) dx^dx" + dy\ (1.7) 



Ml 



PL — 



(1.6) 



With the gauge choice h 



yy 



h 



0, h 



V 



and h^^ = 0, the Einstein's 
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equations at first order in the perturbation are [H [6] : 

/ia/3 = 0. (1.^ 



We may look for a solution of the form h{x'^,y) = ipmiy)^'^^'^^'" (dropping 
tensor indices for simplicity), with = —rn^. The resulting equation for ipm{y) 
can be written in the form 

-d^ - V{z)i.m = -m V{z) = ^1^^ - (1-9) 

where the new variable is 2; = J a^^iv) dy and with the definition tjjm = a~^^'^iprn- 
The problem now reduces to a trivial quantum mechanics' exercise. A particle 
(the 5D graviton) governed by a Schrodinger like equation in a "volcano" shaped 
potential. The eigenvalue m = corresponds the eigenfunction ipoiu) = ^~!jf~- 
This is the "wave function" of the massless graviton of the 4D theory and we 
will see that it reproduces Newton's potential energy between two masses on the 
brane. The other eigenmodes (continuum) are plane waves, asymptotically away 
from the brane, and near the brane their amplitudes are suppressed. Because of 
this they are weakly coupled to the matter on the brane but they will introduce 
small corrections to the standard 4D physics. 

Outside a spherically symmetric source of mass Mq on the brane, the gravi- 
tational potential is given by 

3r' 



V{r)^^^(l + — ]. (1.10) 



This formula is valid in the limit r ^ / and the second term is the first correction 
due to the continuum of Kaluza-Klein (KK) modes. 

When r <^ I, the small scale limit, the potential reveals its pure 5D character 

T/^ N GnMqI /I 1 1 n 

V{r) f» — — — . (1.11) 

The best gravitational experiments available today place an upper bound of 
10~^mm for /. Formulas (11. 4p and (11.61) imply the inequalities 

a > 10^^GeV\ M > W^GeV. (1.12) 
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So, we should not expect to find KK massive modes of tlie graviton in future 
collider experiments. 

In [6] , it was proved, for the one brane case, that we do recover the linearized 
Einstein's equations and not only Newton's theory. When we introduce a second 
brane (negative tension), the KK modes no longer form a continuum but rather 
have a discrete set of values. At linear order, we get a scalar-tensor theory as an 
effective theory on the brane. The physics on both branes is quite different, for 
example, the Newton's constant is 

where (±) refers to the sign of the branes' tension. 

This scalar-tensor theory is a Brans-Dicke (BD) theory with the BD parameter 
given by 



{±) _ 3 (^±2L/l 
BD — 



2 (e^^"/' - 1) • (1-14) 

Observations require ubd ^ ^ [Ej- This is easily achieved on the positive 
tension brane and we have a perfectly acceptable theory of gravity even without 
stabilizing the BD field. In the negative tension brane (where we have an expla- 
nation for the hierarchy problem) the BD parameter is always negative so this 
cannot be our universe unless we stabilize the distance between the two branes. 



1.3 The Einstein's equations of the brane- world 

In the previous section, we have discussed linearized gravity in the brane-world. 
However for example for cosmology, we need to know the non-linear /strong be- 
haviour of gravity. Here we shall obtain the full non-linear projected equations 
for gravity in the brane-world. In the preceding sections, we have used Einstein's 
equations on the brane-world context but we have always assumed particular 
symmetries, in the metric or in the energy momentum tensor. In this section, 
following [8], we shall derive the equations that govern the brane-world dynam- 
ics without assuming any symmetries. We will proceed in this manner as far as 
possible and only then we will use some quite general brane-world's hypotheses. 

For any given 5-dimensional spacetime {V, g^j^y) and any (well behaved) 4- 
dimensional hypersurface (M, q^y) embedded in this spacetime, it is possible to 
write the components of the Riemann tensor of this 4D hypersurface as (Gauss' 
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equation) 



:i.i5) 



where = gixp—n^rti, is the projected metric on the surface, is the unit normal 
vector to M and K^j^y is the extrinsic curvature tensor given by K^j^y = q^^q^V a^p. 
The Codacci's equation reads 



DyKl - D,K = ^'^R.^n'^qP, 



:i.i6) 



where is the covariant derivative with respect to q^u- 
Gauss' equation imphes 



1 



-qflu 



^'^Rp. - ^gpa^'^R) qX + ^'^Rp^n'^n'^qpu 



+ KK^y - K^K,, - -q^y {K^ - K^^^K^p) - E^,, (1.17) 



where 



77 _ (5) pa nPnf^n'^ 

The 5-dimensional Einstein's equations are 



1 

T 



:i.i8) 



:i.i9) 



where T^^ is the 5D energy momentum tensor. Using the decomposition of the 
Riemann's tensor 



J^pauP — ^ {yp[u ^I3]a ya[u J^I3]p) ^yp[vyi3]a + 

where is the Weyl's tensor, we obtain 



pavf3i 



:i.2o) 



2k'q 



pv 



TpM + Tp^nPn'^ - -T q,. 



--q,y {K' - K'^^K^p) - E^ 



pvi 



:i.2ii 



where 



E^y = (^)C 



%pa'^an''q^ql, 



:i.22) 
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is traceless. The Codacci's equation gives 

D,K; - D^K = K^T.^n'^g^. (1.23) 

Equation (11.211) and (I1.23P are completely general and can be used in all cases. 
In order to proceed, we shall assume the following brane-world's symmetries: 

• Gaussian normal coordinates, ds^ = dy"^ + q^ydx^^dx" , where the brane is at 

y=o, 

• T^i, = —Ag^i, + S^^i,5{y), where A is the cosmological constant in the bulk 
and 5*^1, is the distributional energy momentum tensor on the brane, 

• Sfj,y = —Xq^u + T^u, where A is the brane's tension and r^^ is the energy 
momentum tensor with the property r^^'n}' = 0. 

We then impose Israel's junction conditions [18] 



M\y=o = 0, (1.24a) 
[K,,]ly=o = -^l(^S,y-lq,uSy (1.24b) 

The previous equations together with (11.2ip . assuming Z2 symmetry around the 
brane, give us the effective 4D equations that replace the standard Einstein's 
equations in the brane- world context, 

^^^G^^ = -A^q^y + SvrGTvT^,^ + k%ti^„ - E^^, (1.25) 



where 



A4 = jtSfA + i^A'), (1.26a) 

= (1.26b) 

111 1 

+ — rr^^ + -q^uTapr'^^ - ^^a^^^^- (1-26c) 



Important comments about these equations: 

• When the effective cosmological constant of the brane A4 is zero, it cor- 
responds to the RS fine-tuning between the brane's tension and the bulk 
cosmological constant. 
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From Eq. (]1.26bp it is possible to see that A has to be positive, 



The TT^jy tensor is quadratic in r and it will be important at high energies, 

There is the new E^^, term. It is zero for an AdS bulk but not zero otherwise. 
Some, but not all, of its degrees of freedom are constrained by the brane's 
matter distribution as D^^E^^ = k^D^tt^j,. These non-determined degrees 
of freedom are what transforms the system in a non-closed system, i.e., we 
need to solve the full 5-D Einstein's equations and not only the 4-D effective 
ones in order to evaluate E^u, 

A dimensional analysis of (11.25^ suggests that in the low energy limit, i.e. 
when the brane's tension is much larger then the characteristic energy scale 
of the matter, it simplifies to ^^^Gj^i^, ~ —A^q^^ + SttG nt^i^ (Einstein's equa- 
tion). 

From fll.23p and ( 11.24bl) it is possible to show that energy conservation on 
the brane still holds D^rj^ = 0. 



1.3.1 Cosmological scenarios 



As natural, there is a myriad of cosmological models that have been suggested as 
models for our universe. For a good overview of these models see [T9l [20l [2T] and 
references therein. 

Because we are interested in cosmological solutions (isotropic and homoge- 
neous) we will consider the Friedman-Robertson- Walker metric 



dS^^-^ = q^^dx^'dx" = -dt^ + a^{t) 



dr^ 



1 — fcr^ 



:i.27) 



where a{t) is the scale factor of the 3-brane. We assume the perfect fluid energy 
momentum tensor form for t^i, = pt^ty + Ph^'^, where is the fluid's 4-velocity 
and h^'^ = q^^ + t^^f^ is the metric of the spatial surfaces orthogonal to t^. p and 
P are the energy density and pressure respectively. This implies that ir^'^ is 



71' 



111/ 



1 

12 



p [ptH" + (p + 2P) h^""] 



:i.28) 



and its trace is vr = | (p + 3P). A general form for E^y respecting isotropy and 
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homogeneity is 

Ei^ = -diag{-pE,PE,PE,PE). (1.29) 
Efj^i, being traceless implies 

Pe = ^Pe. (1.30) 
It is possible to show D^n^j^i, = 0, hence 

D.Ei^ = 0, (1.31) 

for this cosmological solution. So fll.30p and fll.Sip give us 

3H{pe + Pe)+Pe = 0, (1.32) 

where H = ^. Eq. (11.321) can be integrated to obtain 

Pj^ = Ca-^, (1.33) 

where C is just an integration constant. With the previous results and after a few 
algebraic manipulations of (11.251) we obtain the modified Friedman's equation 



H =-- + ^ + ^ + —ir^ + 1-34 

3 3 3 36 



with Pe = 

The term (comes from the quadratic term in r^^^) will have important 
implications in the very early universe, when the energy density is high. The pe 
term is called "dark radiation" because of its behaviour, and it is related with 
the mass of a black hole that exists in the bulk. 



1.4 Low energy expansion method 

As we have seen in the previous section, the brane-world Einstein equations are 
not easy to solve and only a few exact solutions are known. This is because 
of the term E^,^. This term E^^, in the projected 4D Einstein equations makes 
the system non-closed. As a consequence if one wants to find the 4D dynamics 
one has to determine E^^, by solving the full 5D non-linear differential equations 
[8]. Motivated by this necessity, Kanno and Soda (KS) developed a perturbative 
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method to solve the 5D equations and to find Efj,,^. I shall now describe their 
method which was applied to RS 2 branes model in [TTl [T2] . Because we will 
consider the matter terms in (11. 3p . the branes will not in general be fiat. Conse- 
quently we cannot put both branes at y = and y = I and use Gaussian normal 
coordinates. The line element that captures the problem's symmetries is 

ds^ = e2^(^'"")rf|/2 ^ g^^i^y^ x^)dx^dxr (1.35) 

And so, the physical distance between the A-brane {y = 0) and the B-brane 
{y = /) is ^ 

d{x) = [ e^^y'^^'^dy, (1.36) 



(/)(y, x^) is the so-called radion. The 5D Einstein's field equations are given by 
+4 [\5>B+ fi: -l^i: 1 1 e-^6{y - /), (1.37) 



1.38) 



= -p-f (-4aA+ e-<^5(y) - f (^-4crB+ e-H{y - I), 

V. - (e-'^ir) = 0, (1.39) 

where denotes the covariant derivative with respect to g^^, and the extrinsic 
curvature is K^^, = —^g^^^y. From the singular behaviour of the equations (11.371 - 
11.391) we find the junction conditions 

e-f [Kt^ - 5^;if]|,=o+ = ^ (^-crAS^.+ , (1-40) 



,2 / B 

2 



e-^ [Ki: - Si:K]\y=i- = ( -aBSi:+ Tj^ | . (1.41) 



Using the definitions 

= S^, + -^g^,Q, Q = -e^"^^ In v^- (1-42) 
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The system of equations that we will solve is: 



VaeJ - ^v^g = 0, 

for the bulk equations (obtained from (11.371 11.381 11.39P ). and 



4 



:i.43) 

:i.44) 
:i.45) 
:i.46) 



E^: - -^si^Q 



E(; - ^si^Q 



=0+ 



G 



-^a6^,+ TI: 



-^B6'i+ Tji 



:i.47) 



:i.48) 



y=l- 



for the junction conditions. 

In order to solve a different iable equation using perturbation theory we need 
to identify a suitable small parameter of our equations. The quantities involved 
in the equations have two very distinct characteristic scales. The characteristic 
length scale of the extra-dimension is / ~ 10~^cm and g^v,y ~ The charac- 
teristic brane's curvature length scale is D and is given by R ~ . Then from 
the 4D Einstein's equations we see pi ~ (■i)G^j:,2 equations (11.41 11.61) imply 
~ (4)J^;2 ; SO ^ ~ ("o)^' Astrophyslcal (or cosmological) events have length 
scales D ^ /, so the natural small parameter is e = ^ ~ {15)^ ■ This means that 
the energy density of the matter is small compared to the brane's tension but it 
does not necessarily imply weak gravity. 

The next steps in the derivation are completely analogous to any perturbative 
scheme. We expand the unknown functions in a series 



ht,u{x^)+ {y,x^)+ g% {y,x^) 



+ . 



:i.49) 
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with the boundary conditions on the A-brane reading 

9>^u{y = 0,x'') = h^A^^), (1.50a) 

(n) 

gj,{y = 0,x>') = 0, withn = 1,2,--- . (1.50b) 

Then we plug them back into the initial system and we will obtain successive 
systems for the different orders in e. The 0*'' order system has the solution 

ds^ = e^^^y'^'^dy^ + a\y, x)h^^{x)dx^dx'' , (1.51) 

with a{y,x) = exp [— j dye'^^'^'^^] . To proceed we will assume 4>{y,x) = (f){x) 
thus a{y,x) = exp [— fe"^] . In the I''* order system, the curvature term that has 

been ignored in the 0*'* order calculation comes into play. After solving the set 

(1) 

of equations, it is possible to show that 9 is written in terms of h^u{x), 4>{x) 
and an integration constant Xtiv{x) (traceless and transverse) that basically is the 
Weyl's tensor contribution. The junction condition for the A-brane gives 

^G'e(/^) + xe = ^^^ (1.52) 

where Xu is called generalized dark radiation tensor and contains the effect of the 
bulk and the second brane. Combining (11.521) and its analogy for the B-brane 
one obtains 

= ^t}+^^t}+1{^\:-6^^^\: 



*2 "^1- 2 " 

where ^ = 1 - fi^ = i _ g-s^^ and tu(^) = f^^. The radion field obeys 

„ = 4r±Z (1.54) 

I 2uj + 3 2uj + 3d-^ ^ ^ 

The important point about equations f ll.53p and f ll.54p is that they do not include 
the Weyl's tensor term although, they include the energy momentum tensor 
of the B-brane, KS called this theory "quasi-scalar-tensor" gravity. 

Equations fll.53p and fll.54p are the only equations we need to solve to deter- 
mine the full spacetime dynamics. Let us see what these equations tell us about 
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the cosmology. 

Because we are interested in isotropic and homogeneous cosmological models 
we will choose h to be the FRW metric. We will assume that the matter on 
the branes has a perfect fluid form and so the energy momentum tensor of the 
A-brane and B-brane will be respectively 

l| = dtag{-pA,PA,PA,PA), (1.55a) 
2I = n^dtagi-pB,PB,PB,PB), (1.55b) 

where the Q'^ factor results from the fact that the B-brane's metric is Vt^hp^p. 
The symmetries imply that \1/ only depends on time. Using the physical distance 
between branes instead of ^ (related through ^ = 1 — e"^'^/') we can write (11.531) 
as 

3 (1 - e^'"') (h' + ^)=^{pa + Pse-''/^) +3y [y - 2h] e~''l\ (1.56a) 



_ (1 _ e-^^/i) (2^ + H' + \\ =^(p^ + Pse-'"/') 

\ Ct (X I i 



, 2d £ AHd , , 



and equation (11.541) as 



,-^ + 3J/d = 4'"^"'^^'^'^/"''^"'^""'. (1.57) 

Combining these three last equations we can find 

''^ + H^ + -,=%{Pa-Wa). (1.58) 

If we assume the equation of state Pa = ujpA, the usual energy conservation laws 
in the A-brane imply pA ~ a"'^^^^'^-*. Now ( 11.58^ can be integrated once to give 



This is the well known low energy limit of the modified Friedman's equation that 
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we have already encountered before in subsection ll.3.1[ Once again C appears 
just as an integration constant. 

As a last remark, I should point out that the effective equations of motion for 
gravity on the brane can be consistently derived from the action 



The solutions at higher orders in the method will be left out of this introductory 
chapter, see [HI |T2] for further details. 

In the next three chapters of this thesis we will apply the gradient expansion 
method that we have just presented to obtain the 4D effective theory of gravity 
for three models with different number of extra-dimensions and different physical 
properties. The models that we will study are: a 5D dilatonic two brane model 
that reduces to the Hofava-Witten theory and to the Randall-Sundrum model 
for a particular choice of parameters, a 6D supergravity model and the lOD type 
IIB supergravity. The next three sections of this introduction shall be devoted 
to the presentation of the essential features of these three models. 

1.5 The Hofava-Witten theory 

In this section we will introduce the general features of the Horava-Witten theory 
[21 [22] . Horava and Witten showed that the eleven-dimensional limit of M-theory 
compactified on 3^/7^2 orbifold can be identified with the strongly coupled lOD 
Eg X Eg heterotic string theory with two orbifold fixed planes (branes). This model 
provides a natural setup in which standard model matter fields are localized on 
the branes. A schematic representation of this model can be found in figure 11.21 
Witten also showed that this theory can be consistently compactified so that the 
universe appears as a 5D spacetime with two boundary branes [23]. The six extra 
dimensions were compactified in a Calabi-Yau manifold and the action of the 5D 
effective theory was derived [3], [21]. These references showed that although the 
dimensional reduction of the graviton and the 4-form flux generate a large number 
of fields, it is consistent to retain only the 5D graviton and a scalar field 0. This 
scalar field is related to the volume of the Calabi-Yau manifold as Vcy = The 



S 




(1.60) 
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6D 6D 




Figure 1.2: The IID Horava-Witten theory compactified on a orbifold with 

two end of the world branes. The six small extra dimensions are compactified on 
a Calabi-Yau. 

resulting dimensionally reduced efi^ective action in the bulk is 

S^uik = 7^2 f d'xV^5 (^'^ - dM(pd'U + Vtuik{<l>)) , (1.61) 

where Kq denotes the 5D gravitational constant and Vbuik is the potential of the 
bulk scalar field. It is given by 

V,uik{<P) = (1.62) 

where a is the tension of the 3-brane. 

The topology of the 5D spacetime is S'^/7j2 x M4 where M4 is a smooth 
manifold. We label the coordinate of the extra dimension by y and the two 
endpoints of this space are located at y = and y = r. We place the positive 
tension brane at y = and the negative tension brane at y = r. The branes' 
actions are given as 

S+ = -^a I d^a;v^e^^^ S_ = I ci^^v^e^^^ (1.63) 

where the subscript plus or minus refers to the tension of the respective brane, 
so that for instance 5+ denotes the action of the positive tension brane. 

The cosmological evolution of this model including matter on the branes and 
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arbitrary potentials on the branes and in the bulk have been studied in [25]. In 

|25j . the authors also obtain the 4D effective equations that govern the low-energy 
dynamics of the brane-worlds. Reference [26] showed that this model possesses 
a unique global solution that represents the collision of the branes. They also 
showed that away from the collision the Riemann curvature is bounded and they 
argued that it offers the prospect of modelling the big bang singularity as a brane 
collision in an higher-dimensional spacetime. This scenario is called Ekpyrotic 
universe and early proposals can be found in [271 128] • Some dynamical solutions 
have been found [29], but they are unstable. In these solutions both branes 
collide with a time-like naked singularity resulting in the annihilation of the entire 
spacetime. In the next chapter, we shall investigate the origin of this instability 
using the gradient expansion method. We shall obtain the 4D effective theory 
and we will show that the previous exact dynamical solutions can be reproduced 
within the 4D theory. 

1.6 The self-tuning 6D brane- world 
1.6.1 The 6D Einstein-Maxwell model 

Co-dimension two brane-worlds in a 6D bulk have recently attracted much atten- 
tion as they possess the remarkable property that the brane's tension does not 
curve the 4D spacetime but its only effect is in the extra-dimension, see [30] for a 
review. This property was called self-tuning and was interpreted as a promising 
solution for the cosmological constant problem. 

In Chapter [3] we will study a 6D supergravity brane- world model and we shall 
obtain the low energy effective theory. Before as an introduction let us see what is 
the cosmological constant (CC) problem and how the 6D Einstein-Maxwell model 
tries to solve it. The CC problem consists in the huge difference between the value 
of the vacuum energy density predicted by particle physics, p^ac ~ (TeV)'^ and 
the observed value of pyac ~ {lO^^eV)"^. There is no convincing explanation of 
why the observed value is so small and why does the CC only started to dominate 
the energy density of the universe recently. 

In the 6D Einstein-Maxwell model [31], the action is 



S 
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(1.64) 
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where Mg is the 6D Planck mass, Ag is the 6D cosmological constant and Fab 
is a Maxwell field that is used to stabilize the size of the extra-dimensions, '-^^i? 
denotes the 6D Ricci scalar. The ansatz for the line element is 

ds'^ = rj^ydx^dx" + 'jijdx'^dx^ , (1.65) 

where 77^1, denotes the Minkowski metric and jij denotes the internal metric. The 
gauge field is chosen to have non-zero components only in the compact directions 
and it takes the form 

Fij = y/^BQeij, (1.66) 

where Bq is a constant, 7 is the determinant of 7jj and ejj is the Levi-Civita 
tensor. A static and stable solution is El], 



-fijdx'dx^ = al {de"^ + sin^Odip'^) , (1.67) 

where the radius of the 2-sphere oq and the magnetic field strength Bq are given 
in terms of the bulk cosmological constant as 



ag = M, Bl = 2K. (1.68) 



In order to have a 4D Minkowski spacetime, Bq has to be tuned to the value of 
the bulk's CC. A different value of Bq would induce a de-Sitter or anti-de-Sitter 
4D geometry. This is the CC problem. 

Now we can introduce branes with tension at the poles of the sphere. The 
branes' action is 

Shranes = ~^ j ^^^V^' (1-69) 

where a is the branes' tension. The inclusion of branes will only modify the 
extra-dimensions geometry to that of a "rugby ball" 

-fijdx'dx^ = al [de'^ + ahiTi^Od^'^) , (1.70) 

where 

a = l--^. (1.71) 

27rM| ^ ' 

The internal space can de depicted as a sphere with a wedge (stretching from the 
north pole to the south pole) taken out and the open sides glued together, see 
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Conical singularity 




Figure 1.3: The rugby ball geometry with two branes of tension a at the conical 
singularities. 

Figure [T73l In other words, the branes introduce a deficit angle 6 = 27r(l — a) in 
the sphere. Now, we can easily see that the 4D geometry is independent of the 
branes' tension a. The branes' vacuum energy does not curve the 4D spacetime 
as it would in 4D general relativity. This property was called self-tuning. It does 
not solve the CC problem because we still need to tune the value of Bq, but this 
idea looks promising as it breaks the standard relation between the 4D vacuum 
energy and the curvature of the 4D spacetime. Also, if the self-tuning mechanism 
worked, it would push the value of the CC to zero but we would still have to find 
a way to generate a small CC to agree with observations. 

There are several objections to the idea of self-tuning [32]. The simplest 
argument is as follows. Suppose that we have tuned the initial value of Bq to Bq = 
2Aq and suppose that the branes' tension were to change from an initial value to 
a different final value. This would imply that the area of the extra-dimensions 
(Ana) would change too, because a is a function of a. On the other hand, the 
magnetic flux in the internal space has to be conserved, this necessarily implies 
that the flnal magnetic fleld strength cannot satisfy -Bo/inai ~ ^Ag. Therefore 
the 4D spacetime will no longer be a static Minkowski spacetime but rather a 
time-dependent spacetime. 

1.6.2 Supersymmetric large extra-dimensions 

In the previous subsection, we have seen that the failure of the self-tuning is 
directly related to the tuning of the magnetic fleld strength with the bulk CC, 
which is necessary to obtain 4D Minkowski spacetime. The so-called supersym- 
metric large extra-dimension (SLED) model was proposed to try to address this 
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problem. For excellent review articles see [301 [M] . 

The action of this model is similar to the action of the Einstein-Maxwell model 
with an additional field d>, the dilaton, it reads 



/r A/f ^ 



Mi X ^ 1 ^ 

" - -e-fF^.F'^" . (1.72) 



Comparing this model with the Einstein- Maxwell model it is easy to see that 
there exists a solution with (p = 00, where 0o is a constant. It is also a solution 
of the Einstein-Maxwell model with the rescalings Ag Aqc'"^, Bq B^e"^, 
provided that (j) is at the minimum of its potential 

V'{(t>o) = -\Ble-^^ + A^e^ = 0. (1.73) 

This is the condition to have a flat geometry on the brane SqC"'^'' = 2A6e'^. 
In other words, at the minimum, the vacuum energy is zero and therefore the 
geometry is Minkowski, this was called self-tuning. 

As before, there are several objections to the idea of self-tuning, for example, 
Refs. [32l [35] derived the 4D effective theory using the metric ansatz 



g^^{x)dx''dx'' + Mfe-^^^^^ {dr^ + sm^rdO^) , (1.74) 



and (j) = (f){x). The 4D effective potential is found to be 

Vi^|J, (j)) = M^^e'^'Uiai), U{a,) = ^e'^., _ sM^e"'^^ + 2A6, (1.75) 

where ai = 2^/^ + 0, (72 = 2^/^ — and a is related to the branes' tension as 
a = 1 — 2"m^ • They argued that if we start with a configuration in which 
V{il), (p) = 0, and if we change the branes' tension to a different final value then 
V{ip, (f)) will be different from zero and (T2 will have a runaway potential. The 4D 
spacetime becomes non-static. 

This argument against self-tuning has raised criticisms [3S]- For example, it 
is argued that the metric ansatz fll.74p is restrictive. There is a known class of 
static solutions with warping in the bulk that cannot be described by this ansatz. 
These solutions might play a role if self-tuning is to work. Chapter [3] is dedicated 
to the study of some aspects of this model and we shall derive the low energy 
effective theory including warping. 
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Figure 1.4: A warped throat compactification. 



1.7 Warped compactifications in Type IIB su- 
per gravity 

In chapter m we will present a method to obtain the 4D effective theory for warped 
compactifications including fluxes and branes in the lOD type IIB supergravity. 
In this section we shall introduce the general idea of warped compactification in 
type IIB supergravity (sugra). For further details about type IIB sugra we refer 
the reader to chapter H] or to the extensive literature on the subject [3S1 EH EB] 
and references therein. 

The desire of embedding the Randall and Sundrum (RS) solution to the hi- 
erarchy problem in a string theory compactification was fulfilled with the work 
of Giddings, Kachru and Polchinski (GKP) [13]. They realized that a warped 
string compactification can generate a large hierarchy between the electroweak 
scale and the Planck scale in a similar way to the RS model To generate 
the warping of the extra-dimensions they needed to use not only positive tension 
objects like D3 branes but also negative tension objects. Figure [L4l is a schematic 
representation of their model. 

They start with a lOD spacetime, with six of the dimensions compactified in 
a Calabi-Yau manifold. In some region of the internal space, a throat develops. 
Fluxes and branes curve this space, giving origin to a warp factor h{y). Their 
solution for the metric is 



(1.76) 
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where 7mn and y"^ denote a Calabi-Yau metric and the coordinates of the extra- 
dimensions respectively, denotes the four- dimensional coordinates. The throat 
region of the spacetime is given by the solution found by Klebanov and Strassler 
[39] and it is well approximated by a 5D anti-de-Sitter spacetime with the re- 
maining five extra-dimensions compactified in a small 5-sphere. In this region 
the RS mechanism of solving the hierarchy problem works as before. 

These warped compactifications also provide models of inflation, where the 
infiaton is identified with the position of a brane in the extra-dimension r, where 
r is a radial direction and it corresponds to one of the y"^. Brane-inflation models 
will be discussed in the chapter |5l 



Chapter 2 



Moduli instability in warped 
compactification 

2.1 Introduction 

The idea of using the degrees of freedom of extra spatial dimensions to unify grav- 
ity and the other interactions has attracted much interest. Especially, M-Theory 
has the potential to achieve this ultimate goal. Hofava and Witten showed that 
M-Theory compactified as jZ^ x reduces to a known string theory [E^ x 
Heterotic String) [2]. From a low energy effective theory point of view, this model 
has a 5D spacetime (the bulk) with two 4D boundary hypersurfaces (the branes) 
[3]. The remaining six spatial dimensions are assumed to be compactified. All 
the Standard Model particles are confined to the branes while gravity can prop- 
agate in the bulk. As a consequence of the compactification of the six spatial 
dimensions, a 5D effective scalar field appears in the bulk, which describes the 
volume of the compactified 6D space. 

A novel aspect of this model is that the fifth dimension is not homogeneous. 
The line element for a static solution is given in the form ds\ = h°'{r)dr'^ + 
h^{r)ds1{x), where r is the fifth coordinate, a, (3 are constants and h is called 
the warp factor. In the previous chapter, we showed how this warped geometry 
was used by Randall and Sundrum [5] to address the hierarchy problem. They 
have also showed that it is possible to localize gravity around the brane with the 
warped geometry, providing in this way an alternative to compactification [1]. 

A new class of dynamical solutions that describes an instability of the warped 
geometry has been found. Chen et al. noticed that it is possible to obtain a 
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dynamical solution by replacing the constant modulus in the warp factor hhy a 
linear function of the 4D coordinates |29]. This solution describes an instability 
of the model as the brane will hit the singularity in the bulk, where h = 0. This 
kind of solution exists also in the lOD type IIB supergravity [40] . 

In this chapter, we study the moduli instability in a dilatonic two brane model, 
where the potentials for the scalar field on the brane and in the bulk obey the 
Bogomol'nyi-Prasad-Sommerfield (BPS) condition [H]. For particular values of 
parameters we retrieve either the Hofava-Witten theory or the Randall- Sundrum 
model. We identify the origin of the moduli instability using a 4D effective 
theory derived in Ref. [25] (see Refs. [12] and [13] for different approaches). 
This effective theory is derived by solving the 5D equations of motion using the 
gradient expansion method [121 SH SS], where we assume that the velocities of 
the branes are small compared with the curvature scale of the bulk determined 
by the warp factor. 

Despite the fact that the 4D effective theory is based on a slow-motion approx- 
imation, we will show that the 5D exact solutions can be reproduced in the 4D 
effective theory. In order to understand the relation between 4D solutions in an 
effective theory and full 5D solutions, we revisit the gradient expansion method 
by employing a new metric ansatz. Using this metric ansatz, we can clearly see 
why the moduli instability solution in the 4D effective theory can be lifted to an 
exact 5D solution. 

We also comment on a claim that the 4D effective theory allows a much wider 
class of solutions than the 5D theory [46]. We disagree with that conclusion and 
we show that it is based on the restricted form of the 5D metric ansatz used in 
Ref. [inj. Using our metric ansatz, we provide a way to lift solutions in the 4D 
effective theory to 5D solutions perturbatively in terms of small velocities of the 
branes. 

The structure of this chapter is as follows. In section 12. 2[ the model under 
consideration is described in detail. In section 12. 3[ we identify the solution in 
the 4D effective theory that describes the moduli instability. In section 12.41 we 
revisit the gradient expansion method to derive the 4D effective theory. We 
propose a new metric ansatz which is useful to relate 4D solutions in the effective 
theory to 5D solutions. Using this formalism we explain why the 4D solution for 
the moduli instability can be lifted to an exact 5D solution. In section 12. 5[ we 
comment on the arguments against the 4D effective theory. Section [216] is devoted 
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to conclusions. 



2.2 The model 

Our model consists of a 5D spacetime (bulk) filled with a scalar field. The fifth 
dimension is a compact space with a Z2 symmetry (i.e. identification r —r) 
and will be parameterized by the coordinate r, < r < L. The bulk action takes 
the form 

Sbuik = 7^ c^^xv^ [^R - Vm^V^'^ + Vbuiki^)] , (2.1) 

where kq denotes the 5D gravitational constant. Throughout this work Latin 
indices can have values in {r, t, x, 2}, while Greek indices do not include the 
extra dimension coordinate r. The scalar field potential will be 

V,^M = e-^v^'V^ (2.2) 

where a and h are the remaining parameters of this model. We can retrieve 
either the Randall- Sundrum model (6 = 0) or the Hofava-Witten model (6 = 1) 
according to the value of h. The orbifold /Z2 has two fixed points, at r = 
and r = L, and we can put branes there. The branes' action is 



'branes 



+ (2.3) 



where S+ and are the positive and negative tension brane action, respectively. 
They are given by 

= [ d^xy^.e-^'^a. (2.4) 

The action of our model is 

Stotal — Shulk ~\~ Sbranes- (2-5) 



A schematic picture of our model is shown in Fig. 12.11 
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^ M'^x S/^ 



(|) (x,r) 



^2 ^2 

Figure 2.1: The dilatonic brane- world universe 

2.3 Exact solutions for moduli instability in 4D 
effective theory 

2.3.1 The exact 5D solution 



Exact static solutions for the equations of motion that result from (12.51) were 
obtained in [H]. An interesting dynamical solution was found by Chen et al. 
|29j . They noticed that if we add a linear function of time to the warp factor of 
the solution of [H], it would still be a solution of the equations of motion. Their 
exact 5D solution reads 



dSl 



hr 



dr"^ + hr 



jjiiy doc 



3V2b 
362 + 1 



In hr 



V2 



362 + 1 a ' 



(2.6) 
(2.7) 



where h is an arbitrary constant. 

From the above equations we can read the scale factor of the positive tension 
brane (r = 0) as 



ai(r) = (/ir)36 



+1 



and the scalar field as 



if{r = 0,r) 



3^26 
362 + 1 



ln(/ir). 



(2. 



(2.9) 
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Let us choose h < and r < 0. In this case the proper distance between 
branes, the so-called radion, is decreasing according to 



n 



/ir — y j dr 



(2.10) 



oo —J- „ uu „ 

When r ^ -oo, 7^ ^ L (/ir)^^ and for r = 7^ = l{3b^+l) {L/l)^^^ 

Before the two branes collide, a curvature singularity will appear in the negative 
tension brane (r = L) at r = Ljhl. This singularity will move towards the 
positive tension brane and will reach it at r = 0. This event represents the total 
annihilation of the spacetime. 

It is useful to note that if we drop the modulus sign in (12. 6p . the bulk spacetime 
is a static black brane or black hole solution depending on 6, see [29] for more 
details. If 6 = 1, there is a timelike curvature singularity at hr = \r\/l. From 
the static bulk point of view, the two branes are moving in this static bulk and 
the negative tension brane first hits the singularity. Even if the bulk spacetime 
is static, the existence of the branes, which gives the modulus sign for r, makes 
the spacetime truly time dependent. 

Moduli instability is a serious problem for these types of model. In this work, 
we will try to understand it from a 4D effective theory viewpoint. 



2.3.2 The exact solutions in 4D effective theory 

In [25], Kobayashi and Koyama applied the gradient expansion method to solve 
perturbatively the 5D equations of motion resulting from (12. 5p . Their O*'^ order 
solution reads 

dSl = d{tfe^^'"^^'Ur^ + Tir,t)i^h^,{t)dx^dx'', (2.11) 

^{r,t) = ^lnj^(r,t) + 0(t), J^{r,t) = 1 - V2b'^d{t)a\r\. (2.12) 
v26 

The main result of the paper is a set of equations involving only four dimensional 
quantities that describe the dynamics of the unknown functions d(t), h^y{t) and 
0(t), appearing in Eqs. (12. lip and (12.121) . It can be shown that their dynamical 
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equations can be consistently deduced from the following action 



S. 



ipR{h) 



2(1 + 362) l-ip 



(2.13) 

where here means the covariant derivative with respect to h^j^^, and ip = 

3b^+l 

1 — (1 — d) . This is the action we obtain if we substitute Eqs. (12.111) and 
(12.121) into (12. 5p and integrate over the extra dimension [12]. 

Performing the conformal transformation /i^,^ = e~^^'^ f ^j^^ / ijj and defining new 
scalar fields as 

,2 



ijj = 1 — tanh 

e 



1 + 362 



6 



+ Aln^, A 



we arrive at the action in the Einstein frame 



Se = ^ I d'xy^ [Rif) - ei^ei" - vi>|„vi>i"] . 



(2.14) 
(2.15) 

(2.16) 



It is clear that the moduli fields have no potential and this is the origin of the 
instability. We are interested in cosmological solutions so we choose f^j^p to be a 
fiat FRW metric. The Einstein equations resulting from the action (12.161) are 



a 6 



a \ 1 



(2.17) 



and the equations of motion for the scalar fields are 



+ 2-^' = 0, e" + 2-e' = o, 



(2.18) 



where a(r) is the FRW scale factor and the prime denotes derivative with respect 
to the conformal time r. Equations (12.17112.181) can be easily integrated to find 



6 



362 + 1 



alna2 + 7, e(r) = v^TTsP/Jlna^ + 5, a2(r)=^r + C, 

(2.19) 

«2 + ii±^/32 = 1 + 3^2^ (2.20) 
6 



where a, /3, 7, 5, C, and ^ are integration constants. The solutions for the fields 
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in the original frame can be obtained as 



dir) 



1 - 



(/ir 



\2a 



In 



(/ir + z/)2" + 1 



(2.21) 



(2.22) 



where A (6) is defined in Eq. (12.151) and integration constants fi, u, g are redefined 
from C,, C; 7 and 6. The radion is calculated as 



-d{T)e 



and the square of the scale factor on the positive tension brane is 



(2.23) 



-a\r). 



(2.24) 



We have found the remarkable fact that for a particular choice of the integration 
constants we can reproduce the exact solution found by Chen et al. described in 
the previous section. If we choose integration constants obeying the relations 



1 



a 



2hl 



V 



(2.25) 



we get the same 4D quantities (scale factor and scalar field on the positive tension 
brane) and radion as the Chen et al. solution, Eqs. (12.81 12.91 l2.1Up . 
For this solution, the scale factor in the Einstein frame is given by 



2L 

T 



hr 



21 



(2.26) 



If we take h < 0, t < 0, this corresponds to a collapsing universe, due to the 
kinetic energy of the scalar fields. At r = L/2hl < 0, the universe in the Einstein 
frame reaches the Big-bang singularity. However, we should be careful to interpret 
this singularity. In fact, in the original 5D theory, r = L/2hl does not correspond 
to any kind of singularity. In 5D theory, the negative tension brane hits the 
singularity at r = L/hl and the positive tension brane hits the singularity at 
r = 0. In fact at r = L/2/i/, = and so the conformal transformation becomes 
singular and the Einstein frame metric loses physical meaning. 
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2.4 Gradient expansion method with a new met- 
ric ansatz 

In the previous section, we found that the exact solution derived in [29] can be 
reproduced from the 0*'* order of the perturbative method. In this section, we will 
revisit the gradient expansion method which is used to derive the effective theory 
to clarify the reason why the exact solution derived in [29] can be reproduced 
within the 4D effective theory. 



2.4.1 5D Equations 

In this subsection we study the 5D equations of motion. From the action fl2.5|] 
we derive the Einstein's equations 



'Gab = VAf'^Bf + 7^9 AB [-V^c^VcV^ + H^ifclv^)] 



9,uK^''Be-^^''^5{r) 



9iiv 



5^5^e-^^^5(r - L) 



(2.27) 



Assuming that the 5D line element has the form 



dSl = gAB{X)dX'^dX^ = grrix,r)dr' + g^^{x,r)dx''dx'' 



(2.28) 



{x dependence means dependence of all the other coordinates {x,y,z,t}, except 
the extra dimension r) we can extract the junction conditions for the metric 
tensor from Einstein's equation, 



.=0+ 

r=L- ^ 



r=0+ • 
r=L- 



(2.29) 



where the extrinsic curvature K^^, is defined by K^^i, = —\drgfj,v. The scalar field 
equation of motion is 



'9"e 



r)-y/g" e ^^^6{r - L) , 
(2.30) 
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and so the junction conditions for the scalar field are 



=0+ 



^bae 



-V2hip 



r=0+ 
r=L~ 



(2.31) 



In order to proceed we shall assume that the 5D line element has further 
symmetries, described by 

dSl = e^^^^^''''^H^^{x,r)dr^ + {x,r)g^^dx^'dx'' , H{x,r) = C{x) - -. 

(2.32) 

In this section, we will assume that the position of the second brane is r = L = /. 
For the scalar field we will assume the form 



(p{x,r) 



3V2b 



In if + ip{x, r) 



(2.33) 



This metric ansatz is inspired by the time dependent solution fl2.6|2.7p of Chen 
et al. [29]. Their solution was found by replacing the modulus parameter in the 
static solution by a linear function of time. In the same manner, we introduce an 
X dependence in H through the modulus parameter C{x) in a covariant way . We 
also introduce the function (p for the scalar field moduli. In order to satisfy the 
junction conditions f l2.3ip we must have the exponential factor in the grr metric 
component. The tensor g^y is left completely general. 

After some mathematical manipulations of the Einstein equations we obtain 



„ - ~9b^+i / 36^ — 5 ~ 



-^f2he-^'•^^''H' 



1 



l/l 



3fe2 + l 



362 + T'' 



-3i)^ + l ~ 



-H 36^ + 1 



V2b 



1^ 362 + l"^l" 



V2b^\^^ + (26^ + 1) ^\^^^^^ 



(2.34) 
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1 36^ — 3 o ATt - -9b'^ + l ~ o AT, - -66^+2 / ~ ~ a ~ 

1^362 + 1^1- + 362 + 1^^ ^1" 
-V2bH~^^ + v^6(^l"^|„) = 0, (2.35) 



362 



362 + 1 



362 



362 + 1 



362 



(2.36) 



where K^^, is defined like A'^,^ = —\drg^y, K is the trace of K^u, \a denotes 
covariant derivative with respect to g^^ and R^j,{g) is the Ricci tensor of g^v 
Equation ( 12.30p transforms into 



-36^+5 



\Plb {^^rf - '^,rK + - 



~ 1 3v^6 



/ 362 



1962 



/ 362 + 1 



(2.37) 



>^ 362 + 1 I- ' —^^1"^ ' 352 + 1^ 

With this particular ansatz the junction conditions 02.29112. 311) are significantly 
simplified and read 



r=0+ 

r=r 



=0+ 



= 0, 



0. 



(2.38) 



(2.39) 



It is impossible to solve these equations in general, so in the next section, 
we will solve them using the gradient expansion method up to first order in the 
perturbations. 



CHAPTER 2. MODULI INSTABILITY IN WARPED COMPACT. 



45 



2.4.2 The gradient expansion 
The approximation 

In order to use perturbation theory to solve a system of differential equations 
we need to identify the characteristic scale of the different terms involved in the 
equations and then see if there is a small parameter. 

The derivatives along the extra dimension of the conformal metric Qfj^i, as well 
as the derivative of ip are of order l/l. Typically, we take I to be of order of a few 
tenths of a millimeter. If the characteristic brane's curvature length scale is D 
then d'^g^y ~ 1/ D'^g^^. We assume that variations along the branes' coordinates 
are small in comparison with l/l. This implies that the radion changes slowly or 
that the speed of the branes is small. More precisely, our small parameters will 
be 

fdl{...)^l and (/9,...)'<1, (2.40) 

where . . . represents the conformal metric functions, C or 

As in the usual perturbation method, we expand the unknown functions in a 
series 

(0) (1) 

9i,u{r, x)=g ^y (r, x)+ g (r, x) + ■ ■ ■ , (2.41) 

ip{r, x) =ip (r, x)+ if (r, x) + ■ ■ ■ . (2.42) 
We impose the boundary conditions at the position of the positive tension brane 

(n) 

g ^^{r = 0, x) = 0, for all n > 0, (2.43) 

(n) 

ip (r = 0, x) = 0, for all n > 0. (2.44) 
Other quantities are naturally expanded as 

(0) (1) 

Kf.u =K^, + k^,, + • ■ • . (2.45) 

In practical terms, at zeroth order in the gradient expansion method we ignore 
all derivatives along the branes' coordinates. These terms will only enter the first 
order equations. At zeroth order, the 5D partial differential equations of motion 
reduce to simpler ordinary differential equations on the extra dimension coordi- 
nate. The gradient expansion method has also been used in the 4D cosmological 
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context imilHlllglEolEl]. 

Qth Order (Background geometry) 

The 0*^' order system can be easily integrated with respect to the extra dimension 
coordinate r to get the particular solution 



jo) (0) 
(0) (0) 

(p (r, x) =(p (x). 
This solution clearly satisfies the O^'^order junction conditions. 



(2.46) 
(2.47) 



^st Order 

At I''* order the evolution equations are 
1 „-£!±i f3b'-5 ill 1 



r 



(1)' 



5^ K \ + e 



2 (1) 



, -r — = — C\"'5^ 1 - H~^^ 



(0) (0) (0) 

V2h (^jC + (26^ + 1) <^l/^<^|. 



H 362+1 



362 + 1 

V2h 1 
/ 362 + 1 



, (0) (0) \ (0) 

(36^ + 1) C\,+ + <^l" C\J^, 



_ (0) -962+1 (1) 2 /(O) 



(2.48) 



1362 



/ 362 + 1 



-36^-3 
— H 36^+1 



Q ^ (0) 



2 (1) 



■2^26 



(0) 



2 (1) 



•^H^^ K, 



362 



362 + 1 |a 352 + 1 



6v^63 W 

V9' C|a 



/ 362 + 1 



(0) (0) 



0, (2.49) 
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(0) 



(0) 



I (1) (1) 

-H~^ 3^26 k + (96^ - 5) (p^r 



, 3v/2& /- (0) 962 + 1 (0) 



362 + 1 



0. 



(2.50) 



The junction conditions at this order are 



■ (1) ■ 



r=0+ 
r=l- 



(2.51) 



r=0+ 
r=Z- 



0. 



(2.52) 



In the preceding equations all the indices are raised with the zeroth order metric. 
Combining the trace of equation fl2.48p with equation (12.491) we obtain 



1 ^ (0) (1) 



362 



R 



(0) 



(0) (0)' 



(0) 



^/26 r- (1) 



(2.53) 



Imposing the junction conditions fl2.51|2.52p we get 



'(o)\ (0) W 



(2.54) 



and the equation of motion for the 4D effective scalar field 



(0) 



c\2 + V2bq^ (^1-= 0. 



(2.55) 



Equation (12.531) now reads 



K= -V26 (p,r 
(1) 



Using the decomposition of K^!: in 



(2.56) 



(1) 



(2.57) 
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equation fl2.48p can be easily integrated to find 



r- (°) 362-5 



/ (0) 

(0) (0) 



/ /(0)\ (0) (0)\ / 

-iR^^i^~gj-V2b^\^^- {2b' + 1) ^i.^l'^ j (^Cr - - j 



+ X^ui^), (2.58) 



- traceless 



where Xui^) an integration constant and the subscript [ ]traceiess nieans the 
traceless part of the quantity between square brackets. In terms of S^, the junc- 
tion conditions fl2.5ip are 



r=0+ 



0. 



(2.59) 



From the previous junction conditions (I2.59P we can obtain the 4D effective 
equations of motion 



'(0) 

~9 



C 



1 / (0) (0)- 



/ 1 \ (0) (0) U2 (0) (0) 

+V2b UjC --6i^ + {2b' + 1) ^1.(^1^ --6ii (2.60) 



and 



(2.61) 



Combining equation fl2.56p with the scalar field equation fl2.5UI) and integrat- 
ing it once with respect to the extra dimension, we get (after using the previous 
equations of motion to simplify the result) 



-2V2b 



(0) -3b2 + 5 (1) 



^ (p^,= - Cr 



r2\ I W (0) (0)\ (0) 

' \ i ~\Oi ' I.I I. 



(2.62) 

where H(x) is just an integration constant. 
The junction conditions (12.521) give 



S(x) = 



(2.63) 
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and the equation of motion for the second 4D effective scalar field 



(0) (0) (0) 

+V2h^\^^\^-- 



'1 (0) 



(2.64) 



2.4.3 The 4D Effective theory 

The 4D effective equations of motion are summarized as 



'(0) 

K[ ~9 



^-2 



1 



(0) 



/ (0) 1 (0) \ (0) (0) U2 (0) (0) 

V2h --5^, + (26^ + 1) (2.65) 



(0) 



(0) 

^1 



C\l + V2hC\^ (^l"= 0, 

(0) 



^7 (0) ^'^^ / 1 



-1 (0) 

^1° q^, 



(2.66) 
(2.67) 



and they can be deduced from the following action 



eff 



(0) / 1 

~9 iC-^\e 



r- 



(0) (0) (0)' 
i?(^)-(^|,<^l- 



(0) 



f2.68) 



where \a denotes covariant derivative with respect to g^j,^. We should note that 
this effective action can be derived by substituting in (12. 5p the 5D solutions up 
to the first order and integrating it over the fifth dimension |43j . 

As a consistency check, we see that if we perform the conformal transformation 



(0) 



V(a;) = C3PTT(x) g^^ (x), 



(2.69) 



the previous action reduces to (12.131) if the effective scalar fields of the two theories 
are related through 



^ = 2C"2 - - 



3V2b (0) 

lnC+ ip [x). 



(2.70) 



(2.71) 



362 + 1 

The check consists in seeing that these relations are exactly the ones required so 
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that the two observables (the proper distance between branes and the scalar field 
on the positive tension brane) agree in both approaches. 

2.4.4 The 5D exact solution 

It is straightforward to find a cosmological solution of this 4D effective theory. 
For example we can easily find the following particular solution 

(0) (0) 

gfiu (x) = T]^^, C{x) = ht, (p [x] = 0, (2.72) 

where h is an integration constant. 

Now we are ready to address the question why the above solution can be lifted 
to an exact 5D solution. Let us start by calculating the next order correction ^^^g^i^, 
and ^^^Lp. Eq. (12.621) and the boundary conditions (I2.44p give ^^^(p = 0, if we take 
as O*'* order solution Eqs. (12.721) . We can construct ^^^K^,y from Eqs. (12.561) and 
(12.581) . For the O*'^ order solution (I2.72p this gives ^^''K^u = 0. After imposing the 
boundary conditions (I2.43p . we obtain that the next order correction vanishes, 
^^^(jfMuir, x) = 0. For solution (12.720 it turns out that all the corrections vanish 
and the O*'^ order solution is an exact solution of the non-perturbed 5D Eqs. 
fl2.34tl2.39l) . 

For other solutions of the 4D effective theory, higher order corrections will 
not vanish and therefore they should be taken into account in the reconstruction 
of the 5D metric. Using the gradient expansion method, we can reconstruct the 
5D solution perturbatively. We should emphasize that the choice of the O*'^ order 
metric is quite important in order to reconstruct 5D solutions efficiently. Our 
metric ansatz has the advantage that it is possible to recover the exact solution 
of Chen et al. (12.720 at 0*^ order. Indeed, if we had started with an ansatz like 
Eqs. (12. 11112. 12p we would need an infinite number of higher order terms to obtain 
the exact 5D solution. 

2.5 Validity of 4D effective theory 

In this section we will make comments on a work by Kodama and Uzawa j46j . 
Let us start by briefiy describing their arguments. After deriving the 4D effective 
theory for warped compactification of the 5D Hofava-Witten model (they also 
extend their analysis to lOD IIB supergravity and obtain the same conclusions). 



CHAPTER 2. MODULI INSTABILITY IN WARPED COMPACT. 



51 



the authors show that the 4D effective theory allows a wider class of solutions 
than the fundamental higher dimensional theory. Therefore we should be careful 
in using this effective theory approach, because we may find 4D solutions that do 
not satisfy the equations of motion once lifted back to 5D. 
The authors assume a metric ansatz of the form 

dS'^ = h{x,r)dr^ + {x,r)g^^{x)dx>'dx'' , (2.73) 

where the warp factor has the form h{x,r) = C{x) — r/L. This corresponds to 
taking K^^ = 0, = and 6=1 (for the Hofava-Witten case) in our work. Then 
Eq. (12341) reduces to 

-i/(r,x)"t (^C{x)\''^ + ^C{xpii^ + H{r,x)-'^RXg{x)) = ^ (2.74) 

In order to satisfy this equation for all values of r, we should have 

R^,{g) = 0, q,,, = 0. (2.75) 
They obtain the 4D effective action, by integrating the fifth dimension, as 

Seff oc J d'x^g [c{x) - R{g), (2.76) 

which agrees with our effective action 02.681) . As we have shown, this theory ad- 
mits solutions with R^y{g) ^ (see Eq. fl2.65p ). which do not obey the constraint 
(I2.75P obtained from the 5D equations of motion. 



However, it is clear from our analysis that their metric ansatz is too restrictive. 
If we consider a more general metric as our metric ansatz, we see that the 5D 
Einstein equations contain more terms given by K^^. With the inclusion of these 
new terms, the 5D equations do not necessarily imply (12.751) . Of course, the non- 
vanishing K^y changes the metric fl2.73p and one could argue that the resultant 
4D effective action would be also changed. However, it is shown that even if 
we include the first order corrections ^^^gp.v{x) to the metric, the resultant 4D 
effective action derived by integrating out the fifth dimension does not change 
|43j . Therefore, for 4D solutions that do not satisfy (12.751) . we should include the 
corrections to the metric (12.731) . We have provided this correction perturbatively. 
Using Eqs. ( 12.531) and (12.581) . we can reconstruct the correction to the metric. 



CHAPTER 2. MODULI INSTABILITY IN WARPED COMPACT. 



52 



^^^g^uix), which is necessary to satisfy the 5D equation of motion. 

We should emphasize that the vahdity of the 4D effective theory is based on 
the conditions (12.401) . If the 4D effective theory admits a solution that violates 
the conditions ( ]2.40p . then there is no guarantee that the 4D solution can be 
lifted up to the 5D solution consistently. We should check the validity of the 
4D effective theory by calculating the higher order corrections to ensure that the 
higher order corrections can be neglected consistently. 

2.6 Conclusion 

In this chapter, we have studied the moduli instability in a two brane model 
with a bulk scalar field found by Chen et al.. This model can be viewed as 
a generalization of the Hofava-Witten theory and the Randall-Sundrum model. 
The scalar field potentials in the bulk and on the branes are tuned in order to 
satisfy the BPS condition. 

We used a low energy effective theory, which is derived by assuming that 
variations along the brane coordinates of the metric are small compared with 
variations along the dimension perpendicular to the brane. The effective theory 
is a bi-scalar tensor theory where one of the scalar fields arises from the bulk 
scalar field (dilaton) and the other arises from the degree of freedom of the dis- 
tance between branes (radion). In the Einstein frame, the theory consists of two 
massless scalar fields, and the lack of potentials for these moduli fields was shown 
to be responsible for the instability. 

We found that the exact solution derived in [29] can be reproduced from the 
Qth Qj^(jgj^ Qf ^Y^Q perturbative method, despite the fact that slow-motion approxi- 
mations are used. We revisited the gradient expansion method which is used to 
derive the effective theory, in order to understand why the exact solution derived 
in [29] can be reproduced within the 4D effective theory. We proposed a new 
metric ansatz which is useful to see the relation between the solutions in the ef- 
fective theory and the full solutions for 5D equation of motion. Using this metric 
ansatz, it is transparent why the moduli instability solution can be lifted to a 
full 5D solution. We have also shown that not all solutions in the 4D effective 
theory can be lifted to exact 5D solutions. For these solutions, the solutions in 
the effective theory receive higher order corrections in velocities of the branes and 
we need to find 5D solutions perturbatively. 
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Finally, we comment on the arguments against the 4D effective theory. Ref. 
[16] claims that the 4D effective theory allows a much wider class of solutions 
than the 5D theory. We argued that this conclusion comes from a too restricted 
metric ansatz used in Ref. [46]. Using a more general metric ansatz, we provided 
a way to reconstruct the full 5D solutions from the solutions in the 4D effective 
theory. 

The gradient expansion method can be applied to other warped compactifi- 
cations such as the ones in supergravity models. In fact, there have been debates 
on the validity of the metric ansatz commonly used to derive the 4D effective 
theory in lOD type IIB sugra. Our lOD generalization of the O*'' order metric 
ansatz agrees with that proposed in Ref. |52j and the method presented in this 
chapter will provide a consistent way to reduce the lOD theory to the 4D effective 
theory based on this metric ansatz. This will be the topic of chapter HI Before 
that, in the next chapter, we will obtain the low energy effective theory in a 6D 
supergravity model using the gradient expansion method. 



Chapter 3 



Low energy effective theory on a 
regularized brane in 6D 
supergravity 

3.1 Introduction 

Recently, much attention has been paid to six- dimensional supergravity [531 EH 
[55l [56| [571 [58] . The most intriguing property of six- dimensional supergravity is 
that the four- dimensional spacetime is always Minkowski even in the presence of 
branes with tension. A 3-brane with tension induces only a deficit angle in the 
six-dimensional spacetime and the tension does not curve the four-dimensional 
spacetime within the brane. This feature is called self-tuning and it may solve 
the cosmological constant problem [30], [311 EHl [60] [6l]. This is the basis of the 
supersymmetric large extra-dimension (SLED) proposal [62] . 

There have been several objections to the idea of self-tuning [32 ] [35l [63l [M] [65 ] 
[66], [67]. The self-tuning relies on the classical scaling property of the model. The 
six-dimensional equations of motion are invariant under the constant rescaling 
Qmn ^"^Qmn and e'^~'^, where qmn denotes the six-dimensional metric 

and (p is the dilaton field. Then there is a modulus associated with this scaling 
property. Ref. [32] derived an effective potential for this modulus. This modulus 
is shown to have an exponential potential. Then there must be a fine-tuning of 
parameters to ensure that the potential vanishes in order to have a static solution. 
This is the reason why the static solution always has vanishing cosmological 
constant. However, if this fine-tuning is broken, the modulus acquires a runaway 
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potential and the four-dimensional spacetime becomes non-static. Non-static 
solutions in six-dimensional supergravity have been derived and they are supposed 
to correspond to the response of the bulk geometry to a change of tension of branes 

[MiEniiToiin]. 

However, it is difficult to deal with an arbitrary change of tension with a 
brane described by a pure conical singularity. This is because if we put matter on 
the brane other than a cosmological constant, the metric diverges at the position 
of the brane. Recently, it was suggested that we can regularize the brane by 
resolving it by a codimension one cylindrical 4-brane [721 1131 IZH ESj • These types 
of models may be regarded as a variation of Kaluza-Klein/hybrid brane world 
[76l \77\ [78l [79| [80] . Once the brane becomes a codimension one object, it is 
possible to put arbitrary matter on the brane without having the divergence of 
the metric. Then it becomes possible to study the effect of the change of tension 
on the four- dimensional geometry on the brane. 

There is another interesting issue of whether it is possible to recover conven- 
tional cosmology at low energies in six-dimensional models. Recent works have 
shown that it is impossible to recover sensible cosmology if one derives cosmolog- 
ical solutions by considering a motion of branes in a given static bulk spacetime 
[81], [82]. It was concluded that the time-dependence of the bulk spacetime should 
be taken into account. 

In this chapter, we derive a four- dimensional effective theory for the modulus 
in six-dimensional supergravity with resolved 4-branes by extending the analysis 
of Ref. [83] which studied the low energy effective theory in the Einstein-Maxwell 
theory [HU [85]. Arbitrary matter and potentials for the dilaton on 4-branes 
are allowed to exist. We use the gradient expansion technique to solve the six- 
dimensional geometry assuming that the deviation from the static solution is 
small [m [15]. The gradient expansion method has been applied to various types 
of brane- worlds [861 [HU [Ml [891 [901 [SI]- Using this method, it is possible to 
solve the non-trivial dependence of the bulk geometry on the four-dimensional 
coordinates. By solving the effective four-dimensional equations, we can derive 
the time-dependent solutions and compare them with the exact six-dimensional 
time dependent solutions found in the literature [691 HOI [ZI] • It is also possible to 
study whether we can reproduce sensible cosmology at low energies or not. We 
also study the possibility to stabilize the modulus using the potentials for the 
dilaton on the branes along the line of Ref. [92] . 
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We should mention that we concentrate our attention on a classical dynamics 
in this paper and do not address the quantum corrections. In fact the important 
feature of SLED is that it could also provide a way to address the stability to 
quantum corrections to the cosmological constant [62]. This is because in the 
6D model, the Kaluza-Klein (KK) mass scale is of the order of a fim~^ and so is 
precisely at the energy scale of the cosmological constant. The combination of the 
bulk supersymmetry and scaling solution can maintain the quantum corrections 
to be of order fim~'^ and hence the order of the cosmological constant. 

The chapter is organized as follows. In section 13. 2[ basic equations are sum- 
marized. In section 13.31 we solve the six-dimensional equations of motion using 
the gradient expansion method. In section 13.41 the effective theory on the regu- 
larized branes is derived by imposing junction conditions. Then we derive time 
dependent cosmological solutions in the effective theory and compare them with 
the exact six-dimensional solutions. The possible way to stabilize the modulus is 
discussed. Section 13.51 is devoted to conclusions. 



3.2 Basic equations 



The relevant part of the supergravity action we consider is 



S 



2 2 ^ ^' 



-F^e 



2Zf* 



(3.1) 



is the dilaton, M is the fundamental scale of gravity, := FmnF^'^^ 



^'^^dM4>dN(p and Fmn = QmAn — OnAm is the field strength of the 



where 
(90)2 := g 

gauge field Am- For the moment we are interested in solving the 6D bulk equa- 
tions of motion. In Sec. 13.41 we will add two 4-branes (at positions y = y±) and 
Li denotes the different bulk curvature scales on either sides of the branes, see 
Fig. 13. 1[ We start with the axisymmetric metric ansatz 



gMNdx^'dx"" = Ly'^^^^^ + fe^^^^y'''^-^^^^^f{y)de^ + 2ib^{y,x)dedx^ 



+a^{y)h^y{y, x)dx^dx'', 



(3.2) 



where capital Latin indices numerate the 6D coordinates while the Greek indices 
are restricted to the 4D coordinates. 
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The evolution equations along the ^/-direction are given by 

^ ^Ff.MF''' -IsypAe-^, (3.3) 



M4 



where = e~^\/f /Lj, Ry is the extrinsic curvature oi y = constant hyper- 
surfaces, K is its 5D trace, ^R^ is the 5D Ricci tensor and ^Df^ is the covariant 
derivative with respect to the 5D metric. Here, /i = /i and 9. The Hamiltonian 
constraint is 



m + -K' = — - 



(^FyMFy^'- -J^^ e-^ - 2{nydy<pf + {d<i>f + l^et 

(3.4) 



and the momentum constraints are 

^D, (k/ - 5yk) = ±-J^,,Fy^nye-^ + (3.5) 

where Uy = e^Lij^pf . 

The Maxwell equations are given by 

Vm (e-^F*^^) = 0, (3.6) 

where Vm is the covariant derivative with respect to the 6D metric. The dilaton 
equation of motion is 

VMV"^+^FV*-^e* = 0. (3.7) 

3.3 Gradient expansion approach 

In this section we will use the gradient expansion method [TTl to solve the 
6D bulk equations. We assume that the length scale i is of the same order of L/. 
The small expansion parameter is the ratio of the bulk curvature scale to the 4D 
intrinsic curvature scale. 
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See section 11.41 and subsection 12.4.21 for more details on the gradient expansion 
method. We expand the various quantities as 

V = KAx) + eh^i^Ky. a;) + ■ ■ ■ , ^ = + e^^^~^ + ■ ■ • , 
= + + ■ ■ ■ , Fye =Fye +eFye +■■■ , 

(0) (1) (0) (1) 

K;=K;+eK;+---, K,' =K,' +eK,' +■■■ . (3.8) 
As to the other quantities, we follow [83] and first assume 

^MI/ = ^l/2^Mj;+... ^ F^. = £F^.+--- , (3.9) 
and then will show that all the 0{e^^'^) quantities in fact vanish. Since d^Ag ~ 

(1/2) 

e^^'^dyAg, we have F^e = e^/^ F^^ H . We will show that this C(e^/2) term 

in F^g also vanishes. See appendix El for more details on the 0{e^^'^) quantities. 
The bulk energy-momentum tensor contains terms like Ff^xF"^^ but these do not 
contribute to the low energy effective theory as they are higher order in the 
gradient expansion. The 5D Ricci tensor is given by 

= e^(^R;:[h]-V,V''i>-V,i^V''i>^+-- - , (3.10) 
%' = -e^ {VxV^i, + VxtpV^tp'^ + ■ ■ ■ , (3.11) 

and ^R^ = 0{e^^'^), where ^ := 4'^^^ — ^- ^^(^ T^tJ. respectively the Ricci 

tensor and the covariant derivative constructed from h^y{x). 

3.3.1 Zeroth order equations 

The 9 component of the Maxwell equations at zeroth order reads 

dy(a^e-'*''-"'+^'"' F^^) =0, (3.12) 
while the equation of motion for the dilaton at zeroth order is given by 

is,(aVe*™9,0'»') + \ (-j^ i;„)'e--*'°>-*™ - ie"-*™-*™ = 0. (3.13) 
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The (fiu) and (66) components of the evolution equations are given respec- 
tively by 



4 VM2£ 



1 (0) 

Fye] e 



4 



2/ 



1 (0) 



, 2A-</,(0)-2^{«) _ 1 2A+0W 

4 VM^r^v 4 



(3.14) 



(3.15) 



and the Hamiltonian constraint becomes 

2 



4/ 



3 ( ^) +^ + 29,^(°) 



1 (0) 



— e 



2x+<j>m 



(3.16) 



The solutions for the above equations are obtained as 



aiy) = y/y, f{y) 



y y^ 



A(x) = -<l>(x), (3.17) 



and 



V;W(y,x) 

(0) 

Fye 



(^(x) + aix) 



y.x) 



Iny — $(a;), 



(3.18) 



where /i and q are integration constants. The momentum constraint imphes 



(1/2) 



9^0" = 0, and therefore a = constant. This immediately leads to F^g= and 
hence F^e = 0{e^^'^). In the following, we put o" = without loss of generahty. 
The 6D metric at the zeroth order is given by 



QMNdx dx 



+ a^{y)h^y{x)dx^dx'' . 



(3.19) 



Then we can see that $(x) is associated with the scaling symmetry Qmn ^ ^^9mn 
and e*^ — > e"^'^ . In fact, we will find that a solution for h^y is given by h^j^ = e'^rj^y 
if the brane preserves the scaling symmetry, where rj^y denotes the 4D Minkowski 
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metric. 



3.3.2 First order equations 

At first order, the (fiu) component of tlie evolution equations is given by 



77' 



-4-/2 



(1) 



2 

y 2/ 



(1) 1 . (1) (1) s 



^ (^i?; - V,V^^ - ^P^^I?'^^) - ^e^'^V^')^; + ^^5;, (3.20) 



wfiere 



1 / (0) W (1) W N 



M4 



(3.21) 



Tfie 4D Ricci tensor i?^ does not depend on y because it is computed from /i^j^ 
whicfi is a function of x'^ only and the index is raised by h^^^. 
The 4D traceless part of Eq. (13.201) is found to be 



(3.22) 



(1) (1), 
where we defined K;^ ■=K;^ -{ll^)bl^K^ and 



r := R ^ - \b::R - ( VJ)""^ - 



(3.23) 



where V^^ := h^'^V^.V^^ and (P$)2 := h'^'V^^V^^. The general solution to 
the above equation is given by 



2^/f 



(3.24) 



where the traceless tensor C^(x) is the integration "constant" to be fixed by the 
boundary conditions. 
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The 4D trace part of the evolution equations is 



77 



-$/2 



^yK. +(,+ 2/ 



(1) 



2 

y 



1 

y 



(3.25) 



and the [Qff] component of the evolution equations is 



-*/2 



(1) 



y f 



2/ 



_1 rp2^ + (p$)21 
2y L ^ ^ J 4 4Lj y 



(3.26) 



The Hamiltonian constraint at first order reduces to 



1 

y 



o/7 -$/2 



1 , M 

22/ 2/ 



(1) 



(1) 



i^^^ +- 



Lj y 



9^,0^.(3.27) 



The dilaton equation of motion at first order reads 



d 



+ - + 



2 

y f 



+ (P$)2] 



1 e-* 



/7„-<I>/2 / (1) (1) 

y 



2L2 y " + ^• 



(3.28) 



Now we define convenient quantities 



1 

2 



(3.29) 



and 



/C : 



3 (1) (1) 
4 



:JJ_ -^12 Ml) 

Li V 2/ 2|/; 



(0) 



(1) 



^/ 2/ 



(3.30) 
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The evolution equations for these variables can be derived using Eqs. fl3.25p - 
(I3.28p . With some manipulation one arrives at 



1 



dy (y'^fic) ^ 



,*/2 



yLi 



(3.31) 
(3.32) 



The two equations have the same structure as that of Eq. (13.221) . The general 
solution for each evolution equation contains one integration "constant" which 
will be determined by the boundary conditions. 

In terms of the above variables, the momentum constraint equations are sim- 
plified to 



(3.33) 



3.4 Junction conditions and effective theory on 
a regularized brane 

Our choice of parameters q implies that f{y) vanishes at i/n and ys- These 
points are conical singularities that are sourced by 3-branes. In order to accom- 
modate usual matter on the branes we need to resolve these singularities. We 
will use the regularization scheme of (721 [92]. The conical branes are replaced 
with cylindrical co dimension-one branes at y = y± and their interiors are filled 
with regular caps. See figure 13.11 for a sketch of the model. The geometry of the 
caps and the central bulk is described by the 6D solutions found in the previous 
section, with different curvature scales L_|_ (L^) for the north (south) cap and Lq 
for the central bulk. 

The action of each brane is taken to be 



^1 



brane 



V(0) + -f/(0)(9^S - eAf,){d^J: - eA^) 



(3.34) 



where g^i> is the induced metric on the 4-brane, V{4>) and U (0) are the couplings 
to the dilaton, and is the Lagrangian of usual matter localized on the brane. 
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ys 




'0'-' 


\ ^ — — 


^^yN 








1 — 'ly y 

^-"■""'^ North cap 



South cap 



Figure 3.1: Schematic representation of the bulk spacetime with two regularized 
caps. 



At this stage we assume that the brane matter does not couple to the dilaton 
field. We introduce a Stueckelberg field S, which is obtained by integrating out 
the massive radial mode of a brane Higgs field. The equation of motion for S 
gives the gradient expansion form of the solution as [83] 



E{9, x)=n9 + c(°) (x) + ec^^^ (x) + ■ ■ ■ 



(3.35) 



where n must be an integer because of the periodicity ^ ^ ^ + 27i. 
The jump conditions for the Maxwell field are 



[[n^FMNe-^]] = -eU{dN^ - eA^) 



(3.36) 



while for the dilaton field we have 



(3.37) 



where •= li^^e^o (-^Ij/b+e ~ -^lyb-J- Here and hereafter in this section all 

the quantities are evaluated at the position of the brane under consideration. The 
Israel conditions are given by 



1 



rp V 

' A (tot 



(3.3^ 
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where 



rri V 
^ /i(tot) 



u 



{dfj:-eAi,){d'T.-eA' 



- W + V' 



(3.39) 



and T~^ represents the matter energy-momentum tensor. 



3.4.1 Zeroth order 



At zeroth order in the gradient expansion the junction conditions fl3.36l) - fl3.38l) 
are written as 



Maxwelh 



-j—Fyeye ' 



ef/(o) (n - eAi°) 



(3.40) 



Dilaton: 



Li y 



n — eA 



(0) 



(3.41) 



Israel (/ii/) : 



v7 ^\ $/2 



M4 



y(o) + if/(0) 

2 



n — eA, 



(3.42) 



Israel (99) : 



V-(o) _ ±[/(o) 

2 e-'f 



n — eAr 



(3.43) 



The above conditions relate several parameters with each other, and the detail of 
the parameter counting of the configuration is found in Ref. [92]. In particular, 
the dilaton jump condition fl3.4ip and the Israel condition 03.431) imply 



2Pf \ 2 



#(0) J 



n — eA 



(0) 



0. 



(3.44) 
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The classical scaling symmetry is preserved by the special choice of the po- 
tentials [561 122] 

V{(f)) = ve'^l'^, U{(f)) = ue-'^l'^. (3.45) 

With these potentials the junction conditions (13. 400 - 03. 431) put no constraints on 
$(x) and Eq. (13.441) is trivially satisfied. In this case the first order analysis will 
provide the equation of motion for $(x), as will be seen in the next subsection. 
In the following, we assume that at the zeroth order, the potentials are given 
by (I325D, that is, f/W(0(o)) = n(o)e-<^"'^/2 v^(o)(0(o)) = t,(o)e<^'"V2_ Then we 

expand the potentials as follows: 

^0) = V^(°n0^°^)+^(v^^^n0^°^) + ^0^^^), (3.46) 

W) = f/^°n0^°^)+e(f/^^n0^°^) + ^0^^^), (3.47) 

where V'^^^ [<\)^^^') and f/'^^)(0(°)) stand for the deviations from the zeroth order 
potentials. 



3.4.2 First order 

The 4D traceless part of the Israel conditions at first order is given by 

1 



T*" 



(3.48) 



where := T^ — [Xj^b^^T^ . The 4D trace part of the Israel conditions reduces 
to 



Q (1) (1) 



4M4 ^ 



1 f/(o) 
AV + — -A 



+ 



n — qA 



(0) 



(1) 



(3.49) 



where we defined 



A\/ = V^«(0(°)) + -t/«(0(°)) 



(3.50) 
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and 



e4'^ + (n- eAf^J ^(^)J . (3.51) 
Using the zeroth order junction conditions, Eq. fl3.49p simply gives 

m] = ihlT.' - 4lAK (3.52) 



The (66) component of the Israel conditions is 



(1) 



Tf ^(0) 
1 



A 



M4 



n — eA 



(0) 



Ar + f/«(0(°)) 



n — eA 



(0) 



(3.53) 



and the dilaton jump condition is 



1 

1 d 



A 



M4 



n — eA 



(0) 



(Ay) - -f/« 
M4 rf* ^ ^ 2 



n — eA 



(0) 



(3.54) 



Using the fact that the zeroth order potential have the scale invariant forms 
(13.451) . the above two conditions are combined to give 



m 



MU^^ ' 2M4 



(3.55) 



Therefore, the momentum constraints become 



{flX^ - e^l^AVS;) = Qt/ - A(Ay) _ 1 e^/^v,^, (3.56) 
In terms of the energy-momentum tensor integrated along the ^-direction, 

T; := 27r£v7e*/%^ (3.57) 



this can be rewritten as 



V,T^ =-T,V,^. 



(3.5^ 
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To fix the integration constants completely, we need the boundary conditions 
at the north and south poles. Near a pole with the coordinate y = yp, where 
p = {N,S}, we have f ^ y — Vp- In order for the evolution equations f l3.22p . 
(I3.3ip . and (13.32^ to be regular at the poles, we require 

k;, /C, J<\y-yp\'^'^0. (3.59) 

Now we can determine all the integration constants included in the general so- 
lutions for K^, /C and J'. Since the structure of the evolution equations and 
boundary conditions are identical for these three variables, we summarize the 
procedure to fix the integration constants in appendix [HI and here we focus on 
the resulting effective theory on the brane. 

Using Eqs. fl3.48p and fl3.52p together with the solution for and /C in terms 
of R and $, we end up with the effective equations 



= 4 [t^; - AF^5;) + ^4 (t7 - AV-6;) , (3.60) 



where the 4D gravitational couplings are defined as 



Vn 



4 := ^ith it = i I Ljydy, (3.61) 



ys 



; denotes a covariant derivative with respect to the induced metric g^^, = a^h^^, 
Rfj,^[q~^] is Ricci tensor computed from q^^ and the potential integrated along the 
^-direction is defined as 

AF = 27ii^e^/'^AV. (3.62) 
The first order equations for J' give the equation of motion for $: 
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For simplicity let us ignore the matter energy-momentum tensor and the po- 
tential on the south brane: T^*^ = AV = 0. In the absence of the (66) com- 
ponent of the energy momentum tensor on the north brane, the 4D effective 
equations can be deduced from the action 



S, 



eflf 



2Ki 



(3.64) 



with the Brans-Dicke parameter wbd = 1/2 (see also appendix B of Ref. [70] 



3.4.3 The exact time-dependent solutions in the 4D effec- 
tive theory 

We now consider cosmological solutions in the 4D effective theory and compare 
them with the known solutions to the full 6D field equations [691 EQl 171] . 

Let us assume = AV = and Tg =0. We consider the case where 
the first order potential is scale invariant form. Then AV^ oc e~*/^ and AV~^ = 
const. := A/k^. We go to the Einstein frame defined by 

V = (3.65) 
and then the equations of motion become 

R;:[h]-ls;:R[h] = -Ae-^*^; + 2$|,$l'^ - <|.|,<|.l^5;, (3.66) 
<l>lJf = -Ae-2*, (3.67) 

where | stands for the covariant derivative with respect to the Einstein frame 
metric h^^. 

Taking /i^j, to be a fiat Friedman- Robert son- Walker metric, h^^^dx^dx^ = 
A^(r)(— (ir^ + c/x^), the equations of motion reduce to 

— = -AA^e-^" - h'^ (-X = -AA^e^^" + -<l>'\ (3.68) 

^3 3V^/3 3 ^ ' 

and 

$" + 24^$' = AA^e-^*, (3.69) 
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where ' := d/dr. For A = a solution of the above equations is 

A\t) = A^T + A2, $(r) = ±73 In A{t) + A3, (3.70) 

where Ai [i = 1,2,3) are integration constants. For A 7^ a solution is 

A{t) = $(r) = InA(r), (3.71) 

where = A/2. This indicates that q^^ = e^rj^^y which is expected from the 
scaling symmetry. 

The brane scale factor Ob and the "radion" \E' are given by 

ab = e-*/M, * = e*. (3.72) 

The solution (13.701) gives the same 4D observables (the brane scale factor and 
radion) as an exact 6D solution found by Copeland and Seto (equation (70) 
of [H]). In the same way the solution fl3.7ip reproduces the 4D quantities of 
their equation (76). (The latter solution was first found by ToUey et al. |69j.) 
Thus we show that the solutions to the full 6D equations are reproduced by our 
4D effective theory on a regularized brane with the scale invariant potential and 
with or without additional "tension" A. 

At the zeroth order, the amplitudes v^^'^ and m'-"-* of the potentials are fine- 
tuned. If we change the amplitude of the scale invariant potentials, which is 
equivalent to adding a cosmological constant in the 4D effective theory, we get a 
runaway potential for $ and the 4D spacetime becomes non-static. 

3.4.4 Breaking the scale invariance 

Finally, let us consider the case where the first order potential breaks the scale 
invariance. As the BD parameter is given by 1/2, this model violates the con- 
straints coming from the solar system experiments unless the BD scalar $ is 
stabilized. It is suggested that the potential on a brane can naturally stabilize 
the modulus. For example, if we consider potentials Vi{(j)^^^) = v^^'^ e'^'^^°^ and 
f/i(0'^°)) = u^^'^e^^'^°\ the effective potential is given by 

AF = 2vr^v^e*/2 (^^e"^* + l^^^-^'^''^" - e^"^)^) • (3-73) 
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As we saw in the previous subsection, if we take s = 1/2 and t = —1/2, AV is 
independent of $. However, in general, it is possible to have a potential with 
a minimum by choosing s,t,v^^^ and u^^^ appropriately [92]. Then $ can be 
stabilized and general relativity (GR) is recovered. 

3.5 Conclusions 

In this chapter, we derived the low energy effective theory in the six-dimensional 
supergravity with resolved 4-branes. The gradient expansion method is used to 
solve the bulk geometry. The resultant effective theory is a Brans-Dicke theory 
with the Brans-Dicke parameter given by ubb = 1/2. If we choose the dilaton 
potentials on the branes so that they keep the scaling symmetry in the bulk and if 
we tune their amplitudes then there is no potential in the effective theory and the 
modulus is massless. Thus the static four-dimensional spacetime has vanishing 
cosmological constant. It is also possible to obtain time-dependent solutions due 
to the dynamics of the modulus field and we showed that they are identified with 
the six-dimensional exact time dependent solutions found in [7Tj . 

Even if the potentials preserve the scaling symmetry, it was found that there 
appears an effective cosmological constant in the four-dimensional effective theory 
by changing the amplitude of the potentials. Then in the Einstein frame, the 
modulus field acquires an exponential potential and the static solution is no 
longer allowed. Again, we showed that the cosmological solutions obtained in the 
effective theory can be identified with the six-dimensional exact time dependent 
solutions found in [691 HOI Ej • 

Our effective theory allows us to discuss cosmology with arbitrary matter on 
the brane. As the BD parameter is given by 1/2, it is impossible to reproduce 
realistic cosmology without stabilizing the modulus field. As it was suggested by 
Ref. [92], it is easy to generate a potential for the modulus $ with a minimum 
by breaking the scaling symmetry from the dilaton potentials on the branes. 
Then it is possible to reproduce GR at low energies. However, once we stabilize 
the modulus, the cosmological constant on a brane curves the four- dimensional 
spacetime in the same way as in GR. 

Our result would indicate that it is possible to reproduce sensible cosmology in 
this six-dimensional supergravity model at low energies but it would be difficult 
to address the classical part of cosmological constant problem in this set-up. 
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However, we should mention that our effective theory is vahd only up to the 
energy scale determined by inverse of the size of extra-dimensions. This condition 
is roughly given by Hi^ < 1 where H is the Hubble parameter. If we consider 
scales smaller than gravity becomes six- dimensional and from the table-top 
experiments, i^, is smaller than a few /im. Then for H > 10^^ eV, our universe 
becomes six- dimensional and it is impossible to use the four-dimensional effective 
theory. In order to address the behaviour of the universe at high energies, we 
should deal with time-dependent solutions directly in six-dimensional spacetime. 
This remains an open question. 

Another issue that we did not touch in this chapter is the quantum part of the 
cosmological constant problem, that is the stability to quantum corrections. As 
mentioned in the introduction, there are reasons for believing that the combina- 
tion of the bulk supersymmetry and scaling symmetry can maintain the quantum 
corrections to be of the order of the cosmological constant [62|. In order to ad- 
dress this issue, it is necessary to perform 6D calculations which take into account 
the loops of the KK modes. 

Finally, we briefly comment on the limit where the codimension one branes are 
shrunk to codimension two objects. Our effective theory shows no pathological 
behaviour in this limit as long as the four-dimensional energy-momentum tensor 

integrated along the ^-direction remains finite. However, in this limit, the first 

(1) 

order extrinsic curvature K/iu diverges and then the first order correction to the 
four-dimensional metric diverges. Then it is not clear whether there is a physical 
meaning in this limit. This is related to a deep issue of whether it is possible 
to put ordinary matter on codimension 2 objects [93l [94] and we also leave this 
problem as an open question. 



Chapter 4 



On the 4D effective theory in 
warped compactifications with 
fluxes and branes 

4.1 Introduction 

In this chapter, using the gradient expansion method, we shall derive the 4D 
effective theory for warped compactifications with fiuxes and branes in the lOD 
type IIB supergravity. 

Warped string compactifications with fiuxes and branes have provided a novel 
approach to long standing problems of particle physics and cosmology. Inspired by 
the earher work by Randall and Sundrum [1], Giddings et al (GKP) [13] showed 
that, in type IIB string theory, it is possible to realize the warped compactifica- 
tions that can accommodate the large hierarchy between the electroweak scale 
and the Planck scale. In the GKP model, the warping of the extra-dimensions is 
generated by the fiuxes and the presence of D-branes and orientifold planes. All 
moduli fields are stabilized due to fiuxes except for the universal Kahler modulus. 
The warped compactifications also provide a promising background for cosmo- 
logical infiation. An infiaton can be identified with the internal coordinate of a 
mobile D3 brane moving in a warped throat region [95] . 

So far, most works are essentially based on effective 4D theories derived by a 
dimensional reduction. However, it is a non-trivial problem to derive the 4D ef- 
fective theories including the warping and branes. In fact there have been debates 
on the validity of the derivation of the 4D effective theories. In a conventional 
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Kaluza-Klein theory, the effective theory can be derived by assuming a factor- 
izable ansatz for the higher-dimensional fields. This approach has been applied 
even in the presence of the warping. For example, the lOD metric ansatz of the 
form 

dslo = h~^/\y)e-^<^^g^,{x)dx^dx'' + h^/\y)e^^^^h^ndy'^dy'' , (4.1) 

is often used where x denotes the coordinates of the 4D non-compact spacetime 
and y denotes the coordinates of the 6D compact manifold. The function u{x) 
is identified with the universal Kahler modulus. However, it has been criticized 
that the higher- dimensional dynamics do not satisfy this ansatz [961 EZj- Thus 
it is required to derive the 4D effective theory by starting from a correct ansatz 
and consistently solving the lOD equations of motion. 

These attempts were initiated in Refs. [iQl [521 [16] . Ref [ID] derived an exact 
time-dependent lOD solution that describes the instability of the warped com- 
pactification due to the non-stabilized Kahler modulus. It was shown that this 
dynamical solution cannot be described by the metric ansatz (14.11) . Ref. [52] de- 
rived the potential for the moduli fields by consistently solving the lOD equations 
of motion using the metric ansatz that is consistent with the dynamical solutions 
in the lOD theory. It was found that the universal Kahler modulus is not a sim- 
ple scaling of the internal metric as is assumed in Eq. (14. ip . A similar consistent 
ansatz was proposed in Ref [16] and the 4D effective theory without potentials 
was derived. 

In this chapter, we build on these works to present a systematic way to derive 
the 4D effective theory by consistently solving the lOD equations of motion of 
type IIB supergravity. We exploit the so-called gradient expansion method, which 
has been shown in the previous chapters to be a powerful method to derive the 4D 
effective theory in the context of 5D/6D brane world models [I2l HH HH [HHl l98] . 
Our method also provides a scheme to study the gravitational backreaction to 
the warped geometry due to the moduli dynamics, branes and fluxes. 

4.2 lOD equations of motion 

Let us start by describing the type IIB supergravity based on Ref. [13]. In the 
Einstein frame, the bosonic part of the action for the type IIB supergravity is 
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given by 

Sub = [d^'xV^o(^''^ ^^'^'^''^ ''''' 

2Kfo J \ 



2 (Imr)' 12Imr 4 ■ 5! 

C(4) A (7(3) A G(3) 



^K^Q ./ Imr 



(4.2) 



where the combined 3- form flux is (7(3) = F(3) — 'r_f/'(3) with F(3) = dG{2)^ -f^(3) = 
(ii?(2) and 

r = C(o)+2e-^ (4.3) 

where is the so-called dilaton, C(j) is the Ramond-Ramond potential of rank (j) 
and -8(2) is the NS-NS potential. The bar denotes the complex conjugate. The 
5-form field is given by 

1 1 

-^(5) = -^(5) - 2<^(2) A -^(3) + 2^(2) A -F(3), (4.4) 

with F(5) = (iC(4). The total action in our model is 

S = Sub + Sloe, (4.5) 



where the term Sioc is the action for the localized sources such as D3-branes and 
03 planes; 

Sloe = (--"3 j d^'x^, +T;j d'xG,^ , (4.6) 

where the integrals are calculated over the 4D non-compact space at the point j 
in the compact space and gj is the determinant of the induced metric on a brane 
at the point j (/i3 is positive/negative for D3-branes/03 planes). 

We can derive the equations of motion for the fields from the action (14. 2p . 
The trace reversed Einstein equations are given by 



Re{dATdBr) ^ 1 

J^AB — — — -2 



Re {GacdGq^^) — —GcDEG'^^^'^^^^gAB 



2 (Imr)^ 41mr 
+ ^Fap,...p,F/-^^ + nl, (t^^- - l^'^^gABT^-^ , (4.7) 
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where T^°^ is the energy-momentum tensor for the locahzed sources defined by 

and ^^^^gAB is the metric for the lOD spacetime. 

Following Refs. [IHl [52], we take the lOD metric as 



dS^Q = h y)g^y{x, y)dx^dx'' + (x, y)-^„,n{y)dy"'dy'', (4.9) 

with 

hix,y) = hiy) + ho{x), (4.10) 

where upper case Latin indices run from to 9, Greek indices run from to 3 
(non-compact dimensions) and lower case Latin indices run from 4 to 9 (compact 
dimensions). In the following, all indices are raised by Qfj^i, and 7mn- The function 
ho{x) is the so-called universal Kahler modulus. It should be emphasized that 
this is not a simple scaling of the internal metric. 

Using our metric ansatz, we can calculate the lOD Ricci tensor. The mixed 
component is calculated as 

= -a'^'^K^ppia + - Ih-'h^Kp, (4.11) 

where we defined ^ 

Kfjipp = ~g^iJ,p^ Kp = Kfj^i/p, (4.12) 

and I denotes the covariant derivative with respect to (7^,^, that is, 

K T : — K : — Wr°" K r — Wr°" K ^4 1 '^'l 

where ^^-"r^^ is the Christoffel symbol constructed from g^^. The non-compact 
components are given by 



+h-'Kj - hi~^h^K,g,, - h-'K^,K, + 2h-^K,\Kj, (4.14) 
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where ; denotes covariant derivative with respect to ■jab, that is, 

where ^^^F^^, is the Christoffel symbol constructed from jat- The compact com- 
ponents of the Ricci tensor are 

Rab = ^^^Rab{l) + ^h~'^h.dh''^Jab " ^h~^ h'.^lab " ^h^llab " ^h~'^h.ah;b 

+Ka;b - {htKa + haKb) + ^k-' h''' K.Jab ' K'^^K^fSb- (4.16) 

Note that the self duality of the 5-form field must be imposed by hand 

^(5) = *^{5)- (4.17) 
We shall take the self-dual 5-form field in the form 

F(5) = (1 + *) ^/^dya{x, y) A A dx^ A dx"^ A dx^, (4.18) 

Then the self duality condition fl4.17p is automatically satisfied and the only 
non-zero components of F(5) are 

^{5)afio---fi3 = Oi-a^ ^5^g...^3, (4.19) 
^(5K...„,5 = -/l'«''=^'^£cni...n5, (4.20) 

where ^^''e^Q...^3 denotes the Levi-Civita tensor with respect to g^,^ and ^^^ecm-ns 
denotes the Levi-Civita tensor with respect to jrrm- 

4.3 Gradient expansion method 

In this section, we will use the gradient expansion method to solve the lOD 
Einstein equations (14.71) . 



CHAPTER 4. THE 4D EFFECTIVE THEORY IN TYPE IIP SUGRA 



77 



4.3.1 Gradient expansion 

The gradient expansion is based on the assumption that x derivatives are sup- 
pressed compared with y derivatives 

dl<^dl (4.21) 

See section 11.41 and subsection 12.4.21 for more details on the gradient expansion 
method. Using this assumption, we can reduce the partial differential equations 
with respect to x and y to a. set of ordinary differential equations with respect to 
X. We expand the metric as 

g,,u{x, y) =i°l (x, y)+ (x, 2/) + ■ ■ ■ , (4.22) 

where the first order quantities are of the order d1/dy. Accordingly, K^^^p is also 
expanded as 

(0) (1) 

K =K^up ~l~ K^fiup (4.23) 

4.3.2 First order equations 

We assume that at zeroth-order, the axion/dilaton is constant and the 3-form 
field only has non-zero components in the compact space, i.e. 

G(3) = l^GabMdy" A dy' A dy'^. (4.24) 

Furthermore, we assume that the zeroth-order metric is independent of the com- 
pact coordinates 

g!!^ {x,y) =g)!u (x). (4.25) 
The Bianchi identity/equation of motion for the 5-form field becomes 

ail + 2h-\,a'' = 2/1"^ + 2Kl,h-'npi"^, (4.26) 

121mr 

at zeroth order where *6 is the Hodge dual with respect to Qab and we defined a 
rescaled D3 charge density p^"^ which does not depend on h. We define u as 

00 = a- h-\ (4.27) 
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Then the Bianchi identity is rewritten as 

The non-compact components of the Einstein equations (14.71) can be rewritten as 



- g^u h\^ \ + h K^^^ --h h K^g^, g^. 



t HG(3) - *6G(3) r + ^ Z {T.pT - /i3(y)) ,(4.29) 



96Imr "'^'^ ' " 2 

where '^^^R^u denotes the Ricci tensor constructed from g^y, l^siv) = l^s^iy — 
yi)/y/l which is independent of h and we dropped the non-hnear term in u;. In 
the same way, the compact equations can be rewritten as 

1 |. (1) 1 1 / 1 1 w 

= ^hul^^'yab + ^ {h-aUJ-b + h-b^^.a) 

r - 1 - i - ' 

~^4Imr {^acdGb ) " Y^GcdeG" "'jab - Y^Gpgr *6 G^'^^'^Jab 

+4o - IhhabT^'^' - h-'Tsp^riab^ , (4.30) 

where ^^^Rab denotes the Ricci tensor constructed from 'jab- 
The GKP solution is obtained by taking 

*6G(3)=^G(3), c^ = 0, (4.31) 

with local sources that satisfy fisi^y) = T-sp'f'^. With these conditions, it is straight- 
forward to show that the lOD equations motion are satisfied. Then the warp factor 
is determined by 

~ ^ = li^^'''^"'^ + 2«:?oT3p|,-. (4.32) 

Our strategy is to assume that to, ^^^Rab and the right hand sides of Eqs. (14.291) 
and (I4.30p are first order in the gradient expansion and take into account these 
contributions as a source for the 4D dynamics of gp.v{x) and ho{x). 
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4.4 4D Effective equations 

In this section we solve the lOD Einstein equations to get the 4D effective Einstein 
equations. Hereafter we omit the superscript (z) that denotes the order of the 
gradient expansion. A key point is the consistency condition that imphes that the 
integration of the total derivative term over the 6D internal dimension vanishes 

J d%^v,, = 0, (4.33) 

for an arbitrary v. This condition provides the boundary conditions for the lOD 
gravitational fields and yields the 4D effective equations. 
First let us take the trace of Eq. 04.291) : 

(^)i? + [h-^K, + u,:)'' = HG(3) - *6G(3)r + 2Kih-^ {np'r - My)) ■ 

(4.34) 

Then integrating Eq. 04.341) over the internal manifold, we get 



+2^ / d%Vlh^' {Tspi"' - /^3(y)) , (4.35) 



'(6) 

where V(6) is the volume of the internal space 



Vie) = J d'y^. (4.36) 

On the other hand, by combining Eq. 04.29P with Eq. 04.341) . we obtain the 
traceless part of the equation 

^'^R^,-^g^,^'^R-h-' (^h\^,-^g^,h\'^+h-' (^K^^h- =0. (4.37) 

Integrating this over the compact space we obtain 

'^'^R.u - \9,J'^R = (^h\,, - \9,.h\'^ , (4.38) 



CHAPTER 4. THE 4D EFFECTIVE THEORY IN TYPE IIP SUGRA 



80 



where the function H{x) is defined as 

H{x) = ho{x) + C, (4.39) 

where 

C^-^ [ d%^h{y). (4.40) 

^6) J 

Here hi{y) is obtained by solving Eq. fl4.32p . Finally, we need an equation that 
determines the dynamics of ho{x). Let us calculate the trace of Eq. fl4.30p 

|^G(3) - *6G(3)r + ^4oh-' M - T,pt) . (4.41) 
Combining Eq. (14.411) with Eq. (I4.34p we get 

(6)^ + _ lh\^ +(-K,- hu.^ = UG(3) - *6G(3) I' 

+2t,l^h-^ (/i3 - T^pt) . (4.42) 

Then integrating this over the compact space we obtain the equation of motion 
for H{x) 



I d%^h-'{pM-Tzpt)- (4.43) 



3% 

4.5 The 4D effective theory 

The effective 4D equations are summarized as 

(^)G^, = H-^ {H\^, - g,,H\l - Vg,,) , (4.44) 
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where the potential V{H) is given by 

They can be deduced from the following 4D effective action 

Seff = ^J d'x^ [H^'^R{g) - 2V{H)] , (4.47) 

where k1 = ft;io/^(6)) which can be determined by integrating the lOD action over 
the six internal dimensions. 

Performing the conformal transformation g^^, = H^^f^^ we can write the 
previous 4D action in the 4D Einstein frame as 



R{f,.) - \{V\nHf - 2V{H)H-^ 



(4.48) 



where now V denotes the covariant derivative with respect to the Einstein frame 
metric f^j^^. The 4D effective equations of motion in the Einstein frame are given 
by 

RAf) = \ (V;.lni?) {VMH) + V{H)H-^f^,, (4.49) 

4 9 dV 

V„V" Ini/ = -lv{H)H-' + ^^H-\ (4.50) 
3 6 dH 

By defining p{x) = iH{x), the kinetic term can be rewritten into a familiar form 

(4.51) 



1 

M 



SE,kin — 7^~^ / d Xy —f 



dnod^p 



IP- PI 



We would find the same result for the kinetic term for the universal Kahler 
modulus even if the wrong metric ansatz Eq. (14. ip was used to perform a dimen- 
sional reduction where p(x) = ie^'^^^\ In a region where the warping is negligible 
hi{y) ^ ho{x), H{x) can be identified as e^^^^\ which is the simple scaling of 
the internal metric. However, in a region where the warping is not negligible, the 
original lOD dynamics is completely different between (14.11) and (14.91) [15] . 

The potential associated with the 3-form agrees with the result obtained in 
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Ref. [52]. It should be emphasized that this potential is positive-definite. Nev- 
ertheless, the 3-form contribution to the 4D Ricci scalar ^^^R is negative definite 
[961 [97] in accordance with the no-go theorem for getting de Sitter spacetime 
[99] . The resolution is the kinetic term of the modulus H{x), that is, ^^■'i? is not 
directly related to V [52]. In fact the 4D Ricci scalar is related to the potential 
as 



At the minimum of the potential we have ^^^R = V = dV/dH = 0, but, if we 
move away from the minimum and the modulus H{x) is moving, the potential 



We also notice that D3 branes with fi3{y) = T^p^^'^ do not give any gravi- 
tational energy in the 4D effective theory. On the other hand, anti-D3 branes 
with fis^y) = — Tapg*"^ give a potential energy. This was used to realize de Sitter 
vacuum [lUUj . 

4.6 Discussions 

In this chapter, we presented a systematic way to derive the 4D effective theory 
by consistently solving the type IIB supergravity equations of motion in warped 
compactifications with fiuxes and branes. We used the gradient expansion method 
to solve the lOD equations of motion. The consistency condition that the inte- 
gration of the total derivative terms over the internal 6D space vanishes gives the 
boundary conditions for the lOD gravitational fields. These boundary conditions 
give the 4D effective equations. Once the solutions for the 4D effective equations 
are obtained, we can determine the backreaction to the lOD geometry by solving 



We did not introduce the stabilization mechanism for the modulus H{x). It is 
essential to stabilize this universal Kahler modulus to get a viable phenomenology 
[lUU] . Usually, non-perturbative effects are assumed to give a potential to this 
moduli. The non-perturbative effects will modify the lOD dynamics and the 4D 
effective potential have to be consistent with this lOD dynamics. Most works so 
far introduce the non-perturbative potentials for the modulus field directly in the 
4D effective theory and it is not clear this is consistent with the original lOD dy- 
namics. It is desirable to derive the potential for the moduli fields by consistently 
solving the lOD Einstein equation with non-perturbative corrections. In addition. 




(4.52) 



cannot be read off from ^^^R. 



Eqs. flCTD. fl07|) and fl02D . 
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a mobile D3 brane plays a central role to realize inflation [95]. The D3 brane with 
/i3 = T3 probes a no-scale compactification and the brane can sit at any point of 
the compact space with no energy cost. In fact, the D3 brane does not give any 
gravitational energy in the 4D effective theory. However, once the stabilization 
mechanism is included, this is no longer true. This is in fact an important effect 
which generally yields a potential for the D3 brane that is not enough fiat for 
slow roll. It was also pointed out that the gravitational backreaction of the D3 
brane is essential to calculate the corrections to the potential |101j . 

The method presented in this chapter can be extended to include the non- 
perturbative effects by introducing an effective lOD energy-momentum tensor in 
the lOD Einstein equations. It is also possible to include a moving D3 brane in our 
scheme along the line of Ref . [HH 1102] , which studied the dynamics of D-branes 
with self-gravity in a 5D toy model. Then we can calculate the potential for a 
mobile D3 brane by taking into account the stabilization and the backreaction 
of the D3 brane. In addition, it is essential to study how the motion of the D3 
brane is coupled to 4D gravity in order to address the dynamics of inflation. 
These issues are left as open questions. 



Chapter 5 

Introduction to brane inflation 



So far we have seen how to obtain the 4D effective theories of gravity from the 
more fundamental higher-dimensional theories. It is now important to consider 
applications of these 4D effective theories in cosmology. Specifically we will dis- 
cuss the application of brane-world models to inflationary cosmology know as 
brane inflation. 

This chapter is meant to give a brief introduction to the idea of brane inflation. 
There are several good reviews of brane inflation in the literature, see [H I103[ 
II U4] . We shall discuss the main features of brane inflation, in particular the so- 
called DBI-inflation model, where DBI stands for Dirac-Born-Infeld. The DBI- 
inflation model is rather interesting because the particular and unusual form of 
the kinetic energy of the inflaton can give rise to novel observational signatures. 
For example, it can produce large non-Gaussianities in the distribution of the 
primordial curvature perturbation. The next two chapters will be dedicated to 
the study of non-Gaussianities in rather general models of inflation. The models 
considered will include the DBI-inflation model as a particular case, which was 
our main motivation to study models with non-standard kinetic energy. 

Before introducing the general ideas of brane inflation, in the next section we 
shall briefly introduce a much older model of standard inflation and explain why 
is inflation required in the first place. 

5.1 Standard inflation 

Inflation is a period of accelerated expansion in the very early universe. Dur- 
ing this period the scale factor grew at least thirty orders of magnitude. This 
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extremely fast and brief expansion is needed to solve or aleviate the fine tuning 
problems of the standard big bang theory. For nice reviews on inflation see for 
example [TnHlfTnH] . 

For instance, within the standard big bang model there is no explanation 
for why is the cosmic microwave background radiation (CMBR) so homogeneous 
and isotropic. In fact, regions of the CMBR sky separated by more than one 
degree would have been causally disconnected at the time when the photons 
were emitted. This is called the horizon problem. Inflation solves this problem 
by increasing the particle horizon so that all regions of the present observable 
universe were in causal contact in the past. 

From observations we know that today the ratio of the energy density to the 
critical density is within a few percent of unity [107] . But during the standard 
big bang theory this ratio is driven away from one. This implies that in the early 
universe this ratio would have to be extremely fine tuned to unity. This is the 
flatness problem. The inflationary solution is simple. Because of the accelerated 
expansion, the total energy density is driven towards the critical density to an 
arbitrary precision depending on the duration of inflation. Inflation was originally 
proposed to solve these and other problems of the standard big bang theory [108j . 

Later it was realized that inflation can also explain much more. For example, it 
can be used to explain the observed fluctuations in the CMBR. Inflation increases 
the wavelength of the fluctuations of any light scalar field. When these cross 
the horizon they become classical perturbations that will seed the primordial 
perturbations that will eventually grow into large scale structures, like clusters 
of galaxies that we see today. When the inflationary predictions are compared to 
observations they agree quite well, providing us strong evidence for inflation. 

Accelerated expansion, i.e. a > 0, imphes that the equation of state parameter 
of the fluid is 

c^<-i. (5.1) 

To have quasi-exponential expansion we further require that the Hubble param- 
eter varies slowly, i.e. 

e = « 1, (5.2) 

this implies uj ^ —1, that is the cosmological constant case, e is called the slow- 
roll parameter. In order to inflation to last long enough to solve the horizon and 
flatness problems we need the inflaton field to accelerate slowly. Equivalently, we 



CHAPTER 5. INTRODUCTION TO BRANE INFLATION 



86 



define another slow-roll parameter like 



and we require it to be small. </> denotes the scalar field that is sustaining infiation. 
For a standard kinetic term, single field infiation model with potential V{(f)) the 
previous two conditions are equivalent to require that the slope and the curvature 
of the potential are small. 

Usually infiation is caused by one (or more) scalar field as just described. One 
important open question in infiation is the origin of these scalar fields. Using a 
more fundamental theory, like string theory, brane infiation tries to address this 
problem by finding natural candidates for the infiaton. Natural candidates are 
the scalar fields that describe the positions of a brane in the higher-dimensional 
spacetime. In the next section, we shall introduce the main ideas of brane infia- 
tion. 

5.2 Brane inflation 

The first brane infiation model was proposed in 1998 by Dvali and Tye |109j . They 
pointed out that the interaction energy between a D3-brane and its antibrane (a 
D3-brane) in a higher-dimensional spacetime can give rise to infiation. In this 
model, infiation ends when the two branes collide, annihilating each other. The 
energy released by the annihilation is used to reheat the universe that then follows 
the usual big bang model. Figure ISTTl is a schematic representation of this scenario. 
There are two (or more) branes separated by a distance r. This distance r plays 
the role of the infiaton in the 4D effective theory. When the distance between 
the branes is of order l/Mg, where Mg is the string mass scale, the field becomes 
tachyonic and infiation ends. From a 4D viewpoint brane- antibrane infiation is 
very similar to hybrid infiation. 

For a large separation, the interaction energy between the D3 / D3 pair is [95] 




where T3 is the tension of the D3-brane, is the lOD Planck mass and it is 





(5.4) 
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Antibrane 



Brane 




Annihilation into radiation 
Figure 5.1: Brane-antibrane inflation and reheating. 



given by 



Mf,o) = ^, (5.5) 

where Mpi is the 4D Planck mass and is the volume of the compact manifold. 

There is a problem with this model/potential. If one defines the usual slow-roll 
parameters as 

(5.6) 



V 



77 ' V 



V" 



where a prime denotes derivative with respect to the canonically normalized field 
(j), given by </) = \/T^r, then to have long lasting inflation we need £,?7 <C 1. For 
the potential (15. 4p 

(5.7) 



7] 



-0.3 ( - 

r 



and ?7 ^ 1 is only possible if r > L. To have r > L is impossible in a manifold 
of size L. One possible and natural way to evade this //-problem is to consider 
brane inflation in a warped geometry [95]. We will discuss such models in the 
next section. 



5.3 The KKLMMT model: slow- moving brane 

In this section, we present the Kachru et al. (KKLMMT) model [95] of warped 
brane inflation. They consider a D3-brane and an anti--D3-brane in a 5D anti- 
de-Sitter space {AdS^). In [95], the authors also provide a concrete example of 
warped brane inflation in string theory. In that case the AdS^ is replaced by 
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distance = 
inflaton,(|) 




-e-anti D3 brane 



Figure 5.2: D3 and D3 branes in the KKLMMT model. From [T]. 

the Klebanov-Strassler geometry [39] , but the latter can be well approximated by 
AdS^ X 5*^ in the region far away from the tip of the geometry and this is a region 
of interest to study inflation. Thus, we shall present the AdS^ case because it 
captures the essential features despite being simpler. 

The Z^S-brane is fixed at the infrared end of the geometry (region with small 
values of the warp factor). The D3-brane is mobile and experiences a small 
attractive force towards the D3-brane as in figure 15. 2[ If the D3-brane(s) were 
to be absent then the DS-hiane would feel no force in this background as the 
electrostatic repulsion coming from the 5-form background flux exactly cancels 
the gravitational attraction. 

They consider a string theory compactification on AdS^ x S^. It is a solution 
of the lOD supergravity with a 5-form field strength F^. The metric is 



where the warp factor h is h = R^/r^. R is the characteristic length scale of the 
AdS^ spacetime and it is related to the 5-form charge. We choose to truncate 
the AdS^ to the region where tq < r < r^ax- From (15. 8p one sees that the region 
of small r is the bottom of the gravitational well, this justifies the assumption 
of placing the i53-brane(s) at the infrared end of the geometry. We denote the 
position of the Z)3-brane by tq. The background 5-form field is 



ds'^ = h ^Tj^ydx^dx" + h^dr'^, 



(5.8) 



(5.9) 
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and the 4-form gauge potential C4 is 



= ^. (5.10) 

The location of the D3-brane is denoted by ri. The motion of the D3-brane is 
described by the Born-Infeld action plus the Chern-Simons action as 




S = -T, I d^xv^^^/l - ^gf^-^d^rAn + T, I d^x{C,)^, , (5.11) 



where g'^'^ is the inverse metric along the non-compact directions. T3 is the D3- 
brane tension. The action for a i53-brane is similar to the previous action with 
the sign of the second term changed, because a D3-brane has opposite 5-form 
charge. Now it is easy to see that if the D3-brane is slow-moving and there are 
no i53-branes then the two contributions in fl5.1ip cancel each other partially to 
give 

S = nj d'x^^g^-'d.nd^r,. (5.12) 

This means that at low energies the DS-hiane behaves as a free field. 

Now, if we introduce a single D3-brane this will perturb the metric and the 
F5. The D3-brane will acquire a potential energy as [95] 

V = 2T,^Jl-l'I], (5.13) 



where is the number of D3-branes that are present elsewhere in the manifold 
and create the warping. The previous expression is valid in the limit ri ^ tq. The 
first term is a constant potential energy due to the presence of the i?3-brane, but 
the effective tension of the antibrane is reduced by a factor r^/R^. The second 
term is location dependent and the minus sign indicates attraction between the 
branes. It is slow varying as the inverse fourth power of the radial location of 
the D3-brane. This term is also suppressed relatively to the first constant term. 
Then we can naturally get e,r] <^ 1. 

This model is not totally free of problems, in particular a version of the so- 
called ?7-problem discussed in the previous section is still present due to coupling 
to gravity, see [95] for details. 
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5.4 The DBI model: "relativist ic" regime 

We have seen in the previous section that if the brane is moving slowly ("slow- 
roll" regime) and if the potential is flat enough then inflation occurs. References 
|110l lllll 1112^ I113j noted that inflation is also possible in a different regime of 
the parameter space. They realized that inflation is possible even if the potential 
is steep and the D3-brane is "relativistic" . In their model, the warped geometry 
slows down the inflaton field even on a steep potential. They called this model 
DBI-inflation because the action that governs the movement of the brane has the 
form of the DBI action. It reads 



where /(0) is related to the warp factor and is the rescaled position of the D3- 
brane in the throat. V{(f)) is the potential energy that arises from interactions 
with the background fields. For an AdS-\ike throat, /(0) oc (f>~^. 

There are two different situations that have been studied in the literature. 
The ultraviolet (UV) model [1101 llllj . where the inflaton moves from the UV 
side of the throat to the tip of the geometry (the infrared (IR) side), under the 
potential 



where A is a parameter that depends on the background fluxes. 

The other model that has been studied is the IR model |112[ 1113] , where the 
brane moves away from the tip of the throat towards the Calabi-Yau part of the 
manifold. In this case the potential is 





Mpi 



(5.15) 




m ~ 



H. 



(5.16) 



In the spatially homogeneous case the Lagrangian is 



P{x, 0) = -/(0)-i y/1 - 2Xf{4>) + /(0)-i - v{4>), 



(5.17) 



where X is defined a.s X = —1/29^0(9^0 = 0^/2. The requirement that the 
argument of the square root in the Lagrangian has to be non-negative gives a 



CHAPTER 5. INTRODUCTION TO BRANE INFLATION 



91 



speed limit for the inflaton as 

1 - > 0. (5.18) 
We can define the speed of sound 

If the inflaton's speed 

is close to saturating the previous bound its speed is said to be relativistic. In 
this case, Cg <^ I and ^ ±1/a/ The requirement that e = —jp <^ 
1 implies that the potential energy dominates over the kinetic energy during 
inflation despite the fact that the inflaton is moving relativistically. 

In the next chapter, we shall study linear and non-linear perturbations of gen- 
eral models of inflation that include DBI-inflation as a particular case. We will be 
interested in calculating the higher-order statistics of the curvature perturbation. 



5.5 Multiple field DBI-infiation 

So far in this introductory chapter we have only discussed single field models of 
brane inflation. The inflaton field was identified with the position of a D3-brane 
in a higher-dimensional spacetime, like an AdS^ throat for example. However, 
in more realistic string theory compactifications [131 ES] the throat is a 6D space 
that possesses a radial direction but also five angular coordinates. In general, one 
expects a probe D3-brane to move in all these directions. From a 4D effective 
theory point of view, each of these positions of the brane in the compact space is 
described by a 4D scalar field. Thus, DBI-inflation is naturally a multiple field 
inflationary model. 

Some cosmological implications of including angular momentum in DEI brane 
inflation have been considered in [114[ I115j . Due to some confusion in recent lit- 
erature about the correct action for multi-field DBI-inflation here we shall outline 
the derivation of the multi-field DEI action |116l I117j . 

We consider a D3-brane with tension T3 moving in a lOD warped string com- 
pactification. The lOD metric is 

ds^ = HABdY^dY^ = h-^{y^)g^^dx^'dx'' + h^y^)Gij{y^)dy^dy^ , (5.19) 

where h is the warp factor, (yf^,^ and x'^ denote the metric and the coordinates of the 
4D non-compact space respectively. Gjj and y^ are the metric and the coordinates 
of the 6D internal space respectively. The lOD coordinates are denoted by Y"^ = 
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{x'^,?/^}. The dynamics of the DS-brane is described by the DBI Lagrangian 

P = -n^-det^^,, (5.20) 

where det denotes the determinant and 7^,^ is the induced metric on the brane. 
The induced metric is given by 

= HabO^Y^P^Y;^), (5.21) 

where Y^-^^x'^) is the brane embedding. They are given by Y^-^ = [x^ ^ tp^ {x'^)) . 
Finally the induced metric on the brane can be written as 

= {g,,u + hCijd^^^dy^-^) . (5.22) 
This implies that the DBI Lagrangian is 

p = -Tg/i-^v^x/det {Si^ + hGijd^'^^d^^J). (5.23) 
Introducing the rescaled quantities 

/ = </'' = V^/, (5.24) 

and a potential V (0^) , the DBI Lagrangian is 

P = -fi<p'r' (v^ - 1) - V{^'), (5.25) 

where 

V = det {6i: + fGijd''<f)'d,<f)') . (5.26) 

The potential results from the interaction of the brane with background fields 
and has been explicitly calculated in some particular cases. For the rest of this 
work the potential is left unspecified and general. 
It is useful to define a "kinetic matrix" as 

X^-^ = -]^d^(p^d^'(t)-\ (5.27) 

Noting that 

HpfM'l — pA/32/33/34 /lai 4^2 403 /ia4 (^OR^ 
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where 

= 5| + fBfBj,, Bf = Gijd'^cP', Bl = dpcf>', (5.29) 
and after using the identity 

ea,a.a3«4^"^-"^-^^"^-^* = "(4 - J^l^J^l " " " , (5.30) 

one can show that 

V = l- 2fGijX^-^ + 4fxfxf - SfxfxjX^^ + WfxfxjX^X^K (5.31) 

Because of the antisymmetrization over the field indices the last term vanishes if 
there are at most three fields. For two fields the last two terms of the previous 
expression vanish. For the single field case only the first two terms remain and one 
recovers the result obtained in single field DBI-infiation presented in the previous 
section. It turns out that the non-linear terms in X^"^ of equation (I5.3ip also 
vanish for a multi-field DEI model if one considers homogeneous configurations, 
i.e., configurations in which the fields only depend on time and X^"^ is given by 
X^-^ = l/2<j)^(j)'^ . This fact mislead the authors of |114l I115[ 1118] to ignore these 
terms even in the inhomogeneous case. References |116l I117j pointed out the 
importance of these terms and correctly conclude that they cannot be ignored. 
Because of these terms multi-field DBI-inflation is not a particular case of multi- 
field K-infiation as previously thought. In chapter [TJ we will study linear and 
non-linear perturbations in general multiple field inflationary model. This general 
class of models includes K-infiation and DBI-infiation as particular examples. We 
shall also discuss the non-Gaussianity produced in these models, focusing in the 
three point functions of the adiabatic and entropy perturbations. 



Chapter 6 

Non-Gaussianity from the 
trispectrum in general single field 
inflation 



6.1 Introduction 

The theory of slow-roll inflation generically predicts that the observed cosmic 
microwave background radiation (CMBR) anisotropies are nearly scale invariant 
and very Gaussian. Indeed, the latest observations of CMBR by WMAP5 [107] 
confirm these expectations. This constitutes one of the biggest achievements of 
modern cosmology. 

Despite its successes the theory of inflation still has many open questions. 
For example, we do not know the origin of the scalar field whose energy drives 
inflation, not to mention that we have never detected directly in the laboratory 
these kind of particles. The energy scale at which inflation happened is unknown 
by many orders of magnitude. There are many models of inflation that give 
similar predictions for the power spectrum of primordial perturbations, so which 
one (if any) is the correct one? 

For us to move a step forward in our understanding of the very early universe 
we have to work on two fronts. First the observational side. In the next few years. 



with improved experiments like the Planck satellite l], we will measure the CMBR 
anisotropies to an incredible resolution. For example the observational bounds on 
the bispectrum (the three point correlation function of the primordial curvature 



^http://www.rssd.esa.int/indcx.php?project=Planck 
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perturbation Q will shrink from the present WMAP5 value —9 < Jnl < HI 
|107j to I/nlI ^ 5 |iT9j . where the parameter f^i parameterizes the size of the 
bispectrum. It is because this parameter is constrained to be small that we say 
that the CMBR anisotropies are very Gaussian. The observational bounds on the 
trispectrum (four point function) will also tighten significantly from the rather 
weak present constraint of |rjv£,| < 10^ |12Ut 1121] to the future constraint of 
\tnl\ ~ 560 |122] . where tj^l denotes the size of the trispectrum. These previous 
observational bounds on the non-linearity parameters are for non-Gaussianity of 
the local type. These bounds change depending on the shape of the wave vectors' 
configuration |123] . This is one of the reasons why it is important to calculate 
the shape dependence of the non-Gaussianity. 

In face of these expected observational advances, it is then imperative to push 
forward our theoretical knowledge of our theories and calculate more observa- 
tional consequences of the different inflationary models to make a comparison 
with observations possible. One possible direction to be followed by us and many 
others is to calculate higher order statistics (like the trispectrum) of the primor- 
dial curvature perturbation. These higher order statistics contain much more 
information about the inflationary dynamics and if we observe them they will 
strongly constrain our models. Because these higher order statistics have a non- 
trivial momentum dependence (shape) they will help to discriminate between 
models that have a similar power spectrum (two point function). 

Calculations of the bispectrum for a single field inflationary model were done 
by Maldacena [124] . He showed that the primordial bispectrum is too small (of the 
order of the slow-roll parameters) to be observed even with Planck. Subsequent 
work generalized Maldacena' s result to include more fields and more complicated 
kinetic terms |125lll26lll27j . In |127j . Chen et al. have calculated the bispectrum 
for a quite general model of single field inflation. They showed that for some 
models even the next-to-leading order corrections in the slow-roll expansion may 
be observed. 

In this chapter, we will focus our attention on the calculation of the trispec- 
trum. In [128] , Seery et al. have calculated the trispectrum for slow-roll multi- 
field models (with standard kinetic terms) and they showed that at horizon cross- 
ing it is too small to be observed. But there are models of single field inflation, 
well motivated from more fundamental theories, that can produce a significant 
amount of non-Gaussianity, such as Dirac-Born-Infeld (DBI) inflation [IIH 1129] . 
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In |129j the authors have computed the trispectrum for a model where the infla- 
ton's Lagrangian is a general function of the field's kinetic energy and the field's 
value, their result was obtained using a simple method |130j . that only gives the 
correct leading order answer. In this chapter we will provide the equations nec- 
essary to calculate the next-to-leading order corrections to the trispectrum. We 
argue that for some models these corrections might become equally observable 
in the future. In fact, we will calculate the fourth order action in the uniform 
curvature gauge that is exact in the slow-roll expansion and therefore in prin- 
ciple one could calculate all slow-roll corrections to the trispectrum of the field 
perturbations. 

We will also compute the exact fourth order action for the curvature pertur- 
bation ( in the comoving gauge. For a simpler inflation model (with the standard 
kinetic term) this was recently done in [131] . However |131j did not consider sec- 
ond order tensor perturbations. We will argue that this is an oversimplification 
that leads to erroneous results. The reason simply being that at second order 
in perturbation theory, scalar degrees of freedom will source second order ten- 
sor perturbations and this will give a non-zero contribution for the fourth order 
action and hence the trispectrum. 

There are other reasons why we will perform the calculation in the comoving 
gauge. First of all, in doing so we work all the time with the gauge invariant 
variable C that is directly related with the observational quantities. The comov- 
ing gauge action can also be used for other practical purposes. For example, it 
can be used to calculate loop effects that can possibly have important observa- 
tional consequences. It can also be used to calculate the trispectrum of models 
where the potential has a "feature" (see [1321 1133] for an example of such cal- 
culation for the bispectrum). In the vicinity of the sudden potential "jump" the 
slow-roll approximation temporarily fails and one might get an enhancement of 
the trispectrum. There are well motivated models of brane inflation where the 
throat's warp factor suddenly jumps |134j . 

This chapter is organized as follows. In the next section, we introduce the 
model under consideration. In section [63] we shall study non-linear perturbations. 
First, we compute the fourth order action in the comoving gauge including both 
scalar and second order tensor degrees of freedom. After that we compute the 
fourth order action in the uniform curvature gauge. In section 16.41 we present 
the formalism needed to calculate the trispectrum. In section 16. 5[ we calculate 



CHAPTER 6. THE TRISPECTRUM IN SINGLE FIELD INFLATION 



97 



the leading order trispectrum using the comoving gauge action. We comment 
on previous works and on the observabihty of next-to-leading order corrections. 
Section 16.61 is devoted to conclusions. 

6.2 The model 

In this work, we will consider a fairly general class of models described by the 
following action 

S = ll d^xV^ [Ml,R + 2P(X, 0)] , (6.1) 

where is the infiaton field, Mpi is the Planck mass that we will set to unity 
hereafter, R is the Ricci scalar and X = —■^g^'^d^(j)d,^(j) is the infiaton's kinetic 
energy, g^^ is the metric tensor. We label the infiaton's Lagrangian by P and we 
assume that it is a well behaved function of two variables, the infiaton field and 
its kinetic energy. 

This general field Lagrangian includes as particular cases the common slow- 
roll inflation model, DBI-infiation [1101 1111] and K-inflation [135] . 

We are interested in flat, homogeneous and isotropic Friedman- Robertson- 
Walker universes described by the line element 

ds'^ = -dt^ + a^{t)Sijdx'dx^, (6.2) 

where a{t) is the scale factor. The Friedman equation and the continuity equation 
read 

= p, (6.3) 

p = -3Hip + P), (6.4) 

where the Hubble rate is H = d/a, p is the energy of the infiaton and it is given 
by 

p = 2XPx - P, (6.5) 

where Px denotes the derivative of P with respect to X. 

It was shown in [136j that for this model the speed of propagation of scalar 
perturbations ("speed of sound") is Cg given by 
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We define the slow variation parameters, analogues of the slow-roll parameters, 
as: 

(6.8) 

s^^. (6.9) 

We should note that these slow variation parameters are more general than the 
usual slow-roll parameters and that the smallness of these parameters does not 
imply that the field is rolling slowly. We assume that the rate of change of the 
speed of sound is small (as described by s) but is otherwise free to change 
between zero and one. 

It is convenient to introduce the following parameters that describe the non- 
linear dependence of the Lagrangian on the kinetic energy: 

S = XP^x + 2X^P^xx = (6.10) 

A = X^P^x + ^X^Pxxx, (6.11) 

n = X^Pxxx + 'ix^P,xxxx- (6.12) 
5 

These parameters are related to the size of the bispectrum and trispectrum. The 
power spectrum of the primordial quantum fluctuation was first derived in [136] 
and reads 

. _ 1 _ I 

' ~ 367r2 p + P ~ S^~e' ^^'^^^ 
where it should be evaluated at the time of horizon crossing Cg^k = a^H^. The 
spectral index is 

-1 = — — ^ = -2e-r/-s. 6.14 
ct In K 

WMAP observations of the perturbations in the CMBR tell us that the previous 
power spectrum is almost scale invariant therefore implying that the three slow 
variation parameters should be small at horizon crossing, roughly of order 10~^. 
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6.3 Non-linear perturbations 

In this section, we will consider perturbations of the background ( 16.2p beyond 
linear order. There is a vast literature on second order perturbations that is 
important when one is interested in calculating three point correlation functions, 
see for example [1241 1125t 11271 [T^ . We are interested in non-Gaussianities that 
come from the trispectrum and so we need to use third order perturbation theory. 
For that we need to compute the fourth order in the perturbation action. In this 
section we will obtain the fourth order action in two different gauges. As a check 
on our calculations we will compute the leading order (in slow roll) trispectrum 
in both gauges. We will follow the pioneering approach developed by Maldacena 
|124j and used in several subsequent papers [1251 11261 11291 fT28] . 

For reasons that will become clear later it is convenient to use the ADM metric 
formalism [138] . The ADM line element reads 

ds"^ = -N^dt^ + hij [dx' + N'dt) [dx^ + N^dt) , (6.15) 

where is the lapse function, A^* is the shift vector and h^j is the 3D metric. 
The action (16. ip becomes 

S = ^J dtd^xVhN {^^'^R + 2P) + ^ J dtd^xVhN^^ {EijE'^ - E^) . (6.16) 
The tensor Eij is defined as 

E^J = l{k,-V.N,-V,N?|, (6.17) 

and it is related to the extrinsic curvature by Kij = N^^Eij. Vj is the covariant 
derivative with respect to hij and all contra-variant indices in this section are 
raised with hij unless stated otherwise. 

The Hamiltonian and momentum constraints are respectively 

^^^R + 2P -2tt'^N-^P^x - N-^ {E,jE'^ - E^) = 0, 

V,{N-'Ei) -V,{N-'E) = 7rN~'V.<pPx, (6.18) 

where vr is defined as 

vr = </) - N^Vj4>. (6.19) 
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We decompose the shift vector TV* into scalar and intrinsic vector parts as 

where diN^ = 0, here indices are raised with 5ij. 

Before we consider perturbations around our background let us count the 
number of degrees of freedom (dof ) that we have. There are five scalar functions, 
the field 0, A^, ip, deth and hij ~ didjH, where if is a scalar function and deth 
denotes the determinant of the 3D metric. Also, there are two vector modes 
A^* and hij ~ diXj, where X"' is an arbitrary vector. Both A^* and X'' satisfy 
a divergenceless condition and so carry four dof. Furthermore, we also have a 
transverse and traceless tensor mode '-fij that contains two additional dof. Because 
our theory is invariant under change of coordinates we can eliminate some of these 
dof. For instance, a spatial reparametrization like = x* + 9*e(x, t) + e*^j^(x, t), 
where e and e*^^^ are arbitrary and c^je^^) = 0, can be chosen so that it removes 
one scalar dof and one vector mode. A time reparametrization would eliminate 
another scalar dof. Constraints in the action will eliminate further two scalar dof 
and a vector mode. In the end we are left with one scalar, zero vector and one 
tensor modes that correspond to three physical propagating dof. 

In the next subsection we shall use two different gauges that correctly pa- 
rameterize these dof. Because physical observables are gauge invariant we know 
that both gauges have to give the same result for the trispectrum for instance. It 
seems then unnecessary to perform the calculation twice in different gauges. In 
practice, we will see that both gauges have advantages and disadvantages and one 
is more suitable for some applications than the other. Furthermore, it provides a 
good consistency check on the calculation. 

6.3.1 Non-linear perturbations in the comoving gauge 

In this subsection, we will compute the fourth order action for the general model 
(16. ip in the comoving gauge. In this gauge the scalar degree of freedom is the so- 
called curvature perturbation ( that is also gauge invariant. There are a few works 
on this subject using this gauge, see e.g. |131] . where the authors have calculated 
the fourth order action for a standard kinetic term inflation but they neglected 
second order tensor perturbations. We will show that this is an oversimplification 
that may lead to an erroneous result for the four point correlation function. 



(6.20) 
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In the comoving gauge, the inflaton fluctuations vanish and the 3D metric is 
perturbed as 

50 = 0, 

hij = a'^e^'^hij, hy = 6ij + + ^7ifc7| H (6-21) 

where deth = 1, 7jj is a tensor perturbation that we assume to be a second order 
quantity, i. e. 7jj = 0(C^)- It obeys the traceless and transverse conditions 7- = 
d^'-fij = (indices are raised with 6ij). ( is the gauge invariant scalar perturbation. 
In f l6.2ip . we have ignored the first order tensor perturbations ^■^^'l^jQ^r■ This 
is because any correlation function involving this tensor mode will be smaller 
than a correlation function involving only scalars, see results of |124] . In the 
literature the second order tensor perturbations are often neglected, however they 
should be taken into account. The reason for this is because at second order the 
scalars will source the tensor perturbations equation. Later in this section, we 
will elaborate further on this point. Higher order tensor perturbations, like ^^^'Jij, 
do not contribute to the fourth order action. 

We expand N and A^* in power of the perturbation ( 

iV = 1 + ai + aa + ■ ■ ■ , (6.22) 
Ni = iv/'^ + iv/'^ + ■ ■ ■ , (6.23) 
^ = i^i + ^2 + --- , (6.24) 

where a„, iv/ ^ and ipn are of order 

Some useful expressions for the quantities appearing in fl6.18p . valid to all 
orders in perturbations but for 'jij = 0: 



-2a- V^^ {d,CdX + 2d,dX), 



(6.25) 
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E,jE'^ - E^ 



+a-'e-'< 



d^ipd^ip - 2diNj ( d^CN' + dXN^ ) - ANkdiCd'd"^ 



1 



3« ak„ 



-2diNj (a^C^V + dXd^i^) - Adidj^d^Cd'il) 
+2&C (djCNiN' + 2djCNid'i/j + djCdiipd'ip 



(6.26) 



2[H + (:)diN 



+2djN {d^CNi + diCN^) + 2diNd'^ip 



+N- 



(6.27) 



In the previous equations, indices in the right-hand side are raised with 5ij while 
indices in the left-hand side are raised with hij. In the rest of this section indices 
will be raised with 5 



Now, the strategy is to solve the constraint equations for the lapse function 
and shift vector in terms of C, and then plug in the solutions in the expanded 
action up to fourth order. 

At first order in a particular solution for equations (I6.18p is jl24[ 1125] : 



C 



(6.28) 



Oil = T7' = 0' ^1 = -77 + X, d x = a -^C- 

H H cj 

At second order, the constraint equation for the lapse function gives 
4H 

—d^ijj2 = -2a~^diC {dX + 2Hd'iJi) - 4ai {a^'^didX - 2SC) - 2al (S + 6A) 



(6.29) 
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and the equation for the shift vector gives 

2Hd^a2 - ^a-^a^ivf ^ = -a'^ {dkaid'diiji - d^aid^^Ji + d^d^^jji + did^dk^pi) . 

(6.30) 

Due to the fact that A^* is divergenceless and that any vector can be separated 
into a incompressible and irrotational part one can separate the contributions 

_ (2) 

from a2 and Ni in the previous equation. The irrotational part of Eq. fl6.30p 
gives 

2Ha2 = d^^d'F,, (6.31) 
and the incompressible part gives 



= -d-'F, + d-%d'Fu 



(6.32) 



where Fi is defined as the right-hand side of equation fl6.30p . The operator is 
defined by d~^{d'^Lp) = ip and in Fourier space it just brings in a factor of —1/k^. 

It was shown in [127] that to compute the effective action of order C", within 
the ADM formalism, one only needs to use the solution for the Lagrange multi- 
pliers and A^* up to order C"~^- Therefore in order to calculate the fourth order 
effective action the knowledge of the Lagrange multipliers up to second order is 
required. It is given in equations fl6.28p . (16.310 and (16.320 . 

The second order action is 



^2 



dtd X 



a'l-C,'-ae{dCf 



(6.33) 



The third order action is [1211 [HSl fT27j 



dtd X 



-eaC [dCf - (S + 2A) ^ + ^CC' 



m c 



.(6.34) 



CHAPTER 6. THE TRISPECTRUM IN SINGLE FIELD INFLATION 104 



The scalar fourth order action is 
1 



Si = ^ I dtd^xa^ 



a^hC {dCf + ( 2S + 9A + yH ) - 6Ca? (S + 2A) 



(6.35) 



Here, no slow-roll approximation has been made. The previous action has to be 
supplemented with the action containing terms with one and two tensors 



5*^2 



(6.36) 



5'^^2 = / dtd^x 



Canonical variable for quantization (n 



(6.37) 



In order to calculate the quantum four point correlation function we follow the 
standard procedure in quantum field theory. However there is an important 
subtlety here. The gauge invariant quantity ( is not the correct variable to be 
quantized, because it is not a canonical field . The canonical field to be quantized 
is the field perturbation 6(f), or a convenient parameterization (n defined by 



00 



(6.38) 



where 0o is the background value of the field. We will see that ( is related to (n 
by a non-linear transformation so for the power spectrum calculation both proce- 
dures of quantizing ( or (n give the same answer because the difference between 
these two variables is a second order quantity. However, for the calculation of 
higher order correlation functions (like the bispectrum or trispectrum) (n is the 
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correct variable to be quantized as it is linear in d(f). For the trispectrum, quantiz- 
ing ( gives different results from quantizing (n- We will find the relation between 
( and (n through the gauge transformation equations from the uniform curvature 
gauge (discussed in detail in the next subsection) to the comoving gauge. In the 
uniform curvature gauge the ansatz is 

(f){x,t) = 00 + S(t>{x,t), 

hij = a'^hij, hij = Sij + + -^ikjj H (6.39) 

where deth = 1 and is a tensor perturbation that we assume to be a second 
order quantity, i.e., ^ij = (9(50^). The gauge transformations are 

C(X) = Cn(x)+F2(Cn(x))+F3(Cn(x)), (6.40) 

Itj = lij{t)+fiij, (6.41) 

where F2(Cn),/^ij = Fs^^n) = C'(Cn) the terms coming from the second 

and third order gauge transformations respectively, they can be found explicitly 
in appendix O or in |124l 1131] . 

We now need to find the fourth order action for the variable Cn- Schematically, 
we write the different contributions as 

5C„ = ^4(Cn) + ^,^(Cn) + V(Cn) + ^3(i^2(Cn)) + ^2(i^2(Cn)). (6.42) 

The first three terms come from Eqs. fl6.35H6.3Tj) when we substitute ( with 
(n and '-fij with Eq. (16.411) . Due to the non-linear relation between ( and (n, 
the third order action for Eq. (16.341) . will after the change of variables give a 
contribution to the fourth order action like Ss{F2{(n)) ■ Similarly, the second order 
action, Eq. (16.331) . will also contribute with S'2(-F2(Cn))- In principle, one would 
also need to compute the third order gauge transformation as the second order 
action gives origin to fourth order terms like Cn-^3(Cn) (where F^l^n) is the third 
order piece of the gauge transformation). Fortunately, terms involving F^ can be 
shown to be proportional to the first order equations of motion for (n, therefore 
when computing the trispectrum these terms will vanish and we do not need to 
calculate the third order gauge transformation explicitly at this point. It can be 
easily seen that equations (I6.34p . (I6.35p . (I6.36p . (I6.37P or their counterparts in 
terms of C„, Eq. ( 16. 42^ . have terms that are not slow roll suppressed. However, 
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because in pure de Sitter space ^ is a gauge mode we expect the action (16.421) to 
be slow-roll suppressed (of order e). One can perform many integrations by parts 
to show that the unsuppressed terms of (I6.42p can be reduced to total derivative 
terms and slow-roll suppressed terms given by 



For us to be able to obtain the previous result it is crucial to include the con- 
tributions from the tensor actions (I6.36P and (I6.37p . otherwise the trispectrum 
calculated using the C(e°) of (I6.35P does not vanish, giving the wrong leading 
order result. Neglecting tensor perturbations (sourced by the scalars) for the 
calculation of the trispectrum is not consistent and leads to wrong results. This 
is one of the results of our work. The contribution from that comes through 
'-fij in (I6.36P and (16.371) will result in terms that are already slow-roll suppressed 
and no further integrations by parts are need on these terms (see Eq. (I6.59P of 
next subsection). The final action for is then 



where S*^^ denotes the terms of (I6.42p that are suppressed by at least one slow-roll 
parameter. This final action is slow-roll suppressed as expected and no slow-roll 
approximation was made so it is also exact. 

6.3.2 Non-linear perturbations in the uniform curvature 
gauge 

In order to calculate the intrinsic four point correlation function of the field 
perturbation we need to compute the action of fourth order in the perturbations. 
In this subsection we will obtain the fourth order action in the uniform curvature 
gauge. In this gauge, the scalar degree of freedom is the inflaton field perturbation 




(6.44) 
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54>{x^). There are several works in the hterature where the authors also calculate 
the trispectrum. In |128] . Seery et al. calculate the trispectrum of a multi- 
field infiation model, however the result is only valid for fields with standard 
kinetic energies, i.e. P{Xi, . . . , 0i, . . . , (pn) = Xi + ■ ■ ■ + Xn — V, where 
Xn is the kinetic energy of (pn and V is the potential. In this chapter we will 
generalize their result for an arbitrary function P{X, (p) but for single field only. 
Recently, Huang and Shiu have obtained the fourth order action for the model 
under consideration (16.11) . However the result was obtained by only perturbing 
the field Lagrangian. This procedure gives the right result as long as we are 
interested in the leading order contribution in the small speed of sound limit and 
in the slow-roll approximation, which was their case of interest. In the present 
section we will compute the fourth order action that is valid to all orders in slow 
roll and in the sound of speed expansion. To do that we have to perturb the full 
action ( 16. ip up to fourth order in the field perturbations. The procedure to obtain 
the fourth order action in this gauge is similar to the one used in subsection 16.3. II 
In this gauge, the infiaton perturbation does not vanish and the 3D metric 
takes the form 

= 00 + Hix^t), 

hij = a^hij, hij = Sij + + -^ik^j H (6.45) 

where deth = 1 and is a tensor perturbation that we assume to be a second 
order quantity, i.e., = O(50^). It obeys the traceless and transverse conditions 
7* = 9*7ij = (indices are raised with In the literature the second order 
tensor perturbations are often neglected, however based on our results it should 
be taken into account. We can always use the gauge freedom at second order 
to eliminate the trace and the vector perturbations of hij. The presence of 7jj 
makes the three dimensional hypersurfaces non-fiat so using the name uniform 
curvature gauge might be misleading. 

We expand and A^* in powers of the perturbation 6(f){x, t) 

N = l+ai + a2 + --- , (6.46) 
= iv/'^ + iv/'^ + ■ ■ ■ , (6.47) 
^ = V'l + ^2 + ■ ■ ■ , (6.48) 
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where Ni and ipn are of order 50" and (j)o{t) is the background value of the 
field. At first order in 6(f), a particular solution for equations fl6.18l) is |124l I128j : 



«i = — 0o50Px, iV>^=0, 9^^i = _- -^50 . (6.49) 
2-H cj at \ (j) J 

At second order, the constraint equation for the lapse function gives 
AH 1 

^•9^2 = ^ (S^tAi^^^i - + QH^ (3a? - 2^2) 



f^^^V'i - (300 al - 200 ^2 + 50 - 40o50ai j 



8i/ai 
a 

■\ ( 9,509*50 - 20o9,509Vi ) Px 



+50 P<^0 + 20o50 (^0oai - 50j Px<^ 
2 

of ' ' ' ' ^ 

lOa^0o50ai — 20o9i509Vi + 6a^0o a? 



a2 



-200 a'a2 + 4a'50 - di5<pd'5<i) jP^xx 
-00 (00 + 50 - 20o50ai j P^xxx 

+200^50 (^00"! - 50^ P,xx^ - 0o^50^Px00, (6.50) 
and the equation for the shift vector gives 

2Hd,a2 - ^a-^a^Tvf ^ = AaiHdiai + 0"^ (S.ai^^T/;! - dkaidid''i)i) 

+di6(f) (^50 - 0oai) P,x + 5^5000^ (50 - 0o«i) P,xx 
+diS^4>oScl)P^x<i„ (6.51) 

where P^^ means derivative of P with respect to 0. 

Due to the fact that iV* is divergenceless and that any vector can be separated 
into a incompressible and irrotational part one can separate the contributions 

~ (2) 

from 0^2 and Ni in the previous equation. The irrotational part of Eq. (I6.5ip 
gives 

2ifa2 = d-^d'Fi, (6.52) 
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and the incompressible part gives 

1 .~r(2) 



2a2 



-d-'F, + d-'d.d'Fk, 



(6.53) 



where Fj is defined as the right-hand side of equation (16.511) . 

The scalar fourth order action, where no slow-roll approximation has been 
made, is 

Si = SA + Sb, (6.54) 

where 



Sa 



j dt(fx 


a5(f) 




2 



•^100 ) + a^^o^i + ai{d5(j)y 



P,x<i> - ^ 



6Qa^a2 (^2 — 2a^) + 30Qa^a^ 



-180^a^a?5'/ + 8^la^d'6(p (n^ + (9,^2 ) ( 50 



-40o9i50<9Vi (^(550) - ^roa af + 80oa "i'^^ - 3a^50 J + 120oa ai50 
+40;^a^af(50 - 4 ^gSi^^SViafc^^a'^V^i - 2^la^al{d6(j)y - 4ai 00500^(550) 



-a' 50 -a-'(950)' 



12 L 



60o9j509Vi ( 20oai50 — 50 — 0oai 



+3(950)2 /'20Qai(50 - (/^^a^") _ IG^pyaiScf)"^ + 2402a2a?50^ - 12<pla^alc 



+000^0^ + 3a^ 

+ ((950)2 ^ 20oa,50a>i 



50 - a-\d64>y 



P 



XXX 



1 — 50) — 40oa2Q;i(50 



2(950)' - 2a'(50 + 4a'ai0o50 + 2a'0^a^ + 40o9i509Vi 



a2(50 + 30Qa2a250 

P,XX(f> 

P,XiP(p 



60oaf50 - 40oai50 - 40oQ!:^50 + (p^ai + 50 



+ 24*°"' 

-^0|^50a^ f-50^ + 30oai50^ + <plal - 3<plal6^) P,xxx<t> 



P 



xxxx 



,xx<i)4> 



1 



)o50V ( ai0o - 50) P 



(6.55) 
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dtd X 



1 

'2a 



ala^6(f)P^^ + -a^a?(50^P<^0 + -a^ai6(f)^P^ 



{di6(f)d'ipiy - 2ai6(j)a'di6(f)d'ilJi + aia'{d6(f)) 



24 



44> 



+2a^di6(j) (^ivf ^ + diij2^ (s'cj) - ai</)o) + <pla^a2 (02 - ^af) 



+ 



4a 



A~ri2) ^ ^7 (2) 



+ 12a^H'^a2 (aa - 2al) + Aaid'^ipid'^ip2 



(6.56) 



The previous actions should be supplemented with the pure tensor terms and the 
tensor-scalar coupling terms: 



S = ^j dtd'x [a%,f' - adk^.jd'f^ 



S 



dtd X 



(6.57) 



(6.58) 



This constitutes the main result of this subsection. It is a good check for our 
calculation to see that the previous action f l6.54p reduces in some particular cases 
to previously known results present in the literature. 

For example, if we restrict our model to the standard inflation case, i.e., 
P{X,(j)) = X — V{(j)), where V{(f)) is the infiaton potential, then all the terms 
in the scalar action (16.551) vanish and the only contribution to the fourth order 
action comes from (16.561) . These terms exactly reproduce the result of Seery et al. 
|128j . their equation (36), restricted to single field. However, in the total fourth 
order action there are also the tensor contributions ( 16.57P and (16.58^ . In general, 
to proceed one has to calculate the equation of motion for the second order tensor 
perturbations 7^^ from Eqs. (I6.57p . (16.581) to get 



% + 2-% - d%, 
■' a 



2PM(t)dM + 4Px(l)odj6(j)d^tlj, 



TT 



(6.59) 



where TT means the transverse and traceless parts of the expression inside the 
parenthesis (see appendix [D] for details of how to extract the TT parts of a ten- 
sor) and then solve this equation to obtain as a function of 6(j). One can 
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immediately see that at second order the scalars will source the tensor perturba- 
tion equation as it was previously shown by others [1391 1140] . At this order in 
perturbation theory, equation (16.591) should also have a source term quadratic in 
the first order tensor perturbations, ^^^lijcw We neglect these terms because we 
expect that any correlation function where ^'^^JijQ-^^ enters, which is sourced by 
the first order tensor modes squared, must be smaller than a correlation function 
with only scalars, see Ref. |124j for an example. In Fourier space, the source 
term of (I6.59P is suppressed by k"^, where k is the wave number. Once we have 
the solution of in terms of 6(j) we can plug back the result in (I6.57P and (16.581) 
to get the total fourth order scalar action. 

6.4 The general formalism to calculate the 
trispectrum 

6.4.1 The trispectrum of (n 

Now we shall give the basic equations needed to calculate the trispectrum [124[ 
1141] . First we need to solve the second order equation of motion for (n (obtained 
from (I6.33P ). Defining new variables 

a\/2e 

Vk = zuk, z = , (6.60) 

Cs 

where the Fourier mode function Uk is given by 

Uk = J c/=^xCn(t,x)e-'''-", (6.61) 
the equation of motion for (n is 

< + clPvk - -Vk = 0, (6.62) 



where prime denotes derivative with respect to conformal time r. This is also 
known as the Mukhanov equation. The previous equation can be solved, at 
leading order in slow roll and if the rate of change of the sound speed is small 
[T27] - to give 

Uk = u{r, k) = ^^£= (1 + ikc,T) 6-'"'^^ (6.63) 
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We do not need to impose any constraints in the sound speed and it can be 
arbitrary. Only its rate of change is assumed to be small. The next-to-leading 
order corrections to the previous solutions are also known and can be found in 
|127j . In the general case, we would have to solve Eq. (]6.62p without assuming 
slow roll. This can be done numerically. 

In order to calculate the Cn correlators we follow the standard procedure in 
quantum field theory. is promoted to an operator that can be expanded in 
terms of creation and annihilation operator as 

C„(r, k) = M(r, k)a(k) + M*(r, -k)a\-k). (6.64) 

The standard commutation relation applies 

[a(ki),at(k2)] = (27r)3(5(3)(ki - ka). (6.65) 

The vacuum expectation value of the four point operator in the interaction picture 
(at first order) is [HI] 

{n\ C„(t, ki)Cn(t, k2)Cn(t, k3)Cn(t, k4)|l^) 

= -l [ dm[Cn{tM)Cn{tM)Cn{t,ks)Cn{t,U),Hj{i)]\0), (6.66) 
J to 

where to is some early time during inflation when the inflaton vacuum fluctuation 
is deep inside the horizon, t is some time after horizon exit. \Q) is the interacting 
vacuum which is different from the free theory vacuum |0). If one uses conformal 
time, it is a good approximation to perform the integration from — oo to because 
T ^ —{aH)~^. Hi denotes the interaction hamiltonian and it is given by Hj = 
7r(n — L, where vr is defined as vr = ^ and L is the lagrangian. In this work, we 
will only calculate the contribution for the four point function that comes from 
a part of the interaction hamiltonian determined by the fourth order lagrangian 
Hj = — L4, where L4 is the total lagrangian obtained from the action fl6.44p . We 
should point out that the other terms that we do not consider here in the fourth 
order interaction hamiltonian are indeed important to obtain the full leading 
order result (see equation (16.751) ) as was recently shown by [ 129j . 

Of course in the end we are interested in the four point correlation function 
of C and not of (n- At leading order in slow roll these two correlation functions 
are equal but they will differ at next-to-leading order. 
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6.4.2 The trispectrum of ( 

In this subsection we calculate the relation between the trispectrum of (, on large 
scales, and the trispectrum of Cn calculated using the formalism of the previous 
subsection. This relation also involves lower-order correlation functions of (n 
present in the literature. The variables ( and Cn are related up to third order by 

C(X) = Cn(x) + F2(Cn(x)) + F3(Cn(x)), (6.67) 

where -^2(0 = C^(Cn)' ^siCn) = C^iCn) terms coming from the second and 

third order gauge transformations respectively. F2 can be found in appendix [Cl 
it is 

where (3 is given in Eq. (]C.6p . In the large scale limit (super-horizon scales), we 
can ignore P as it contains gradient terms. F3 was calculated in [131] and reads 




3H^o H(j)o 3 3 J \2H(j)o ) H 2H^ 

+ /a(Cn) + /b(C,7.,), (6.69) 

where fa denotes terms that contain gradients (it can be found in [131j ). fb is 
the part of the third order gauge transformations that contains 7ij. The explicit 
form of fb is to the best of our knowledge still unknown. To find out the explicit 
dependence of these terms on one would have to solve the equations of motion 
for 7jj, equation fl6.59p . We do not do this in this work. We believe that these 
terms will vanish in the large scale limit and therefore do not contribute to our 
calculation. 

A field redefinition like ( = Cn+ciiCn\n^ +a2Cn\n\ where (^^''^'^''^^ denotes one 
of Cn, Cn, Cn, givcs after using Wick's theorem a relation between both trispectrum 
like 

(C(xOC(x2)C(x3)C(x4))c = {T} + {PB} + {PPP} + 0(Pf )^ (6.70) 



where 



{T} = (Cn(Xl)Cn(x2)Cn(x3)Cn(x4)), 



(6.71) 
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{PB} = 

+ (C(")(xOCn(x4))(enXl)Cn,(x2)Cn(x3)) + (C^''nXl)Cn(x2))(Ci'^nXl)Cn(x3) 

+ (CfnXl)Cn(x3))(eHxi)Cn(x2)Cn(x4)) + (eHxi)Cn(x4))(enXl)Cn(x2)a(^^ 

+3 perm. 



+ a2 



c,b d) + 3 perm. 



(6.72) 



{PPP} = a?[(e)(Xi)e)(x2))((Cf (Xl)Cn,(x3))(e)(x2)Cn,(x4)) + (x3^X4; 

+ (Ci'^nXl)enX2))((Cf (xi)Cn(x3))(e)(x2)Cn(x4)) + (X3 ^ X4 
+ (Cf (xi)enX2))((d"nXl)Cn(x3))(d^nX2)Cn(x4)) + (X3 — - X4 
+ (Cf (xi)Cf (x2))f(d'^nXl)Cn(x3))(eHx2)Cn(x4)) + (X3 ^ X4 

+5 perm 

{a c,b ^ d) + 5 perm 



+ al 



+ 2aia2[(C^'^)(xi)d^Hx2))((d'Hxi)Cn(x3))(C^'^Hx2)Cn(x4)) + (x3— X4) 
+ (Ci'^Hx0d'^)(x2))((C(^)(x0Cn(x3))(Cl^)(x2)Cn(x4)) + (X3 ^ X4 
+ (enXl)enX2))((enXl)Cn(x3))(e)(x2)Cn(x4)) + (X3 ^ X4 
+ (d''nXl)d'^Hx2))f(eHxi)Cn(x3))(d^nX2)Cn(x4)) + (X3 — - X4 



+5 perm. 



(6.73) 



where "perm" means the other permutations of the preceding terms and (9(P^)^ 
denotes terms that are suppressed by successive powers of the power spectrum. 
{a —>■ c,b ^ d) means terms equal to the immediately preceding terms with 
a, b replaced by c, d respectively. (X3 < — > X4) means a term obtained from the 
preceding term with X3 and X4 interchanged. 

If the field redefinition contains third order pieces like C = Cn + &lCn"^Cn'^Cn'\ 
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they contribute with additional terms as 

(C(Xi)C(x2)C(x3)C(x4))e = (Cn(Xl)Cn(x2)Cn(x3)C„(x4)) 

+ h [(e)(Xi)Cn(x2)) ((d^)(Xl)Cn(x3))(Cl^Hxi)Cn(x4)) + (6 — c)^ 

+ (enXl)Cn(x2))((e)(Xl)Cn(x3))(enXl)Cn(x4)) + (« — c) 

+ (d'=Hxi)Cn,(x2))((d'^nXl)Cn(x3))(d^Hxi)Cn(x4)) + (« ^ 6) 

+3 perm.l + C(Pf )^ (6.74) 

where for example (6 < — > c) means a terms obtained from the term immediately 
preceding with b and c interchanged. To the best of our knowledge the expec- 
tation values involving operators containing derivatives of (n have not yet been 
calculated in the literature. However, once the mode function equation ( 16.62^ is 
solved, one has all the ingredients needed to calculate these expectation values, 
including the interaction Hamiltonian. 



6.5 Calculation of the leading order trispectrum 

In this section, we will use the formalism of the previous section and the fourth 
order exact interaction Hamiltonian of subsection 16.3.11 to calculate the leading 
order trispectrum, under the assumption that the "slow-roll" parameters (16.7116.91) 
are always small until the end of inflation. 



6.5.1 The leading order trispectrum of (n 

To calculate the leading order trispectrum of Cn in slow roll, we need to evaluate 
Eq. (16.661) where Hj is read from the order e terms of the action (16.441) . The 
interaction Hamiltonian we get contains terms with 7jj. Fortunately it can be 
shown that to compute the leading order trispectrum we do not need to know 
the solution for and the knowledge of its equation of motion (Eq. (16.591) ) is 
sufficient. At this order we use the solution for the mode functions Eq. (I6.63p . 
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The integrals in Eq. (16.661) can then be performed analytically to give 

H^ 1 



(f}|a(ki)C„,(k2)Cn(k3)Cn(k4)|fi) = (27r)^5^(ki + + + 



X 



1 (ion + 3A) - 7^ fsA - ^ + Hh) 4-A, - l^Aa 



(6.75) 



where the momentum dependent functions Ai are defined as 
A\ - 



_ klklik^ ■ k4) / , 3(A:3 + ki) , I2k^ki 

— — : i H — -— — 



A 



\ K 

(ki-k2)(k3-k4 



H 77^5 — + perm. 



3 



K 

+perm., (6.76) 

and "perm." refers to the 24 permutations of the four momenta. Note that 
the quantities of Eq. (I6.75P are evaluated at the moment r* at which the total 
wave number K = Et, fc. exits the horizon, i.e., when Kc.^ = ajl.. This 
leading order result agrees with the result of Huang and Shiu |129j □ that did 
their calculation in the uniform curvature gauge and using a simpler method 
that it is only valid to calculate the leading order contribution for models with 
<^ 1. In fact, we can compare our uniform curvature gauge result (I6.54p with 
the result of Huang and Shiu [129] . We see that the last terms of the fourth, 
sixth and ninth lines of equation (16.551) are exactly the ones obtained by [129j . 
their equation (15), using the method of just expanding the field Lagrangian as in 
|130l I142j . For a model with a general field Lagrangian these terms are the ones 
that give the leading order contribution for the trispectrum, in the small sound 
speed limit, equation (I6.75p . The contribution coming from the tensor part will 
be of next-to-leading order in this case. 

For standard kinetic term inflation, H = A = and = 1 and Eq. (16.750 
vanishes, the leading order is then given by the next order in slow roll. 



■^The full leading order result for the four point function can be found in the revised version of 
|129j which takes into account all the contributions for the fourth order interaction hamiltonian. 
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6.5.2 The leading order trispectrum of ( 

It is well known that if the slow-roll conditions are satisfied until the end of 
inflation and we can ignore gradient terms then the gauge invariant curvature 
perturbation ( remains constant on super-horizon scales to all order in pertur- 
bation theory. In this subsection, we will see that this fact greatly simplifies the 
relation between the trispectrum of ( and (n- 

In the large scales limit, Eq. fl6.67p simplifies to give 

C = Cn + aC + ^Un + 0{C), (6.77) 



where a is defined as 



+ (6.78) 
2 2H^o 



Using the fact that C = on super-horizon scales and the equation resulting from 
a time derivative of Eq. ( I6.77P one can show that 

C^ = -dC + 0{0. (6.79) 

This equation has a simple interpretation. The variable Cn is not constant outside 
the horizon, only the gauge invariant quantity ( is. This is the reason why the 
term jjCnCn in the second order gauge transformation cannot be ignored when 
one is calculating the trispectrum of (. Substituting Eq. (16.791) in Eq. (16.671) and 
taking the large scale limit we get 

C = Cn + aC^ + feC' + ■ ■ ■ (6.80) 

where ■ ■ ■ means cubic terms that contain at least one time derivative of and 
that will only give a contribution to the five point function. The variable b is 
defined as 

6n echo eri a 



3i/20o H(t)o 3 3 H 



6i720g 2H^o '^H^(PI 3 6 



2 ^2 

(6.81) 



-0 , ^ 



We shall now compare fl6.80l) with the result given by the 6N formalism [143^ 
I144[ 11451 1146] . In the SN approach ( is expanded in series in terms of the field 
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perturbation as 



( = N'S(j) + lN"6(j)^ + lN"'S(j)'^ + 0(5(j))\ (6.82) 
2 6 



where is the number of e-folds and a prime denotes derivative with respect 
to (j). Now comparing Eq. (16.801) with the previous equation and observing that 
Cn = —-^S(f) we expect 

2 02 6 03 

We verified that this is indeed the case. 

Using Wick's theorem one can now relate the connected part of the four point 
correlation function of ( with the four point correlation function of Cn calculated in 
the previous section |147l 1148] . This relation also involves lower order correlation 
functions of Cn, like the bispectrum (C„(xi)C„(x2)Cn(x3)) (the leading and next- 
to-leading order in slow roll bispectrum was previously calculated in |127j ). The 
relation is 

(C(X1)C(X2)C(X3)C(X4))C = (Cn(Xl)Cn(x2)Cn(x3)Cn(x4)) 
+2a (Cn,(Xi)Cn,(x2))(Cn(Xl)Cn(x3)Cn(x4)) + (Cn,(Xl)Cn,(x3)) (Cn(Xl)Cn(x2)Cn(x4)) 

+ (Cn(Xl)Cn(x4))(Cn,(Xl)Cn(x2)Cn(x3)) + 3 perm 

( Cn (Xl ) Cn, (X2 ) ) ( Cn (Xl ) Cn (X3 ) ) (Cn (X2 ) Cn (X4 ) ) 

+ (Cn(Xl)Cn(x2))(Cn(Xl)Cn(x4))(Cn(x2)Cn(x3)) + 5 perm 

(Cn(xi)C„(x2))(Cn(xi)C„(x3))(C„(xi)C„(x4)) + 3 perm] +0(Pf)^ (6.84) 



+4a 



+66 

where "perm" means the other permutations of the preceding terms and (9(P^)'^ 
denotes terms that are suppressed by successive powers of the power spectrum. 
Now, one can easily see that at leading order in slow roll the trispectrum for C 
and Cn are equal, this is because the constants a and b are slow-roll suppressed. 
These terms will only contribute to the next-to-leading order corrections. 
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6.5.3 The next-to-leading order corrections for the trispec- 
trum 

In subsection l6.5.1l we showed that for standard kinetic term inflation the leading 
order result, equation (16.751) . vanishes and in fact in this case the leading order 
of the trispectrum of C„ is of order (the next-to- leading order is the leading 
order). To obtain these leading order contributions it is easier to perform the 
calculation using the uniform curvature gauge action Eq. (I6.55I) - (I6.58I) . Eq. 
(16.551) vanishes exactly for standard kinetic term inflation. The action (16.561) is 
exact in the slow-roll approximation but it is instructive to determine the slow-roll 
order of the different terms. One can see that the leading order contribution comes 
from terms of order 0{e^), as pointed out in |128] . If we take P{X, 0) = X — V{(j)) 
then the leading order (in slow roll) source of Eq. (I6.59P will be of order 0{e^). 
We therefore do not expect to be slow-roll suppressed and the actions (I6.57p . 
(16.581) will contain unsuppressed terms of the same order as the leading order 
term of the action (I6.56p . These tensor contributions were absent in the analysis 
of |128] and we have shown that they are of the same order as the fourth order 
action considered in |128j . our Eq. (16.561) . It is still an open question how these 
new contributions will change the trispectrum result of Seery et ai. 

For the general Lagrangian case, the leading order trispectrum was given in 
the previous subsection and in [129] . Contrary to the method of [129] , our method 
of obtaining the fourth order action (16.541) does not rely on any approximation 
and therefore the action (16.540 is valid to all orders in slow roll and in the sound 
speed expansion and it can be used to study the next-to-leading order corrections. 
Depending on the momentum shape of these next-to-leading terms they might 
become big enough to be observed in the next generation of experiments. A simi- 
lar argument applies for the next-to-leading order corrections for the bispectrum, 
as it was shown in [127] . For example, for DBI inflation, |129] showed that the 
leading order non-Gaussianity parameters t^l scales like tnl ~ 0.1/Cg (for a spe- 
cific momentum configuration) and /tvl ~ V^s- They argue that if ~ 0.1 then 
fisiL is still inside the value range allowed by observations but tnl ~ 10^ could 
be detected with the Planck satellite CMBR experiment. Therefore, assuming 
that the slow-roll parameter e is of order e ~ 0.01 (at horizon crossing) these 
next-to-leading order corrections for the trispectrum could possible be observed 
with the Planck satellite. A more careful and systematic study of the momentum 
dependence of these new terms is required and it is left for future work. 
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6.6 Conclusion 

We have computed the fourth order action for scalar and second order tensor 
perturbations in the comoving gauge. Our resuh is exact in the slow-roll (SR) 
expansion but practically it is useful to study the SR suppression of the different 
terms. We were able to show that after many integrations by parts the unsup- 
pressed terms contained in the previous action can be reduced to total derivatives 
terms plus corrections that are SR suppressed. The resulting action has the cor- 
rect order in SR. It is suppressed by e as it should be, because in pure de Sitter 
space the curvature perturbation is a pure gauge mode. An important lesson from 
our work is that in order to obtain the correct SR order for the action, the sec- 
ond order tensor perturbations cannot be ignored as assumed in previous works 
|128j and [131j . We found the explicit form of these tensor perturbations in the 
comoving gauge by using the gauge transformations from the uniform curvature 
gauge. Fortunately, for a general inflation model like (16. ip . we showed that we do 
not need to solve the equations of motion for the tensor perturbations if we are 
interested in calculating the leading order trispectrum. However, to calculate the 
next-to- leading order corrections to that result, or to calculate the leading order 
trispectrum for standard kinetic term inflation, we do need to solve explicitly 
the equations of motion for the tensor perturbations. This will be left for future 
work. 

Using the comoving gauge action we have calculated the leading order in SR 
trispectrum of (. We compared our result with the result of [129] . obtained using 
the uniform curvature gauge, and we found an agreement. 

For the uniform curvature gauge action, that is also exact in the SR expansion, 
we identified the terms that will contribute to the next-to-leading order correc- 
tions to the trispectrum. We pointed out that depending on the model and on the 
momentum configuration, some of these corrections might be observable with the 
Planck satellite. After taking particular limits, the previous action nicely reduces 
to previously known results |128[ I129j with the caveat that the above mentioned 
works ignore tensor contributions. 

Finally we have obtained the relations between the trispectrum of ( and 6(f) 
(on large scales) using the third order gauge transformations and compared the 
result with the 6N formalism. 

To conclude, we have provided the necessary equations (fourth order action 
and the relation between ( and 6(f)) to calculate the trispectrum for a fairly 
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general model of inflation that are also valid for models where SR is temporarily 
interrupted, i.e., around a "step" in the inflaton's Lagrangian |134j . In this case, 
it is impossible to apply the 6N approach and it is required to evaluate the n- 
point functions numerically |132] (see |149j for a different approach). We leave 
this more practical application of our results for future work. 



Chapter 7 



Non-Gaussianity from the 
bispectrum in general multiple 
field inflation 



7.1 Introduction 

In this chapter, we shall study the non-Gaussianity of the primordial curvature 
perturbation produced by a phase of multiple field inflation. In particular, we 
shall concentrate on the bispectrum or three point correlation functioijl]. 

The simplest single field inflation models predict that the non-Gaussianity of 
the cosmic microwave background (CMB) fluctuations will be very difficult to be 



detected even in future experiments such as Plancj!|. The detection of large non- 
Gaussianity would mean that the simplest models of inflation are rejected. There 
are observational hints I152j of some deviation from perfect Gaussianity in 
the WMAP3 data |153] but the WMAP5 data is consistent with Gaussianity 

There are a few models where the primordial fluctuations generated during 
inflation have a large non-Gaussianity. In the single field case, if the inflaton 
field has a non-trivial kinetic term, it is known that the non-Gaussianity can be 
large. For example, in i^-inflation models where the kinetic term of the inflaton 
field is generic, the sound speed of the perturbations can be much smaller than 



^This chapter in based on our work |150j . While we were writing up that work, similar 
results appeared on the arXiv [117j . 

■^http:/ /www.rssd.esa.int/index.php?project=Planck 
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1 [135[ 1136] which leads to large non-Gaussian fluctuations. Additionally, Dirac- 
Born-Infeld (DBI) inflation, motivated by string theory, can also give large non- 
Gaussian fluctuations [1101 lllll 11121 1113] . The inflaton is identifled with the 
position of a moving D3 brane whose dynamics is described by the DBI action. 
Again, due to the non-trivial form of the kinetic term, the sound speed can 
be smaller than 1 and the non-Gaussianity becomes large [1271 I155j . In the 
previous chapter, the third and fourth order actions for a single inflaton fleld 
with a generic kinetic term have been calculated by properly taking into account 
metric perturbations. Three and four point functions have also been obtained 
|127l 11291 1156] . For the detailed observational consequences of single- fleld DBI- 
inflation see [1571 [M [TSU [1621 [T63lfT6ll . 

Multi-fleld inflation models where the curvature perturbation is modifled on 
large scales due to the entropy perturbations have been also recently studied. 
In the case of the standard kinetic term, it is not easy to generate large non- 
Gaussianity from multi-fleld dynamics |165l I166[ 11671 1168[ 11691 1170] (see how- 
ever Ref. |171j and the curvaton scenario e.g. [1721 11731 1174] ). In the DBI- 
inflation case, the position of the brane in each compact direction is described 
by a scalar fleld. Then DBI- inflation is naturally a multi-fleld model |114j . The 
effect of the entropy perturbations in the inflationary models based on string the- 
ory constructions in a slightly different context is also considered in |175l 1176] . 
Recently, Huang et al. calculated the bispectrum of the perturbations in multi- 
fleld DBI-inflation with the assumption that the kinetic term depends only on 
X = —G^'-'d^(j)^d^(ly^/2 where 0"^ are the scalar flelds (/ = 1,2,...) and Gjj is 
the metric in the fleld space, as occurs in K-inflation [115] . They found that in 
addition to the usual bispectrum of adiabatic perturbations, there exists a new 
contribution coming from the entropy perturbations. Then they showed that the 
entropy fleld perturbations propagate with the speed of light and the contribu- 
tion from the entropy perturbations is suppressed. This property can also be 
conflrmed by the analysis of a more general class of multi-field models where the 
kinetic terms are given by arbitrary functions of X |118l 1177] . However, Langlois 
et al. pointed out that their assumption cannot be justified for the multi-field 
DBI-infiation [116] . Even though the action depends only on X in the homo- 
geneous background, there exist other kind of terms which contribute only to 
inhomogeneous perturbations. They find that this dramatically changes the be- 
haviour of the entropy perturbations. In fact, it was shown that the entropy 
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perturbations propagate with the same sound speed as the adiabatic perturba- 
tions. Ref. [116] also calculated the bispectrum by generalizing the work of Huang 
et al. [129] . 

In this chapter, we study a fairly general class of multi-field inflation models 
with a general kinetic term which includes K-infiation and DBI-infiation. We 
study the sound speeds of the adiabatic perturbations and entropy perturbations 
and clarify the difference between K-infiation and DBI-inflation. Then we calcu- 
late the third order action by properly taking into account the effect of gravity. 
We continue to obtain the three point functions at leading order in slow-roll and 
in the small sound speed limit. We can recover the results for K-infiation and 
DBI-infiation easily from this general result. 

The structure of the chapter is as follows. In section 17.21 we describe our 
model and derive equations in the background. In section 17.31 we study the 
perturbations using the ADM formalism. Additionally, the second and third order 
actions are derived by properly taking into account the metric perturbations. 
Then we decompose the perturbations into adiabatic and entropy directions and 
write down the action in terms of the decomposed fields. In section 17.41 we study 
the sound speed in several models including K-infiation and DBI-infiation. It 
is shown that in general, adiabatic and entropy sound speeds are different and 
both can be smaller than 1. In section 17.51 the third order action at leading 
order in slow-roll and in the small sound speed limit is obtained in terms of the 
decomposed fields. Then the three point functions are derived for a generalized 
model which includes K-infiation and DBI-inflation as particular cases. Section 
17.61 is devoted to the conclusion. 



7.2 The model 

We consider a very general class of models described by the following action 

S = \j d'x^g [Ml^R + 2P{X'', 0^] , (7.1) 

where cj)^ are the scalar fields (/ = 1,2, ...,A^), Mpi is the Planck mass that we 
will set to unity hereafter, R is the Ricci scalar and 

X" ^ (7.2) 
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is the kinetic term, (yf^j, is the metric tensor. We label the fields' Lagrangian 
by P and we assume that it is a well behaved function so that derivatives of P 
with respect to X^"^ can be defined. Greek indices run from to 3. Lower case 
Latin letters ...) denote spatial indices. Upper case Latin letters denote field 
indices. See section 16.21 for the single field case with non-standard kinetic term. 
The Einstein field equations in this model are 

G^. = Pg^u + P,x'^d^(P'd,<p' = T^,, (7.3) 

where P^x^ denotes the derivative of P with respect to X^"^ §. The generalized 
Klein- Gordon equation reads 

g^''{Pxud,^').^^ + Pj = 0, (7.4) 

where ; denotes covariant derivative with respect to g^^i, and Pj denotes the 
derivative of P with respect to 0*^. 

In the background, we are interested in fiat, homogeneous and isotropic 
Friedman- Robert son- Walker universes described by the line element 

ds"^ = -dt^ + a^{t)5ijdx'dx\ (7.5) 

where a{t) is the scale factor. The Friedman equation and the continuity equation 
read 

= Eo, (7.6) 

Eo = -3H{Eo + Po), (7.7) 

where the Hubble rate is H = d/a, Eq is the total energy of the fields and it is 
given by 

Eo = 2Xo"Po,x^/ - Po, (7.8) 

where a subscript zero denotes background quantities, Xq"' = 1/20Q0Q. The 
equations of motion for the scalar fields reduce to 

Po,x'^4>o + {^HPo^x^j + Po,xu) 4 - Po,j = 0. (7.9) 



^Strictly speaking, we adopt the symmetrized derivative, P^x'-' = 5 ( + aX^l 
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7.3 Perturbations 

In this section, we will consider perturbations of the background ( 17.5P beyond 
linear order. For this purpose, we will construct the action at second and third 
order in the perturbations and it is convenient to use the ADM metric formalism 
[m [TM [m [m im im [T77j . The ADM llne element reads 

ds"^ = -N^dt^ + hij [dx' + N'dt) [dx^ + iV^dt) , (7.10) 

where is the lapse function, A^* is the shift vector and hij is the 3D metric. 
The action (17. ip becomes 

S = ^J dtd^xVhN + 2P{X^\ 0^)) + ^ y dtd^xVhN^^ {EijE'^ - E^) . 

(7.11) 



The tensor Eij is defined as 



E^, = ^ (/i., - V^N, - V.AT,) , (7.12) 

and it is related to the extrinsic curvature by Kij = N~^Eij. Vj is the covariant 
derivative with respect to hij. X^-^ can be written as 

1 Ar-2 

X^J = --h'Jdi(f)^dj(f)^ + -^v^v^ , (7.13) 

where f ^ is defined as 

= <j)^ - N^Vjcj)^. (7.14) 
The Hamiltonian and momentum constraints are respectively 



+ 2P - 2N-^P^x'JV^v^ - {EijE'^ - E^) = 0, 

V, {N^^Ej) - Vi {N~^E) = N^^P^x'-'V^Vi 



/V7 



(7.15) 

We decompose the shift vector A^' into scalar and intrinsic vector parts as 

Ni = Ni + diilj, (7.16) 
where diN'^ = 0, here and in the rest of the section indices are raised with 6ij. 
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7.3.1 Perturbations in the uniform curvature gauge 

In the uniform curvature gauge, the 3D metric takes the form 



7 _ 2 r 



(j)'{x,t) = <j)',{t) + Q'{x,t), 



(7.17) 



where denotes the field perturbations. In the following, we will usually drop 
the subscript "0" on 0q and simply identify (p^ as the homogeneous background 
fields unless otherwise stated. 

We expand and A^' in powers of the perturbation 



N 



1 + «! + 02 + 



ijj = 1pl+1p2^ , 



(7.18) 
(7.19) 
(7.20) 



where a„,, Ni and tpn are of order (Q^)'^. At first order in , a particular 
solution for equations fl7.15p is: 



ai 



0, 



2H 

2„, , o r^K llljp 



a 

2H 



+ ( PxiJ + ^''^'''PxLMxi-^ ' ^'^^ 



The second order action is calculated as 



(7.21) 



(2) 



dtd X — 



X^'^X^^^Px'Jxi^M + Px'jQ^Q^ - a-'^P^x'jd'Q^diQ^ 
+2P^x'.jkQ'' + ¥Q') + P,klQ''Q'' 



where 



(7.22) 



(7.23) 
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After integrating by parts in the action and employing the background field equa- 
tions, the second order action can be finally written in the rather simple form 



1 



-^P^jjdiQ'd'Q' - MuQ'Q' + AfijQ'Q'] , (7.24) 



with the effective squared mass matrix 



M 



ij 



-P.IJ + 



K 



H 



{PxJkPjxlm + P^x'^P^jxi 



._X^'''xPQP^xmnxpqPxikP,x^l<P''^'^ 



1 d 



^Pj,rKPxj4''<P'' 



Ml 



IJ 



2 P 



X 



MN 



H 



-P 



XlK XMN P X-JL 



L \ iK 



(7.25) 
(7.26) 



In the same way, the third order action is given by 



S'(3) = / dtd^xa^ 



ai 



+ 



+ 



1 



P,xijkQ 



K 



+ -^KlJ^ V + 77-^1 -^1 -^1 ^,XiJXLMxQR + -y^i y^i i^,X'JXLMj^Q 
2d 2 



Ixi'P^x'^LMQ'^Q'' + IpjjkQ'Q'Q'' 

2 D 



(7.27) 



7.3.2 Decomposition into adiabatic and entropy pertur- 
bations 



We can decompose the perturbations into the instantaneous adiabatic and entropy 
perturbations, where the adiabatic direction corresponds to the direction of the 
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background fields' evolution while the entropy directions are orthogonal to this 
|178j . For this purpose, following |118j . we introduce an orthogonal basis e^(?T, = 
1,2,..., TV) in the field space. The orthonormal condition is defined as 

P,x'reiei = Snm, (7.28) 

so that the gradient term P xudiQ^ d'^Q'' is diagonalizecl^. Here we assumed that 
P^xi-^ is invertible and it can be used as a metric in field space. The adiabatic 
vector is 

, (7.29) 

which satisfies the normalization given by Eq. (17.28^ . The field perturbations are 
decomposed on this basis as 

Q' = Quel . (7.30) 
We defined the matrix Z^n which describes the time variation of the basis as 

~ ^m^rnn i (7-31) 

which satisfies = — ^nm — Px'-'^L^i as a consequence of (P^xiJ^n^m) = 0. 

In terms of the decomposed fields, the second order action fl7.24p can be 
rewritten as 

If 1 

S(2) = 2 dt(f^<^^[^mn{DtQm){DtQ„) - -^SmndiQmd'Qn 

—■MmnQmQn + J^mnQn{DtQm)~\ , (7.32) 

where 

DtQm ~ Qm ~\~ ZjYinQn ; (7.33) 

JCmn = 5mn + {P,xM^^"'^'')P,x^KxJLe{e'^eie';,, (7.34) 
Mmn = Mijeiei, (7.35) 
AC,„ = Uijelei (7.36) 



*It is worth noting that the normahzation given in [116|,lll7|lll8j is Gijel^e'^ = Snm, unless 
Px'-' = Gij, that is, the case with the canonical kinetic terms, our results written in terms of 
Qn are different from the results of |116[ 11171 1118] . But it can be shown that when written in 
terms of w„ which are the canonical variables for the quantum theory, our results are consistent 
with the results of [m [TTflfTTg) . 
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Prom the constructions, fCmn, M.mn and Nmn are symmetric with respect to m and 
n. The exphcit form of the effective squared mass matrix in this representation 
is 



Mr 



H 



4if2 
1 d 
dt 



'P 



aK1L\3 



XKL( 



— {P^XMN 



In M N P Q\i; s 



Px'KPx-j^e^e^ 



-AP,xPQ<p''<p'^)\PxKLxMNe^e^e^'e^)6ni 



H 



where P^n = P,ije'^ei and P^nX'^ = P^x'^e:^ 
The equation of motion is obtained as 

1 d 



I^O {^2}CYnr D tQ rn ~^ -^mrQ m)~\ i^^mn^ nr ~^ ■^rar-) DtQr 

2 



a3 dt 



^d^Qr = . 



(7.37) 



(7.38) 



(7.39) 



7.4 Linear perturbations 

In this section, we study the hnear order perturbations using the second order 
action derived in the previous section. 



7.4.1 K- inflation 

Let us consider K-inflation models where P{X^'^ , (p^) is a function of only the 
trace X = X^^ Gi,j{(f)^) of the kinetic terms where Gij{(f)^) is a metric in the 
field space: 

P(X",0O = n^,0')- (7.40) 
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The derivatives of P can be evaluated as 

P^x'-J = GjjP^x, (7.41) 

Pi = ^GjK,i<p'<p''Px + Pi, (7.42) 

P,xi-'xKL = GijGklP,xx , (7.43) 

P,x''K = Iglm,k<P''<P^'GuP,xx + Gij,kP,x + GuPxk, (7.44) 

Pij = \Gkl,iGmnA''^''^^'^''P,xx + ^Gxi,,j0^0^Px 

+^(/.*^0^(GMiv,jPx/ + GuN^iP^xj) + P/j , (7.45) 

and the sound speed is defined as 



cl ^ li^^ . (7.46) 

In terms of the decomposed field, the second order action can be written as 



5'(2) ~ 2 ^ dtd^xa? 



Smn + ( 4- - 1 ) ^Im^ln } {P>tQm){P>tQn) 



cl 



a? 



(7.47) 



where we do not show the explicit forms of M.mn and Mmn- 

The fact that the sound speeds for the entropy perturbations are unity has 
been recognized in Ref . [118] . This is because the non-trivial second derivative of 
P only affects the adiabatic perturbations. 

7.4.2 DBI-inflation 

An interesting class of models is the DBI-inflation which describes the motion of 
a D3 brane in a higher dimensional spacetime. The DBI action is given by [179j 

S = - j d'xj^^^-det[g,, + f{(^^)Gij{(^^)d,4>^dAJl (7.48) 

where det denotes a determinant, G^^ is a metric in fleld space and (j)^ corresponds 
to positions of a brane in higher dimensional spacetime. Recently it was pointed 
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out by Ref. [116] that multi-field DBI-infiation is not included in multi-field K- 
inflation discussed in the previous subsection. Indeed, P(X^'^) is not a function 
of X, but it is given by 

P(X",0O = ^(^,0'), ^=^^^' (7-49) 

where 

V = det{5'j-2fX'j) 

= 1 - 2fGijX^-^ + 4/2x/xjl - 8/=^x/x/xJl + IG/^x/x^/xf xfl 

(7.50) 

where [ ] denotes the anti-symmetric bracket and we used G^^ to raise the indices. 
In the background, X = X. However, this does not mean that the full action is 
a function of X only. The DBI action takes a specific form of P 



P{X, <!>') = -J (\/i-2/x - 1^ - 



(7.51) 



where we allow for a potential V{(f)^). The sound speed is defined as 

cl = . (7.52) 

The derivatives of P evaluated on the background can be calculated as 

I dx\ I dX \ ~ I d^X \ 

= ^-^ [lxn)[lx^^)^ («Twfr j • (^-5*' 

where 

dX 



■clGrj + 2fXjj, (7.55) 



dX'J 

d^X _ / „ „ 1 „ „ 1 



dX^-'dX^^ ^2/ ^GuGkl — -^GikGjl — -GilGjkJ + 0{X ) . (7.56) 
Here we do not explicitly write down the higher order terms in X^^ in the second 
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derivative as they will not contribute to the final result. In the following, we will 
omit these terms. We can also show that 



Pi = -GjK,i^'^''P^<i + Pi, (7.57) 
Pij = -,GKmGMN,A''4>''4>''^''P^xx + ^GkljA'^'-P^x 



1 

-( 

2 

\Gkl,iGmnA''^^''^''P^xx + \^ 
+ ^0^^0^(Gm^,jP^, + GmnjP^xj) + PiJ . (7-58) 

P,X'-'K = {{^ — '2fX)Gjj^K ~ '2f,KXGjj — 2fGLM,KX^'^^Gjj 

'2f,KXij + 2fGjL,KX^j + 2fGjM,KXj^^)Px 

+ {clGu + 2fXu) (^IglmM''<P''P,xx + Pxk^ ■ (7-59) 

It is worth noting that even though P^x^k seems to be a bit complicated, we can 
show that 

P,X^'K^' = \GLM,K^''^^'hP,XX + ^lP,XK + GlJ,K^'P,X , (7-60) 

which is just the same form as K-infiation case. We can also show that 

Pjf7j0^0^ = 2XP^ . (7.61) 

The orthonormality conditions for the basis give 

1 1 1 - 
e-nemi = — -K.n " ^; • (7.62) 

P,xC^s P,X 

Using these results, the second order action can be written in terms of the de- 
composed perturbations as 



S{2) = ^ / dtd^xa^ 



(7.63) 



Unlike K-infiation models, all field perturbations have the same sound speeds 
as was pointed out by Ref. |116j . In order to understand the difference between 
K-inflation and DBI-infiation, we will consider a generalized model where both 
cases are included. 
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7.4.3 Generalized case 

Let us consider models described by 

P{X'-\4>') = P{Y,(P'), (7.64) 

where ^ 

Y = Gjj{^^)X^-' + ^(X^ - X/Xj) . (7.65) 

The functional form of Y is chosen so that Y = X = GjjX^'^ in the background 
as in the DBI-infiation model. This model includes as particular cases K-inflation 
model for 6 = and the DBI-inflation for b = —2/ and if P has the DBI form. 
This might be surprising as the DBI action contains additional terms of order 
and p in X (see equations (17.491) and fl7.50p ). but it turns out that these terms 
do not contribute to the second order action and the leading order third order 
action. 

Following a similar procedure to the previous subsection, the second order 
action can be written in terms of the decomposed perturbations as 



S'(2) = ^ I dtd^xa? 



Smn H -p^ ^imSln + ^ ^ (^nl^ml "~ ^mn) 



a? 



(7.66) 



Now we are in a position to explain the difference between K-inflation and 
DBI-infiation. As in K-infiation case, the non-trivial second derivative of P affects 
only the adiabatic perturbations. On the other hand, the non-linear terms of X^^ 
in Y only affects the entropy perturbations. This can be seen from the fact that 
the sound speed for adiabatic perturbations c^^ and for entropy perturbations Cg^ 
are given by 

cl,^ ~ , cL = l + 6X, (7.67) 

and they are independently determined by Pyy and (PY / {dX^'^dX^^) respec- 
tively. Thus in general they are different. Let us derive the condition under which 
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the two sound speeds are the same, i.e., c^^ = c^^. This condition is given by 

2X^ = (7.68) 
Py l + bX ^ ^ 

Then we find that the DBI action is a solution of this equation where b = —2/ 0, 
although we should note that it has not been proved that multi-field DBI infiation 
is the only case where the sound speeds are the same. 



7.5 The leading order in slow-roll three point 
function 

In this section, we will calculate the leading order in slow-roll third order action 
for the generalized model of the previous subsection and then we shall calculate 
the leading order three point function for both adiabatic and entropy directions. 
Finally we will obtain the three point function of the comoving curvature pertur- 
bation. 



7.5.1 Approximations: slow-roll 

In order to control the calculations and to obtain analytical results we need to 
make use of some approximations. We will use the slow-roll approximation, where 
we define a set of parameters and assume that these parameters are always small 
until the end of infiation. We define the slow-roll parameters as 

' = -Jp = ^^ V=^^ (7-69) 

Xad = 77, Xen = 77- (".70) 

It is important to note that these slow-roll parameters are more general than the 
usual slow-roll parameters and that their smallness does not necessarily imply 
that the fields are rolling slowly. Assuming that the parameters Xad and Xen are 



^Wc thank Gianmassimo Tasinato for a discussion on this point. While we were writing up 
this work, there appeared a paper in the arXiv [117j . It is argued that the speciahty of the DBI 
action comes from the fact that it describes the fluctuations of the positions of the branc in the 
higher-dimensional spacetime and they propagate at the speed of light. From the point of view 
of observers on the brane, the sound speeds are smaller than 1 due to the Lorentz factor 1/c^. 
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small implies that the rates of change of the adiabatic and entropy sound speeds 
are small, but the sound speeds themselves can have any value and we assume 
that they are between zero and one. 

It is convenient to define a parameter that describes the non-linear dependence 
of the Lagrangian on the kinetic term as 

A = '^X'P,YYY + X^P^YY. (7.71) 

We will also assume that the rate of change of this new parameter is small, as 
given by 

(7.72) 

At the end of this section, we will show that the size of the leading order three 
point function of the fields is fully determined by five parameters evaluated at 
horizon crossing: e, X, H and both sound speeds. 

It turns out that the equations of motion for both adiabatic and entropy per- 
turbations at first order form a coupled system of second order linear differential 
equations, see appendix [E] for details. In general, the coupling (denoted by ^ in 
equation (lE.lOl) ) between adiabatic and entropy modes cannot be neglected but 
in this work we will study the simpler decoupled case, where we assume that ^ 
is small when the scales of interest cross outside the sound horizons, i.e., we will 
assume that ^ ~ 0{e). With these approximations the adiabatic and entropy 
modes are decoupled and the system of equations of motion can be solved ana- 
lytically. For simplicity, we will also assume that the mass term present in the 
entropy equation of motion is small, i.e., /i^/iiT^ <^ 1 (refer to appendix [El for 
more details). When calculating the leading order three point functions, we as- 
sume that the quantities related to the time derivatives of the basis vectors given 
by Zmn are also slow-roll suppressed. Finally, the calculation of the three point 
functions in the next subsections is valid in the limit of small sound speeds. Our 
results will also include sub-leading terms of 0{1) but these terms will in general 
(for small sound speeds) receive corrections coming from terms of the order of 
e/c^, that we have neglected. 
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7.5.2 Third order action at leading order 

At leading order in the previous approximations, we can neglect terms containing 
«!, ipi and derivatives of P with respect to the fields. Then the third order action 
for the general model (17.11) is calculated as 

3(3) = 11 dx^dta' PxiKxJL^^'Q'^^Q-'Q'' - ^Px^KxJL^^'Q'^^diQ'd'Q'^ 
2 J I ' 

+^Px-x-xMiv0(^Q^)0(^Q^)0(^^g^)l . (7.73) 

After decomposition into the new adiabatic/entropy basis the third order action 
can be written as 



5(3) = J dx^dta? 
where we define the coefficients 'E.nmi and Tnmi as 



„^nmlQnQmQl „ r,'^nmlQnidiQm)id QA 

2 2a^ 



(7.74) 



^nmi = P,xIKxJL^J P^xMN<pM^^e['efJeief^^ 

+ ^Px-x-x--(PA-Qrt'^)^/^e(^ef)ei-^e^)e(^^ef), (7.75) 

Tnmi = P,xiKxJL^Jp^xMNr'^''e['e^hie^ . (7.76) 

We shall now give some useful formulae of the previous quantities for the different 
infiationary models considered in this work. 

K-infiation 

For K-inflation model we have 

SnmZ = (2XPx)^2 I '^^ 5l{nSml) + ^^^^ SnlSml^ll | , (7.77) 

T^^i = {2XPx)--^^SniSmi. (7.78) 
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DBI-inflation 

For the DBI-inflation scenario they are given by 

1-c? 



Enmi = (2XP^)-2^^5i(„<5w), (7.79) 

where should be understood as the sound speed defined in Eq. (17.521) . 
Generalized case 

For the generalized case of subsection 17.4.31 'Enmi and Tnmi can be written as 

'AX^P,yyy (l-cL)(l-cL) 



3 Pv c^^c^ 



,2 N 



'^ad'^en 



(7.81) 



~ -i /I - 2(1 - Cg„) \ 

T'nmi = (2XPy) 2 I __2_^^-|^(5^^ _ — — 5n{m5l)l) j ,(7.82) 

and it is obvious that the DBI-infiation is a specific case of the general model 



with cl^ = 4^ = cl 



7.5.3 The three point functions of the fields 

In this subsection, we derive the three point functions of the adiabatic and entropy 
fields in the generalized case and at leading order in slow-roll and in the small 
sound speeds limit. We consider the two-field case with the adiabatic field a and 
the entropy field s. 

The perturbations are promoted to quantum operators as 

Qn(r, x) = ^ J d'kQ^ir, k)e^''■^ (7.83) 

where 

Qnir, k) = Un{T, k)a„(k) + <(r, -k)at (-k). (7.84) 
a„(k) and a|^(— k) are the annihilation and creation operator respectively, that 
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satisfy the usual commutation relations 

[an{k^),alik2)] = {27rf5^^\k^-k2)5^„,, 

[a„(ki),a^(k2)] = [at(ki),aj„(k2)] =0. (7.85) 

At leading order the solution for the mode functions is (see appendix [E] for details) 

n„(r, k) = An-^ (1 + tkCnT) 6-""""^ (7.86) 

where Cn stands for either the adiabatic or the entropy sound speeds. 
The two point correlation function is 

(0|Q„(r = 0,ki)g„(r = 0,k2)|0) = (27r)35(3)(k^ + k2)PQ„^5„„, (7.87) 
where the power spectrum is defined as 

14 |2 tt2 tt2 

and it should be evaluated at the time of horizon crossing Cjii^ki — a^H^: 

The vacuum expectation value of the three point operator in the interaction 
picture (at first order) is [1241 1141j 



(fi|gKi,ki)g^(t,k2)g„(t,k3)|fi) 

dm [Qi{t, ki)g^(t, k2)g„(t, ks), Hi{i)\ |0), (7.89) 



-I 



to 



where to is some early time during inflation when the field's vacuum fluctuation 
are deep inside the horizons, t is some time after horizon exit. \Q) is the inter- 
acting vacuum which is different from the free theory vacuum |0). If one uses 
conformal time, it is a good approximation to perform the integration from — oo 
to because r ~ —{aH)~^. Hj denotes the interaction Hamiltonian and it is 
given by Hj = —L^, where L3 is the Lagrangian obtained from the action fl7.74p . 
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At this order, the only non-zero three point functions are 

(n|g.(o,ki)Q.(o,k2)Q.(o,k3)|n) = (27r)V3)(ki + k2 + k ^2cad|A.r 1 1 

6cL(C3 + C,)^§^ - C,k\\.-, . k3 (l + + 2^ ) + 2 cychc terms 



{n\QAO, ki)a(0, k2)Q.(0, = (27r)35(3)(ki + k2 + ka)^ 



(7.90) 



AJ^U.I^ 1 1 



Cdklk, . k^ (i + + + ^Si^^^^hh'^ + ^ 

+4C3cLcL^ - 2(C, + C,)cl,klk, . k3 (l + c,.^ + 2cL^ 

(7.91) 

where K = ki + k2 + k^, K = Cadki + Cen{k2 + /cs), cychcterms means cychc 
permutations of the three wave vectors and {k2 ^ k^) denotes a term hke the 
preceding one but with k2 and k^ interchanged. The pure adiabatic three point 
function is evaluated at the moment r* at which the total wave number K exits the 
horizon, i.e., when Kcad* = a^H*- Because of the different propagation speeds, 
the adiabatic and entropy modes become classical at different times, however at 
leading order we assume that the background dependent coefficients of ( 17.9ip do 
not vary with time and so they can also be evaluated at the moment r^,. 
The different constants Cat are given by 



Ci = (2i/2g)-5^_2i4, C2 = -2{2Hh 



1 „Z 1 „z 

2,N-i '^ad _ o^o 117-2 ^en 



'^ad ^en 

Cs = {2Hh)-'^^-^, C,= {2Hh)-U^-^-^). (7.92) 

7.5.4 The three point function of the comoving curvature 
perturbation 

In this subsection, we calculate the leading order in slow-roll three point function 
of the comoving curvature perturbation in terms of three point function of the 
fields obtained in the previous subsection. 
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During the inflationary era the comoving curvature perturbation is defined as 

""-WTp/^ (7.93) 

in the uniform curvature gauge where 6q^i = 6T^. Using the constraint equation 
6q = -2Hai fTmi) and Eq. ([72D, we get 

7^ = --^. (7.94) 
^\/^ 

It is convenient to define the entropy perturbation S as 

S = --^,[^, (7.95) 



so that Vs, — Vn,, where the subscript * means that the quantity should be 
evaluated at horizon crossing. 

In this work we will ignore the possibility that the entropy perturbations 
during inflation can lead to primordial entropy perturbations that could be ob- 
servable in the CMB. But we shall consider the effect of entropy perturbations 
on the final curvature perturbation. We will follow the analysis of Wands et al. 
|180j (see also [178[ 1181] ). where it has been shown that even on large scales the 
curvature perturbation can change in time because of the presence of entropy 
perturbations. The way the entropy perturbations are converted to curvature 
perturbations is model dependent but it was shown that this model dependence 
can be parameterized by a transfer coefficient T-jis |180] as 

Tl = Tl^+ T-jisS^ = AaQa* + AsQs*, (7.96) 

with 

X=|— ^1 , A = T„J— ^./^l . (7.97) 

Using the previous expressions we can now relate the three point function of 
the curvature perturbation to the three point functions of the fields obtained in 
the previous subsection. The three point function of the curvature perturbation 
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is given by 

(7^(kl)7^(k2)7^(k3)) = Al{Q^{k^)Q,{k2)QM) 

+A^Al {{QAki)Qs{k2)Qs{ks)) + 2perms) . (7.98) 

For the DBI-inflation case the previous equation can be simphfied and the 
total momentum dependence of the three point function of the comoving curva- 
ture perturbation is the same as in single field DBI [116j . For our general model 
this is no longer the case, i.e., the different terms of the previous equation have 
different momentum dependence. Once again one can see that DBI-infiation is a 
very particular case and more importantly it provides a distinct signature that 
enables us to distinguish it from other more general models. 

7.6 Conclusion 

In this chapter, we studied the non-Gaussianity from the bispectrum in gen- 
eral multi-field inflation models with a generic kinetic term. Our model is fairly 
general including K-inflation and DBI-inflation as special cases. We derived the 
second and third order actions for the perturbations including the effect of gravity. 
The second order action is written in terms of adiabatic and entropy perturba- 
tions. It was shown that the sound speeds for these perturbations are in general 
different. In K-infiation the entropy perturbations propagate at the speed of 
light. DBI-infiation is a special case where the sound speeds for the entropy and 
adiabatic perturbations are the same. 

Then we derive the three point function in the small sound speed limit at 
leading order in slow-roll expansion. In these approximations there exists a three 
point function between adiabatic perturbations Qa and entropy perturbations 
Qs, {Qaik.i)Qs(k2)Qs(ks)), in addition to the pure adiabatic three point func- 
tion. This mixed contribution has a different momentum dependence if the sound 
speeds for the entropy and adiabatic perturbations are different. This provides 
a possibility to distinguish between multi-field models and single field models. 
Unfortunately, in multi-field DBI case, the sound speeds for the entropy and adi- 
abatic perturbations are the same, and the mixed contribution only changes the 
amplitude of the three point function. This could help to ease the constraints on 
DBI-infiation as is discussed in Ref. |116j . 



CHAPTER 7. THE BISPECTRUM IN MULTIPLE FIELD INFLATION 143 



In order to calculate the effect of entropy perturbations on the curvature 
perturbation at recombination, we need to specify a model that describes how the 
entropy perturbations are converted to the curvature perturbations. In addition, 
even during inflation, if the trajectory in field space changes non-trivially, the 
entropy perturbations can be converted to the curvature perturbation. In this 
work, we modelled this transition by the function T-jig. It would be interesting 
to study this mixing in specific string theory motivated models. 



Chapter 8 

Conclusions and Discussions 



Brane-worlds in higher-dimensional spacetimes have been proposed to explain 
long standing and crucial problems in physics, such as the hierarchy problem or 
the cosmological constant problem. Although it was shown that some of these is- 
sues were not (yet) completely solved by the brane-world proposals, these higher- 
dimensional models provide a way to tackle these difficult problems from a new 
perspective, which might suggest new solutions. 

To tackle these unsolved questions, brane-world models always introduce new 
physics in the higher-dimensional theory. These new ingredients will always in one 
way or another give new contributions to the 4D effective theories. In most cases, 
these new 4D effective degrees of freedom will have observational consequences 
on our universe. 

In this thesis, we gave a step towards calculating the 4D observational pre- 
dictions implied by several higher-dimensional theories. These predictions are 
confronted to observations of our universe, for example, the cosmic microwave 
background and table top experiments on Newton's law at short distances. It 
is extremely encouraging that we are getting strong observational constraints on 
these higher-dimensional models which aim to provide self-consistent microphys- 
ical descriptions of the early and late time universe. 

The main goal of this thesis was the study of the 4D effective theories of 
several higher-dimensional proposals and their observational consequences. In 
order to achieve it we used a powerful technique called the gradient expansion 
method. 

For example, in chapter [2] we applied this method to the 5D dilatonic two 
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brane model that reduces to the Randall and Sundrum model and to the Hof ava- 
Witten theory for particular choices of our parameters. We were able to show that 
exact solutions of the 5D theory, that describe an instability of the model, can be 
reproduced within the 4D effective theory obtained with the gradient expansion 
method. We have looked at this instability from a 4D effective theory perspective 
and we have shown that it arises because in the Einstein frame the theory consists 
of two massless scalar fields, the dilaton and the radion (a degree of freedom that 
describes the distance between the branes). We have argued against some claims 
that the 4D effective theory allows a wider class of solutions than the 5D theory 
and we have provided a way to obtain the full 5D solutions from the solutions of 
the 4D theory. 

In chapter [3], we derived the low energy effective theory in 6D supergravity 
with resolved co-dimension one branes. We showed that the 4D effective theory 
is a Brans-Dicke theory with the Brans-Dicke parameter equal to one half. We 
have calculated the 4D effective potential for the Brans-Dicke field. Then it was 
easy to see that one can tune the potentials on the branes so that this modulus 
is massless and the static 4D spacetime has a vanishing cosmological constant. 
We have obtained dynamical solutions due to the evolution of the modulus field 
and we showed that these solutions can be identified with exact 6D solutions. If 
we do not tune the potentials then, in the Einstein frame, the Brans-Dicke field 
has an exponential potential and static solutions are not allowed. Once more 
we found that these solutions can be lifted back to 6D exact solutions, showing 
that the gradient expansion method is a powerful method. Because the Brans- 
Dicke parameter is one half, this theory is immediately ruled out by solar system 
constraints on this parameter. However, we showed that a 4D potential with a 
minimum can be engineered and 4D realistic cosmology can be obtained. The 
drawback is that once the modulus is stabilized, the brane cosmological constant 
will curve the 4D spacetime in the same way as in general relativity and the 
self-tuning mechanism will not work. There are still some hopes that this model 
will provide a quantum solution for the cosmological constant problem, see [62], 
but that discussion is outside the scope of this work. 

The final higher-dimensional theory that we studied in this thesis was the 
lOD type IIB supergravity. In chapter HI we presented a systematic way to de- 
rive the 4D effective theory in a warped compactification including fluxes and 
branes. We have obtained the 4D effective potential for the volume modulus 
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but this potential does not stabilize the volume. If one wants to have a viable 
model, the volume modulus has to be stabilized. This is usually done by in- 
troducing non-perturbative effects [100] in the 4D effective theory. It would be 
desirable to obtain the 4D effective theory from the lOD theory including lOD 
non-perturbative effects. This is an area of intense research at the moment and 
it is still an open question. This problem is related to the question of whether 
there are explicit working models of slow-roll warped D-brane inflation or not. 
As discussed in the main text, the position of a mobile D3 brane in these com- 
pactifications can be identified with the infiaton. At first sight, the potential 
generated in the compactification seems to be fiat enough to produce long lasting 
slow-roll inflation. However, it was argued in [95j that volume stabilization by 
non-perturbative effects introduces potentially fatal corrections to the infiaton's 
potential that will halt inflation. 

Only recently some of these corrections have been calculated [lOH 11821 \183\ 
1184] . It turns out that in the explicit calculated examples |182[ I184j these cor- 
rections do not cancel and the infiaton's potential is always too steep to sustain 
inflation. There is an explicit model |183j where slow-roll inflation can occur but 
the model has to be very fine tuned. The problem of constructing a fully explicit 
model of brane inflation in string theory still remains challenging due to the lack 
of understanding of all of the corrections to the 4D effective potential [183] . It is 
our hope that the gradient expansion method will help because it can be extended 
by including moduli stabilization effects as an effective lOD energy-momentum 
tensor and it can also allow for a moving D3-brane coupled to gravity. Then 
one can calculate the potential for the D3-brane in a stabilized compactification 
including the backreaction of the brane. 

To study observational consequences of brane inflation, in chapters [6] and [TJ 
we took a more phenomenological approach. We started with the 4D effective de- 
scription of the movement of a D-brane given by the DBI action and we computed 
several observational predictions for our inflationary theory. 

In chapter El we computed the leading order in slow-roll trispectrum for a 
general single field model of inflation that contains the DBI-inflation model as 
a particular case. To achieve this, we had to develop third order perturbation 
theory to obtain the fourth order interaction Hamiltonian including scalar and 
second order tensor perturbations. We pointed out that in general second order 
tensor perturbations cannot be ignored in the calculation. This fact will have 
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important implications for existing results of the trispectrum for standard kinetic 
term inflation |128|, I131j . Because our interaction Hamiltonians are exact in 
the slow-roll expansion, one can also in principle study the next order slow-roll 
corrections. This is left for future work. 

As discussed in chapter [6l for single field DBI-inflation, the size of the bispec- 
trum {/nl) is of order of l/c^ and the size of the trispectrum {t^l) is of order 
of 1/c^. So if the sound speed is Cs ~ 0.1 then /tvl of order of one hundred is 
still allowed by observations and tnl will be in the detection range of the Planck 
satellite. With the present observational constraints (and with the expected im- 
provements in the near future) we are already significantly constraining models of 
inflation constructed from higher- dimensional theories, such as DBI-infiation. Re- 
cently, several works [1591 11611 11621 1185] showed that single field ultraviolet DBI- 
infiation is under severe pressure from observations. Reference |161] in particular 
showed that there exists a microphysical upper bound on the tensor to scalar 
ratio (r < 10~^) that is incompatible with a lower bound (r > 0.1(1 — n^)) that 
applies to models that generate a red spectral index (as favored by observations) 
and a large detectable non-Gaussianity. The infrared model of DBI-infiation is 
not excluded by observations |186j . 

DBI-inflation is naturally a multiple field inflationary model. So in chapter [71 
we studied the non-Gaussianity in general multi-field inflation models that have 
DBI brane inflation as a particular example. Specifically, we derived the three 
point function of the curvature perturbation for a generalized model in the small 
sound speed limit at leading order in slow-roll. At this order, we showed that in 
addition to the pure three point function of adiabatic perturbations there exists 
a mixed three point function of adiabatic and entropy perturbations. This mixed 
contribution has a different momentum dependence if the sound speeds for the 
entropy and adiabatic perturbations are different. In multi-field DBI inflation, 
because the sound speeds are equal, this mixed contribution only changes the 
amplitude of the final three point function, and the momentum dependence is 
the same as in the DBI single field case. It has been shown |116j that this 
mixed term will in fact reduce the size of the bispectrum and this could help 
to ease the observational constraints on the ultraviolet model of DBI-inflation. 
To find out how big this reduction is, one has to know the details of how the 
entropy perturbations are converted to the curvature perturbations. This point 
was not answered in this thesis and will be left for future projects. An important 
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outcome of chapter [7] is the fact that we now have a theoretical prediction for the 
bispectrum for quite general models of multi-field inflation, such as multi-field 
DBI-infiation and K-infiation. 

We shall conclude by saying that the future looks promising because we have 
been able to calculate more and more observable predictions from our fundamen- 
tal theories and observations are becoming better and more diversified. From 
the meeting of these two we are seeing that some models are being ruled out, 
some still prevail. More importantly, we are progressing in our knowledge of the 
universe. 
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Appendix A 

quantities 



In this appendix, we analyze the equations of motion present in chapter [3] that 
are of order e^/^ in the gradient expansion. We will show that all the 0{€^^'^) 
quantities vanish. 

The /i component of the 0{e^^'^) Maxwell equations reads 



a, 



(1/2) \ 



^4 py^J ^ 



and thus we have 



r 



The 0{e^^'^) evolution equation reduces to 

e-*/2 I 



Li 



(1/2) 



1 (0) (1/2) 

CI 



which can be integrated to give 



(1/2) 



(1/2) (1/2) 



(A.l) 



(A.2) 



(A.3) 



(A.4) 



The 0{e) evolution equations contain terms like F^yF^y(x h^xCiCi/ f{y). In 

to avoid the singular behaviour at the 

in the cap 



the cap regions, we thus require = 

(1/2) 

poles. Further, the regularity of Kg" at the poles imposes 
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regions. To fix the integration constants in the central bulk, we use the Maxwell 
jump conditions and Israel conditions at each brane: 



^'/^^ ^(0) 



(1/2) 



-ef/(9'^S-eA^)^^/^^ 
U 



M4 



(n - e4°^) (a'^S - eA^") 



(1/2) 



(A.5) 
(A.6) 



(1/2) (1/2) 

Combining these two equations and noting that F'^^= =Kq in each cap, we 
obtain two linear algebraic equations for the bulk values of Ci and C2 : 



(1/2) 



eM4 



ij(n-e4°')".f-®--*"" 



(A.7) 



2/±=F<: 



Therefore, = C2 = in the bulk. Now we also see that 



(1/2) 



on the branes, 



and = C(53/2). 



Appendix B 



Solving the bulk evolution 
equations 



The structure of the bulk evolution equations and boundary conditions of chap- 
ter [3] are identical for the three variables K^, K, and J (denoted here by the 
generic variable K). In this appendix, we summarize the procedure to solve these 
equations and to find the integration constants. 

All of the key evolution equations in the main text have the form of 




(B.l) 



subject to the boundary conditions 



K{yN,x) = K{ys,x) = 0, 



(B.2) 
(B.3) 



In the south cap we have the solution 



K 




(B.4) 



In the bulk the solution can be written as 



K 



e"/' [y'LpR + Cjx)] 



(B.5) 
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where the integration constant C{x) is determined by the condition (]B.3|) as 



C{x) = [-y'_Lo + {y'_ - yl)L.] R + 2./! v^e-^/^T" 



(B.6) 



In the north cap we have the solution 



K = yl^e-'/'L^R, 



(B.7) 



and the boundary condition (1B.3P requires 



[{y% - yl)L^ + ylL,] R + C= -2y^ v^e^^/^T^ 



fB.8) 



In the above we defined f± := f{y±). Substituting Eq. flB.6|) into Eq. fIB.SI) we 
obtain 



( r Ljydy) ■ R = -Y^yj^.e^'^'T. 

\Jys / i=± 



(B.9) 



Appendix C 

Gauge transformations up to 
second order 



In this appendix, we will find the change of variables that one needs to perform 
to go from the uniform curvature gauge (16.451) to the comoving gauge fl6.2ip . 
A similar result can be found in |124j . Alternative but equivalent methods to 
find the gauge transformations up to second order and beyond can be found in 
[T87lfT88lfT89lfT90] . 

In order to go from the gauge (16.451) where the field fluctuation is not zero 
to the gauge (16.21 1) where 50 = we need a change of variables that satisfy 
0(t + T(t)) + 50(t + T(t)) = 0(t). 

At first order in perturbation theory we only need to do a time reparametriza- 
tion. Let t and t be the time coordinates in the gauges ( 16.211) and ( 16.451) respec- 
tively. The time reparametrization is t = t + T. At first order 

T=4=« C = 4**. (C.I) 

At second order the time reparametrization is 



T = --:-~^^^ + -^. (C.2) 
00 200 ^ 

At this order we also need to perform a spatial reparametrization given by = 
X* + e*(x, t), where e* is of second order in the perturbations. The metric in the 
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gauge fl6.2ip becomes 

dT dT ^Ji)dT ,Ai)dT ... ^ f de. deA 

^„2^2ffT+«T» + ^,_^.(() + ^ (C.3) 

where nI^"* is the first order shift vector in the gauge (16.451) . If the vector e* obeys 
the equation 

da ' da ■ 

+ 5^ + 5^ = + (C.4) 

with /Xjj being a transverse and traceless tensor and 5hij being defined as the first 
four terms of Eq. (IC.Sp . then the gauge transformation equations are given by 

HT^ 

C = HT + ^ + (5, 
la = iijit) + ^ij. (C.5) 

To obtain the quantities (3 and /ijj it proves to be useful to decompose in 
e* = (9*e + al with dial = 0. After a few mathematical manipulations of equation 
flC.4p one can obtain 



(3 = ^ {6hi - d~^d'&5Kj) , (C.6) 

+]^5i,d-^d'dHhik + ]^d-^d,d,5hl + ]^d-^d,d,d'dHhi^ , (c.r) 

where 5hij can be written explicitly as 

11 • 

^hij = -J^diCndjCn + {diCndj^l + SjCn^i^l) + JjCnlij, (C.8) 

where ^/'i is from the uniform curvature gauge and we have used the variable Cn 
introduced before, Eq. fl6.38p . As is of second order now, the term Cnlij is of 
third order. We kept it in Eq. (1C.3P for the sake of comparison with the result 
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of [124] . For ( we have 

C = Cn + f C + -f^C + UnCn + (C.9) 



Appendix D 

Extraction of TT part of a tensor 



In this appendix, we provide a way to extract the transverse and traceless (TT) 
part of a tensor . This will be used to obtain the TT part of the source term of 
equation (I6.59p . 

Let Tij to be a given 3D symmetric tensor, as it is the case for the source of 
equation f l6.59p . Then it can be decomposed into a trace part and a traceless part 
as 

T., = |5.,+T,,. (D.l) 
The traceless part (5 degrees of freedom) can be written like 

fij = DijX + diXj + djXi + Xiv (D-2) 

with Dij = didj — l^ijd'^, d'^Xi = and d'^Xij = = XiJ where indices are 
raised by 6ij. The equation d'^Tij = \djT + ^djd^x + d'^Xj can be solved using a 
similar method as the one we used to solve the second order momentum constraint 
previously. We then find 

X = ^d-'d^F„ X, = d-'F^ - djd^d-'F,, (D.3) 

where Fi = d^Tij - \diT. And 

T 

Xij = Tij - -6ij - DijX - diXj - djXi- (D-4) 

In conclusion, given a tensor Tij, Eq. (]D.4p defines its transverse and traceless 
part. 
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Let us see how this works at the action level for the particular case of the 
tensor perturbations described in the main text. In the action (I6.58p . the source 
for is of the form 

S = [ dtd^xf^Tij. (D.5) 



with Tij being quadratic in S(f). We can see that because 'jij is transverse and 
traceless we are allowed to replace T^j in the previous action with Xij defined in 
(]D.4p . If we calculate the equations of motion by varying the resulting action we 



get as a source Xij ^ot simply Tij (see Eq. (16.591) ). ensuring that both sides 
of the equations of motion are transverse and traceless. 



Appendix E 

Equations of motion for the 
fluctuations 



Here, we derive the equations of motion for linear perturbations for the general- 
ized model introduced in I7.4.3I In terms of the field space "covariant quantities" 
|144j which are given by 



VtQ' = Q' + T 
VjVjP = Pjj 

^^KLJ = ^^KJ,L ~ r KL,J 



rf,p 



I -pi -pM -pi -pM 



(E.l) 
(E.2) 
(E.3) 
(E.4) 



{Tjj^ denotes the Christoffel symbols associated with the field space metric Gjj), 
the second order action can be expressed as 



(2) = 2 ' dtd^xa? 
1 ^ 



PyG 



P. 



y<j/j -t- r vYVlVJ 



--Py [(1 + hX)Gu - hXij] d.Q'd'Q' - MuQ'Q' + 2PyjM'VtQ' 

(X 



(E.5) 



with the effective squared mass matrix 
Mij = 



V p 

-^{Pyj^pj + Pyj^j) 



'-'ad 



2H^ 



•^i(pj 











'-'ad/ 





(E.6) 
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It is worth noting that except for the coefficients of the kinetic term and the 
gradient term, this action is the same as K-inflation case and DBI-inflation case 
which are derived in [118] and |116j . respectively. 

From now on we will derive the equations of motion for the fluctuations. For 
simplicity, let us now restrict our attention to the two field case (/ = 1,2). Then, 
the perturbations can be decomposed into = Qa^^i+Qs^i, where = e{ and 
is the unit vector orthogonal to e^. As in standard inflation, it is more convenient 
to use conformal time t = J dt/a{t) and define the canonically normalized fieldi^ 



Cad 



'Qcr ) 



-Qs 



(E.7) 



From the similar calculations with K-inflation and DBI-inflation cases ana- 
lyzed by |118j and |116j . respectively, we find the equations of motion for tv and 

Vs as 



'AY. 



z 



V a iiAvs 

a z 



where the primes denote the derivative with respect to r and 



f4 



CrP^YCad 

p 



:i + cl)P,s-cla'P,Ys 



with 



a 



a& 

CadH 



P\ 



KdX P'y 

a\l Py 



P,YsP,i 



P. 



p 



.Ys 



P^YjeU PyC, 








P,, = {V,VjP)eieiP,YC 



„2 



(E.8) 
(E.9) 



(E.IO) 
(E.ll) 
(E.12) 



(E.13) 



and R denotes the Riemann scalar curvature of the field space. 



^ Since the definitions of and Qs in this work are different from the definitions in |118[ 
I116[|117] the relations between Qcr and Do-, and Q„ and Vg are also different from the analogous 
relations in those papers. 
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If we assume that the effect of the couphng ^ can be neglected when the scales 
of interest cross the sound horizons the two degrees of freedom are decoupled and 
the system can be easily quantized. If we further assume the slow-roll approxi- 
mations, the time evolution of H, Cad, and & is small with respect to that of the 
scale factor and the relations z"/z ^ and a" /a ~ hold (see subsec- 
tion [72]T] for these approximations). The solutions of flE.Sp and (IE.9P with the 
Bunch-Davies vacuum initial conditions are thus given by 

~ I f^-ikCadT i I ^ (E.14) 

V2kCad \ kCadTj 

vsk ^ -=^e-'=^^- (l - , (E.15) 

when /i^/if^ is negligible for the entropy mode. 

Therefore, the power spectra for Q„ and Qs are obtained as 

Vq^ - , Vq^ - , (E.16) 

which are evaluated at sound horizon crossing. Here for adiabatic perturbations, 
the sound horizon is determined by Cad and for entropy perturbations, it is de- 
termined by Cen- The ratio of the power spectra for the adiabatic and entropy 
modes is thus given by VqJVq^ = Cad/cen- 



